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ADDITIONAL MEETING OF igs56 JULY 10 


in the Physics Lecture ‘Vheatre of the University of Bristol 
Professor Sir Harold Jettreys, President, in the Chair 


Che President announced the death of Walter Sydney Adams, an Associat 
f the Society, and paid a brief tribute to his memory, the Fellows standing 


fhe President announced that this Additional Meeting was the eighth to 
take place outside the rooms of the Society in London. He expressed the 
thanks of the Society to the Vice-Chancellor and Senate of the University, at 
vhose invitation the meeting was being held at Bristol, and to Professor M. H. L 
Pryce for the use of the Physics Lecture ‘Theatre Ihe thanks of the Society 
vere also extended to the Warden and staff of Wills Hall, where accommodation 
had been provided for Fellows and guests participating in the various activitie: 
irranged in connection with the visit, and to Dr W. J. Bates, who had undertaken 
much of the work involved in the local organization of the meetings and associated 
ictivities 
One hundred and thirty-four presents were announced as having been 
received since the last meeting, including 
J. de Sacrobosco: Fr. lvnctini commentaria in sphaeram loannis de Sacrobosco 
accuratissima (in duo priora capita) (presented by Sir Eric Miller) 
P. A. Moore: The Planet Venus (presented by the publishers) 
Benjamin Gompertz 1 manuscript mathematical notebook and analysis and 
notation applicable to the value of life contingencies, with the author's 
manuscript notes (presented by Mrs W. M. Archibald); and 


P.W. Merrill : Lines of the chemical elements in astronomical spectra (presented 
by the Carnegie Institution of Washington). 


Dr J. Opik then spoke on changes that had recently taken place in the 
ADH organization; papers were read by Dr D. E. Blackwell, Mr B. Eelsmor 
and Mr J. E. Geake; and an account was given by Professor H. Bondi of two 
papers by Mr C. B. Haselgrove and Mr IF. Hoyle jointly 

Other meetings and activities in connection with the Society’s visit to Bristol, 
included a symposium on Schmidt Optics, held in the morning of July g, and 
on the evening of the same day, a Public Lecture, entitled The /saac Newton 


{ 
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Telescope Proposals, by the Astronomer Royal. ‘Trips to places of interest 
included a number of excursions to points of scenic and archaeological interest 
in Somerset; a visit, kindly arranged by Mr Jack Miller, to the H. H. Wills 


‘Tobacco Factory at Bedminster; and an inspection of the work in progress at 


the University Physics Laboratory. On ‘Tuesday evening, July 10, Fellows 


were the guests of the Lord Mayor and Council of the City at a reception in the 


Council House, and the programme came to an end with a dinner later in the 


evening at Wills Hall 


MEETING OF 1956 OCTOBER 12 
Professor Sir Harold Jeffreys, President, in the Chair 


he President announced the death of Gregory Shajn, an Associate of the 


Society, in whose memory the Fellows stood for a few moments in silence. 
The election by the Council of the following Fellows was duly confirmed : 


Colin H. Barrow, 87 Park Hill Road, Bexley, Kent (proposed by L. A. Brown); 

John Robert Brace, 1228 Pacific Highway, Pymble, N.S.W., Australia 
(proposed by J. H. Finch); 

John Avery Crawford, University of California, Berkeley 4, Cal., U.S.A. 
(proposed by L. D. Henyey); 

Richard Rayman Doell, University of ‘Voronto, Canada (proposed by J. 
Jacob ) 

Richard Mead Goody, Imperial College, London, 5.W.7 (proposed by R. O. 
Redman) 

‘Thomas Heath Helling 2 Poultney Road, Radford, Coventry (proposed by 
H. Gs. Miles) 

John Heywood, Norwood ‘Technical College, Knight’s Hill, London, S.W 
(propo ed by W H Marshall); 

Geradus Franciscus Wilhelmus Mulders, c/o U.S. Office of Naval Research, 
Keysign House, 429 Oxford Street, W.1 (proposed by C. W. Allen); 

Francis R. Openshaw, 5 Jointon Road, Folkestone, Kent (proposed by H. | 
Hitchins) 

Robert Paul Randall, 23g Minehead Court, Alexandra Avenue, South Harrow, 
Viiddlesex (proposed by A. Armitage) 

Samuel Shenton, 22 London Road, Dover (proposed by H. Bondi); 

Nevil Harvey Silverton, 238 Lidgett Lane, Leeds 17, Yorkshire (proposed by 
I. G. Cowling): and 

Piyal Weerasinghe 219 Peradeniya Road, Kandy, Ceylon (proposed by 
S. S. Ayer) 





THE MOON’S OBSERVED LATITUDE FROM OCCULTATIONS 


1932 1953 
Murray 
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(Received 1956 


Summar\ 


The effects on the Moon’s latitude from occultations, when compared 
with Brown’s Tables, of (i) the change of ellipticity of the Earth, from Brown's 
value (1/294) to the International value (1/297), (11) the known error of the 
FK3 equinox, (iii) the error in Brown’s Tables, discovered by Woolard, 
are derived, and are shown to be fully confirmed by the occultation lunation 
solutions from 1932-1953 

After removing these terms, the observations are analysed in two series, 
1932-1942 and 1943-1953, for periods corresponding to the Moon’s mean 
anomaly and argument of latitude 

The existence of a term in the latitude with the period of the mean anomal) 
has been well established by other investigators, and has been attributed to 
the systematic effect of limb irregularities, varying with the libration in 
longitude : the re is no dynam al foundation for such a te rm The olutions 
for the two series confirm the character of the term, but the coefficient derived 
from the first series is only 75 per cent of that derived from the second. ‘The 
phase of the term is very close to that of the libration in longitude, and it is 
suggested that it may arise through a real ellipticity of the lunar surface 

The solutions for the two series for the period of the argument of latitude 
are nearly identical and may be interpreted as corrections to the inclination 


and longitude of the node 


1. Introduction.—\t is well known that the Muvoon’s observed latitude, a 
deduced from lunation solutions of occultations, shows large systematic fluctua 
tions with periods of the order of a year; that more than one such period is 
present is evident from inspection of the results (1). It is generally acce pted that 
these periods are produced by a beat effect, arising from the fact that the mean 


epoch of the observed occultations occurs at nearly the same age of the Moon in each 


lunation; errors with true periods of the nodical, tropical and anomalistic month 
will thus appear as errors with periods of 11°74, 12°37 and 13-94 lunations respec 
tively. It will be assumed throughout this paper that any observed period is duc 
entirely to an error with the corresponding shorter period 

We investigate first certain known causes which give rise to terms in the Moon’s 
observed latitude as deduced from occultations and compared with Brown’s 
Tables. ‘Vhe sign throughout ts of that which would appear in the residuals O-C, 


whether as tabular or observational corrections. 





1. Murray, The Moon’s observed Vol. 116 
In the usual notation let 


|, = mean longitude of the Moon 
‘=< mean longitude of the Sun, 
mean longitude of the Moon’s perigee, 
mean longitude of the Moon’s node 
mean elongation of the Moon from the Sun, 
mean anomaly of the Moon, 
mean argument of latitude, 
2y = coefhcient of sin F in latitude. 
‘To a first approximation, a correction to 2y is the same as a correction to the 
inclination of the Moon’s orbit to the ecliptic. 
We assume that the mean epoch of the occultations, within each lunation, 
occurs at the instant 5:5 days before full moon, i.e. when D= 112 
2. Figure of the Earth.—'Vhe ellipticity of the Earth adopted by Brown in 
constructing his Tables was 1/294 (2); that given by modern geodetic and gravity 
data 1s 1/297 (3). ‘This latter value is that of the International Spheroid which is 
used for computing geocentric coordinates of observers for the reduction of lunar 
observations. ‘This inconsistency in the adopted system of astronomical constants 
has been pointed out in a previous paper (4). Assuming that the true value is 
1/297, we can compute the corrections required by the Moon’s tabular ecliptic 
coordinates. 
The departure of the Earth’s gravitational field from spherical symmetry 
introduces the following corrections to the elements of the Moon’s orbit, (§), (6) : 


Ellipticit 297 1/294 
oL ; 25s ; 7°261 sin 4 
Om 20348 s 2°076 sin 44 
0 4) 93°95 sind 95°96 sin 44 


4°229 COS ¢ 4°318 cos 44 


If A, denotes the correction to a tabular quantity required in changing the ellipticity 
from 1/294 te 1/297, the only significant changes in the ecliptic coordinates ar 


\,A \,L + 0-1098(A,L — A,a@)cos/ 


! 
A,B = 2A,y sin F + 0°0895 (A, L — A,S.) cos F 
or, with the values given above, 


o”-149 sin, —o”-o1osin(/+$.)+0”-o10sin(/ 


hal (2) 
\ 


A 
B= +0171 sin( f+.) + 0":007 sin( F ). 


For the purpose of the present investigation, only the component of A,A 
perpendicular to the plane of the orbit will be required. ‘The second and third 
terms in A,A are accordingly negligible. 

3. Equinox of the FK3.—The star positions used in the reduction of occulta 
tions are taken from the Zodiacal Catalogue (Z.C.) (7) whose system is that of the 
IF Ky. Modern meridian observations have shown that the equinox of the F K3 
requires a correction \ positive equinox correction E means that the tabular 
right ascensions of the stars are too small by FE. _ It follows that the observed right 
ascension of the Moon when referred to the stars is also too small by E. 
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\ small change in the right ascension, A,x, say, may be converted into ecliptic 
coordinates by the transformation 


\.A 15A\,% cose (1 — tan 8 tane sin A) 


: (3) 

Ae ) 15\,%sine cosa P 
where € 1s the obliquity of the ecliptic It is sufficient in the present case to neglect 
the second term in A,A and to substitute + ( L)forAinA,ps 


\s In a previous paper (4, p. 95) we adopt F \ 0 0*-022, from the latest 
Washington meridian observations, whenc« 


\,A 


(4) 
\,8 ' 


$. Error in Brown's Tables Ihe investigations of Woolard (8) established 
that, in order to correct an error in Brown's Tables, the tabular ecliptic latitude must 

be increased by 
\,B = +0’ (2D 7 (5) 


Putting D 113 we obtain, 


0”-073 sin(F 0-108 cos (F \+o0"-105 sin (/ 


oO 070 co (/ ). (6) 


Combination of corrections. Combining (2), (4) and (6), and neglecting the 
mall terms in A,A, we obtain 


O *303 oO *14gsin 
\,8 
44 sin(F 0":023 cos (fF o”-112sin(/ 
076 cos(/ 


lhe occultations give corrections to the orbital latitude, B readily seen 


that, with sufficient accuracy, 
\B = AB — 2) A, (8) 
hence trom \7) 


2 ( 2 ‘+025 sin (/ 
Oo I! in (Ff 0 (/ ) (Q) 


hat should appear in the lunation solution on ount of the causes hers 
( msidered 
6. Observational material Ihe observations discussed are the lunation 
lutions tor orbital latitude 6B hose from 1932 1942 have been taken fron 
lable | of a paper by Brouwer and Watts (g) Phe solutions from 1932-40 have 
been reduced to the system of the Z.C. by the application of the corrections listed 
column ‘‘d”’ of the ‘Table For 1941 and 1942 the solutions have been taker 
unchanged lhe solutions from 1943-1953 have been taken from the annual 
discussions publi hed by McBain (10) 
Since IQ4l il] stal place u ed in the occultation reduction have been taken 
trom the Z.¢ lhe whole series has thus been reduced to the same system of 
tal places. 


Only dark limb di appearance a\ ia uded in the adopted lunation 


witiol 
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7. Analysis.—Vhe original purpose of this investigation was to test the occulta- 
tions forthe f+ $2 termsin(g). ‘The results for these terms were very satisfactory, 


as will be seen below ; it was then decided to carry out the analysis for the arguments 
land F. 


On the assumption that each lunation solution refers to the epoch when 


D= 113, the required arguments are given by the following expressions 


48°80 + 2g°107 (n 
sO oo 
210°S8 + 25°38 112) 


49°40 + 30°670 (n— 112) 


vhere nis the lunation number in Brown's series. (Lunation 112 began on 1932 
Jan. 7, and is the first lunation used in the present discussion.) 
lhe observations, 6B,, were grouped in 30 ranges of each argument, centred on 
30, etc., and the twelve mean values, 68, for each 30° range were solved 
by harmonic analysis for a first harmonic only. Each lunation solution was given 
unit we ight, and the small difference in weighting, due to the fact that the number of 
lunation solutions in each 6B was not quite the same, has been neglected. 


lhe values of 5B for the arguments F 4 are given in ‘lable I. 


PAaBLe | 


Ihe first harmonics deduced from these are 


( 


”-26s8in(/ ) "O04 CO (f ) Oo °s2 in( ) o”-1ocos(F ) (11) 


Ihe agreement between (11) and the corresponding terms in (g) is all that could be 
desired. ‘lhe expected corrections to take account of the change of ellipticity, 
the error of the F Kz equinox and the error in Brown's Tables, are satisfactorily 
confirmed 
A term in sin (Ff +§2,) could also arise from an error in the obliquity of th 
ecliptic; the agreement between (9) and (11) suggests that any such error must be 
very small. 
he error in Brown’s Tables has been corrected in the /mproved Lunai 
Ephemeris (11). ‘The ellipticity used in that work is, however, still 1/294; the 
inconsistency referred to above ts thus perpetuated If, as is generally accepted, 
the true value of the ellipticity IS 1/297, observations of the Moon’s longitude 
referred to the new ephemeris, as well as those referred to Brown’s Tables, will 
contain an erroneous term o”’-1s sin $2. his term, having a period of 18-6 
ars, will not be eliminated from annual means, and will thus affect the derived 


bE. J { 1. difference 
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sefore analysis for the arguments / and F, the observations were corrected 
for the terms in (g) with arguments / +S). ‘The whole series was then divided 
into two parts, each of 136 lunations. Series | extends from lunations 112 to 247, 
ind Series I] from lunations 248 to 383. 


Che values of 5B for arguments / for each series are given in ‘Table II 


] ABLI I] 
(1) (11) y / /(11) 


oka o°'6ob ) < o°75 


0°53 21° o’61 
0°30 
oO'll 
0°go 
150 oO Q5 0 
the harmonic analyses of these give 
Series | 36 sin /+0"-04 cos / 


Series I] »”-47 sin 1+ 0"-0g cos / 


mean o"-42 sin /+0°:06 cos / 


Ihe values of 6B for argument /, Series I and I] combined, are plotted in Fig. 1 
together with the term — 0-64 —0”-42 sin (/—8 ) 
Before analysis for argument F, the observations were corrected tor (12). 


‘Table II] gives the values of 5B for the two series for argument / 


‘Tasie IT] 
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lhe harmon analy es of these give 


Series | 0": *+0"-07 cos F 


Series I] oO" *+0”":08 cos I 

mean o”: + 0"-08 cos F. 
I he mean latitudes are , from ‘Table II], 

Peries | 0” 64 

Series I] "+63 


mean 0"-64 (14 


8. Discussion lhe existence of a term with argument / in the latitude was 
first pomted out by Brown (12). Watts (13) has shown that about half of it could 
be explained by systematic irregularities of the limb varying with the libration in 
longitude I'he main terms in the libration in longitude are 


0°66 sin 2D +6°% 2g sin /—1°+27 sin (/—2D) (15) 


or with D 


o°5+7 ‘2sin(/—7>). (16) 


‘The observed term (12) is thus almost exactly in phase with the libration in longi- 


tude lig. 1 shows that the agreement between the twelve values of 5B corres 
ponding to this argument and a pure sine term ts very striking. If the term were 


due to the fortuitous distribution of features round the limb, we should not expect 


such a smooth variation with the libration. ‘Che observations suggest that the 
term may be due to the overall figure of the Moon; an apparent ellipticity of the 
limb, varying with the libration, due to an ellipsoidal figure of the Moon such as 
that discussed by Hopmann (14), could account tor the term. 

Before discussion of (13), the first term of (g) must be removed. ‘lhe terms 


which require interpretation are thus 

0”:26 sin f'+0"-os cos F. (17) 
Formally we may identify (17) with 

2A) sin F—AS cos kf (13) 
where 2Ay and \ are corrections to the inclination and node. ‘Thus 

2Ay 0:26; AS o”'6, (19) 


‘The correction to the node is acceptably small and 1s in the same sense as that 
derived by Spencer Jones in his revision of Newcomb’s occultations (15). There 
is, however, some difficulty in ascribing the coefficient of sin F in (17) toa correction 
tothe inclination. ‘lhe value of this constant used by Brown was derived from the 
Greenwich meridian observations, and that derived by Spencer Jones trom 
occultations was in almost exact agreement with it. ‘he mean epoch of both series 
was about 1875, so that the occultations and meridian observations during the 
last century both yield a value which does not represent the modern occultations. 
\ possible cause of this discrepancy is that the modern occultations are confined to 
dark limb disappearances, whereas about 30 per cent of the occultations included 
by Newcomb and Spencer Jones were reappearances. It is very probable that, if 
limb corrections are not applied, the values of the inclination obtained from the 
two limbs will differ, as the principal term in the libration in latitude varies with 


Sith I 
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‘The values of the mean latitude from Newcomb’s occultations are given by 
pe ncey Jone 5 aS 


disa ypcarances 0 61 
PI : (20) 
reappearances O 25. 


I hat derived in the present discussion (14) agrees well with the corresponding 
phase in (20) 

g. Comparison with the observations.— Vhe full line curve in Fig. 2 is the sum of 
(g), (12), (14) and (17). ‘The plotted points are the lunation solutions, 68, 
reduced to the system of the Z.C. ‘The mean residuals from the smooth curve for 


on point are 
Series | 
Series I] 


‘The larger value for the first period is due probably to fewer observations and more 
uneven distribution within each lunation, rather than intrinsic inferiority of the 
individual observations. 

Series | is co-terminous with that analysed by Brouwer and Watts (loc. cit.). 


l hey solved for the three arguments F , F, and / by least squares, weighting 


each 6B, according to its probable error, and applying limb corrections deduced 


from Watts’ statistical discussion of Hayn’s charts (13) heir solution for the 
/ and / terms wa 


ro 35 sin (/ ) o) o2cos(/ ) o”-20sin fk 0”:22 cos F 


‘The difference between the coefhcients of sin (F ) and cos F’, and the corres- 
ponding values derived in this paper, while due partly to the differences of 
weighting, must chiefly be due to the relatively short extent in time of Series I. 
During this pe riod, regressed from o to 1 45 sin , therefore, had the same 
sign for most of the period The separation of F from F was thus not 
perfect The unsuspected term in F , introduced by the error in Brown’s 
lables, complicated the separation of these terms still further. 

In order to compare the representation of the original observations, reduced 
to the system of the Z.C., by the two solutions over their common range, the 
curve of Brouwer and Watts was obtained by subtracting their residuals v, from 
5B,+‘'d’’ (loc. cit lable 1) I‘his curve is the broken curve in Fig. 2 lhe 
agreement between the two curves indicates that no serious error has been intro- 
duced by the simplified analysis adopted in this investigation 

10. Conclusions Ihe results of this investigation show that the observations 
of the Moon’s latitude confirm 


(1) the known error in Brown's Tables, 


(ii) the generally adopted value of the Earth’s ellipticity, 


(111) the known error of the FK3 equinox, 


and that these account completely for the terms with arguments F+. 

‘The analyses for the two series are similar for argument /, and nearly identical 
for argument /' and the mean latitude. ‘lhe existence of these terms is confirmed, 
but their interpretation is doubtful owing to the possible systematic effect of limb 


irregularities or figure of the Moon 
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\ good deal of the scatter of the points from the curve in Fig. 2 ts most likely 


due to the simplifying assumptions that the lunation solutions can be treated as 


individual observations, equally spaced in time. It seems to be highly desirable 
that a rigorous analysis of individual occultations should be carried out when the 
Washington limb survey is completed. When Brown’s Tables become obsolete 
in 1960, they will have been in use for 37 years, or almost two revolutions of the 
node. For half this period the places of the occulted stars have been taken trom 
the same star catalogue. ‘lhe occultations during this period form the most 
homogeneous series of lunar observations of comparable duration ever obtained. 
\ discussion of these could be expected to yield valuable information on the 
adopted lunar constants. 

My thanks are due to Dr R. d’E, Atkinson for many discussions on the subject 
of this paper, and to Mrs F. M. McBain Sadler tor communicating the occultation 


results for 1953 in advance of publication 


Royal Greenwich Observatory, 
Herstmonceux ( astle, 
Susse x 
956 May 11. 
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MEAN DAILY AREAS AND HELIOGRAPHIC LATITUDES 
OF SUNSPOTS IN THE YEAR 1955 


Royal Greenwich Observatory 


inicated by the Astronomer Roya 


(Received 1956 September 20 


‘The following results are in continuation of those for the years 1953 and 
1954 in M.N., 115, 577, 1955. ‘They are derived from the measurement at 
Herstmonceux of photographs taken at the Royal Greenwich Observatory, at the 
Royal Observatory, Cape of Good Hope, and at the Kodaikanal Observatory 
India Plates for two missing days were supplied by the Superintendent 
Mlount Wilson Observatory, California, U.S.A., and the Director, Fraunhofer 
Institut, Freiburg, Western Germany. 


TABLE I 


Mean daily areas 


A 


Projected * Corrected for foreshortening? 
No Rotation Da , A 
“ commenced photo- Umbrae whok Faculae LU mbrae Whole Faculae 
rotation UT graphed spots spots 
1954 
Dec. 2 3 , 694 281 
Jan ; 2! 340 4058 
Feb i 197 364 
Vilar i 25 1 69 263 
Ap 2 , 115 414 
May , 28 72 402 486 


June ,*4 y 661 587 
June 2” 53 300 749 
July 27°80 2" 788 723 
Auy 2 3 4 OR 875 
Pel | pt 472 1070 
Oct 2007 gio 
Nov ae ) ra) 1905 


Dec 


Ii2t 
1139 1350 


sed in millionths of the Sun’s disk 


ed in millionths of the Sun’s hemisphere 


ase Il 


Viean dail area 


A 


Corrected for foreshortening} 


No. of dav Projected* 


A 


Photo Witho ole le 
- itl = mbrac WI Fac ulae { mbrae de ; I ac ulae 
graphed spot spot spot 


305 } 9 $9 117 6 35 


36 ' 125 15 684 of 


* Expressed in millhionths of the Sun’s dis} 
i:xpressed in millionths of the Sun’s hemisphere 
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lable I gives the mean daily areas of umbrae, whole spots and faculae for 
each synodic rotation of the Sun in the year 1955. ‘lhe numbers of the rotations 
are in continuation of Carrington’s series. ‘lhe annual means for the year appear 
in ‘Table I1, which also includes similar figures for the minimum year 1954. 

lable 111 gives for each rotation in 1955 the mean daily area of the whole 
spots (corrected for foreshortening) and mean heliographic latitude of the spotted 
areas separately for the northern and southern hemispheres of the Sun. ‘The 
mean heliographic latitude of the entire spotted area and the mean distance of 
all spots from the equator are also tabulated. 


Tas.e III 
Spots north Spots south 
of the equator of the equator Viean 
A Mean distance 
Mean Viean latitude from 
No. Rotation Viean Mean . 
helio- helio of entire equatot 
of commenced daily daily 
gg graph graphy spotted of all 
rotation fe oF area area 
latitude latitude area spots 
195475 
1355 Dex 21°63 455 29°97 $2 23°70 21°79 29°03 
1356 Jan 17°96 169 29°53 122 22°73 7°78 20-5 
1357 Feb 14°30 95 20°02 61 21°81 3°63 20°72 
1358 Mar. 13°63 44 23°80 s 24°96 16°18 23°99 
1359 Api 9°93 62 29°75 43 31°30 75 30°43 
1360 May 7°18 221 26°14 63 25°64 14°59 6°03 
1361 June 3°39 69 29°11 407 23°52 16°71 R424 
1362 June 30°59 124 29°22 103 29°88 2 38 9°52 
363 July 27°80 260 18°69 312 22°07 2-63 0°53 
364 Aug 24°03 395 27°73 139 22°75 14°66 20°44 
1365 Sept. 20°28 319 23°16 66 22°13 15°35 22 °Gd 
1 366 Oct 17°57 S29 24°75 782 4°17 1°O! 24°47 
1367 Nov. 13°87 882 24°51 523 15°06 9°79 099 
368 Dex 11°15 490 24°00 357 23°32 4°06 23°71 
"TABLE IV 
Spots north Spots south 
of the equator of the equator Viean 
A Vlean distance 
Mean Viean latitude from 
Mean Mean 
Days helio- helio- of entire equator 
, : daily daily 
Year photographed graphi graph potted of all 
area area 
latitude latitude area spots 
954 366 12 28 34 23 I 0 I°s5 116 
‘ J Old Cycle I 9°25 12 S24 7°17 K-30 
1954 95 | New Cycle if 30°27 11 27°60 1°36 25 °Q6 
1955 365 325 25°39 225 22°35 71 $14 


Corresponding annual mean values of all these quantities for 1954 and 1955 
re given in ‘Table IV. 

[ables I] and IV are in continuation of those in Monthly Notices tor the 
years 1874 to 1888, 49, 381, 1889: 1589 to 1902, 63, 465, 1903: 190T to 19g14, 
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76, 402, 1916: 1913 to 1924, 85, 1007, 1925: 1923 to 1933, 94, 870, 1934: 
1933 to 1945, 110, 501, 1950: and 1944 to 1954, 115, 577, 1955- 

The chief features of the record for 1955 are as follows: 

(1) ‘The rise in sunspot activity was the steepest since the Greenwich records 
began in 1874. ‘The mean area for the second six months was more than double 
that of the first half-year. 

(2) ‘The largest spot group crossed the Sun’s central meridian on November 
14°g, its mean latitude being 24°-0 N., and mean area 1026 millionths of the 
hemisphere. In addition to this group there were six others with mean areas 
exceeding 500 millionths. 

(3) On 45 days no sunspots were seen and on 21 days faculae were also absent ; 
ror 1954 the figures wert 237 and 180 respectively. 

(4) The ratio of mean corrected areas of faculae/sunspots was 1°44, and that 

f mean corrected areas of umbrae/whole spots 0-174. 


(5) As in 1954, the northern hemisphere was more spotted than the southern. 


lhe number and distribution, northern and southern hemispheres, of spot 


oup ot 
(a) two days’ duration or longet 
(b) one day’s duration 


vere a follow 


(a) 
Northern spots 96 
Southern spots 65 


‘Total 161 


(6) ‘Vhe following table gives the mean daily areas of sunspots and facula¢ 
corrected for foreshortening and expressed in millionths of the Sun’s hemisphere, 


for each calendar month. 


Vionth, 1955 Spots Faculac 


January 5 4 400 
lebruary 28 520 
Viarch 7 311 
April ; 443 
Via sig 
June 716 
July $31 
\ugust 57. 769 
September 1126 
October 1163 
November 1273 


December 77 1445 


Royal Greenwich Observatory, 
Herstmonceux Castle, Sussex 
1956 September 19 
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Summary) 


In an attempt to determine the vertical distribution of the velocity field 
from its centre-hmb variation, 274 line-of-sight velocities obtained in an 
earlier paper (from Fe line Doppler displacements) are analysed. ‘They are 
divided into 3 groups, whose successive positions are separated by 11 -o of 
heliographic longitude, and a relation is sought between the velocity dispersion 
displayed by a group and its position on the disk. No definitive result is 
obtained 

A positive result is, however, derived from a comparison of velocity 
measurements on strong (the two Na D) lines and on weak (one Ni and one 
T1) lines having, respectively, average Rowland intensities 25 and 2:5. ‘The 
Doppler shifts are found from microphotometer tracings at 40 positions on 
each of 4 of the spectra earlier used for the Fe line measurements At each 
of the 40 positions the four measurements are combined, thereby reducing the 
error in the velocity at any point to + 0°05(8)+0°03(7) km sec! for the 
Ni-Ti sequence and to + 0°06(0) + 0°03(8) km sec ' for the Na sequence: 

After correction for effects of measurement error, the mean velocit 
amplitude is +0°16(2) km sec”! for the Ni-T'l sequence, compared with 
+ 0'12(0) km sec for the Na 4 statistical discussion shows that the 
velocity sequences are similar, but that the difference in mean amplitude ts 
significant. From a simple model of line formation and taking account of 
redistribution of light by the apparatus function of the spectroscope, the 
weak and strong line velocities are shown to refer to approximate optical 
depths 0°42 and 0°32, respectively, as measured in the continuum. ‘The 
amplitude of the velocity field, therefore, diminishes extremely rapidly with 


increasing height in the photosphere 


Introduction.—Vhe presence of a velocity field in the equatorial region of 
the solar photosphere has been demonstrated in a previous paper (7). It was 
shown there that the field takes the form of irregular deviations from the mean 
i 


rotational velocity, having line-of-sight components of up to +0°5kmsec ' and 


extending over distances of the order of 1ro4km. ‘hat this velocity distribution 
retains its identity for a considerable period of time was shown by the marked 
stability observed over a period of an hour. However, betore the motion can be 
completely specified and the physical mechanism producing it identified, several 
more facts need to be known. Of these, the distribution of the field with depth 


in the photosphere will clearly be of paramount importance ; it forms the subject 


of the investigation described in the present paper 

Ihere are two ways in which the variation of velocity with depth can be 
determined. ‘lhe centre-limb variation of the velocity field is one and 
comparison of the velocities found from weak and from strong lines is the other 
In the first section the former method is considered; the two sub equent section 


contain a discussion of the second method. 
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Papers in the series ‘‘ Motions in the Sun at the Photospheric Level” wall 
be referred to in the text by their serial numbers. 

1. The centre-limb variation of the velocity field.—'\he material for an analysis 
of the variation of the velocity field across the solar disk is contained in Paper V1* 
(7). It consists of velocity measurements in the equatorial region, made at the 
beginning and end of an hour’s interval, along an arc of solar surface extending 
inwards from the limb over 34°-5 of longitude. ‘These velocity values, F, and F, 
respectively, are given in VI, ‘lables IV (i) and (ii) in the form of residuals from 
the mean rotational velocity; each of the two sequences contains 140 velocities. 
Corresponding entries in the tables refer to measurements made at the same disk 
positions in the two cases (VI, Section 2.3), but owing to the solar rotation the 
tabulated heliographic longitudes differ by approximately o <5. 

lhe procedure 1s first to divide the velocity residuals into groups and then to 
see if there is any relation between the velocity dispersion displayed by a group 
and its position on the disk. It was shown in VI, Section 4.1, that the averags 
linear extent of a velocity “ oscillation’’ in the field is 2°6 x 10*km, which is 
equivalent to 2-1 of longitude. In order that the mean velocity dispersion 

hould be well determined, each group must contain several oscillations. 
\ccordingly, the measurements were divided into three groups each covering 
11 -o of longitude (and containing between them 274 of the measured velocities). 
Position on the disk 1s defined by the value of (@+6,), 6 and @, being the angles 


subtended by the sub-terrestrial point and the mid-point of the range at the 


centre of the Sun and at the observer, respectively. ‘lable 1 below gives the 
values of cos(@+6,) for the three ranges. ‘The column headed L shows thx 
relevant longitude values for entry into VI ‘Tables IV (i) and (11), and D 1s the 
root-mean-square velocity deviation calculated from the m corresponding velocity 
values. It 1s useful to have some indication of the limits between which D may 
be expected to lie, if the m velocities are regarded as a sample from a hypothetical 
infinite population of velocity measurements for the group. ‘his is provided 
by the standard error, D/[2(n—1)]"* (Paper LV (6), Section 4.1), of D, given in 
the last column of ‘lable I. 

If the emergent light from the Sun can be considered to originate at a fixed 
distance ¢ in the line-of-sight, below the surface (10a), then the radial depth of 


‘TABLE I 


dD D/[2(n—1)]*!* 


cos (0 0,) 
km sec km sec 


(1) 224°5~235°3 + 0°13(1) 
(11) 224°0-234°5 
(1) 235 240°3 o'13(9) 2 | 0°01(4) 
(11) 235 245°8 


(1) 2460°6 2 
5 Oro 
(11) 246°0 0'15(5) o'o1(o) 


| 


* References to ections, table et in Paper Vi, * Large-Scale Motions in the | juat 
Region ’’, are preceded by VI 
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observation is ¢cos(#+4,). Thus the variation of the velocity dispersion with 


cos(#+4,) gives the variation with depth in the atmosphere. ‘The values of 
cos(#+4,) and D given in ‘Table | are plotted in Fig. 1, the vertical line through 
each point showing the estimated standard error of D. No regular change of 
dispersion with height is clearly revealed. A solution for a linear variation, 
D=mcos(6+6,)+k, was therefore made by the method of least squares; in so 
doing, equal weight was given to the three regions. Values of 
m= +0°07(7)+0°16(0) km sec"! and k=0-11(8) + 0°08(6) km sec 

were obtained, and the resulting line is shown in Fig. 1. However, the relative 
sizes of m and its error show that the solution is, unfortunately, entirely formal. 





‘ - ¥ T ; i. T 


ariation of velocity disper 


lbscissae are proportional to the dept) 
Ordinates are values of the oity d 
hus no information concerning the distribution of the velocity field with 

height can be derived from the existing data. ‘lhe reason for this is clear. Owing 
to the large linear extent of a velocity oscillation, each region contains only a small 
number of complete oscillations. ‘The values of the velocity dispersions ar 
therefore very much more uncertain than would appear from the number of 
measured velocity residuals in each group. Observations at a number of well 
eparated epochs are necessary if the centre-limb variation of the field is to be 
used to determine its distribution with depth. Moreover, if line-of-sight velociti 
are used in the analysis, the centre-limb variation due to a change in velocity with 
height must be distinguished from that due to any anisotropy of the field 


2. Velocity measurements on weak and strong lines 

2.1. Plan of investigation.—With the failure of the centre-limb variation of 
the field to give its vertical distribution, the possibility of a comparison between 
veak and strong line velocities was considered. ‘The first investigation of thi 
kind was made by W. S. Adams (1), who in 1906 measured the velocities displayed 
by the Hz line of hydrogen and A 4227 of Ca and compared them with measur« 
ments made on lines of Rowland intensity 1 and 2 in the wave-length region 
4195-4295 A. As, however, the measures were differential cross-wire ones mad 
on two juxtaposed lines, the shifts found for the wide strong lines were hab! 


to a considerable ubjec tive erro! 
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‘The velocities obtained in Paper VI and analysed in Section 1 above, were 
found from measurements on six Fe lines of Rowland intensity 2-6 in the 
wave-length region 5g00 A. Consequently, the same (two) plates provide ideal 
material for the new investigation, if the velocities exhibited by the Na D lines, 
Rowland intensities 20 and 30, are also measured. As there are only two solar 
lines from which to obtain the new velocities, the determination of their Doppler 
shifts with sufficient accuracy presents some difficulty. ‘The width of the lines 
excludes the possibility of micrometer measurements ; to make microphotomete 
tracings at low magnification and then to measure these seemed the best procedure. 
ven so, it is unlikely that the accuracy will equal that of the Fe line measurements. 
Hence it was decided to limit the measurement to a region where the previous 
results had shown that the largest velocity deviations were likely to occur 
Inspection of VI, big. 4, shows that the region between heliographic longitudes 
225 ‘yg and 238 -3 on the second of the two plates is the most suitable ; the velocity 
values appear in VI, ‘Table IV (i). 

Lhe selected plate, obtained at Oxtord on 1953 May s, carries four spectra 
15°5 mm in height extending inwards from the East equatorial limb on a solar 
image 85 mm in diameter; the four spectra were exposed within a few minutes 
of each other (VI, Section 1). Each of the spectra was measured at exactly the 
same disk positions (VI, Section 3.2) and it is the mean of the four velocities at 
each disk position which appears in VI ‘Table [V (ii) and is shown graphically 
in VI, Fig. 4 

‘Tas.e I] 


errestrial line Solar line 


Wave-lengt! ) Wave-lengtl Row |. factor 
Element 


\ int > (Tl 
Na(D.) SO 
Ni 
Na(D,) 
I) 


In the wave-length range covered by the water-vapour lines needed to provide 
a reliable standard of zero velocity (16) for the new measurements, one Ni and 
one ‘li line also occur The intensities of these lines are closely comparable 
with those of the Fe lines previously measured (VI, ‘Table I1), so that a comparison 
of the velocities derived from the two sets of lines will provide a most stringent 
test of the reliability of the new results, since only one water-vapour line is 
common to both series of measurements. ‘The wave-lengths recommended by 
the International Astronomical Union (2, §) and the Rowland intensities (20) 
for the four solar and five water-vapour lines are given in ‘lable II. 

2.2. Photometry.— Photometry of the spectra was performed with the Oxford 
Moll-type microphotometer (14). Light from a constant-voltage lamp is focused 
on a slit, whose height and width can be varied and which is in turn focused on 


’ 


the plate being measured. The plate is driven across the beam at a uniform rate 


and the transmitted light falls on a Moll quick vacuum thermocouple. ‘This is 
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connected to a Moll microgalvanometer, the deflection of which is recorded, 
by a spot of light, on a piece of bromide paper moving in a direction perpendicular 
to the excursion of the light spot. Simultaneously with the registration of the 
plate transmission, a coordinate system is recorded on the paper. The horizontal 
lines of the grid are impressed continuously throughout the recording. A series 
of teeth, engaging with a contact-maker at equal fractions of a zevolution of the 
screw driving the plate carriage, actuate the device providing the vertical lines. 

Measurements of line positions made earlier by Miss Adam, from micro 
photometer tracings, had shown large discrepancies between tracings made with 
the plate travelling in opposite directions past the analysing light beam. ‘This 
torward-—reverse discordance presumably arises from the inability of the galvano- 
meter to respond sufficiently quickly to sharp changes in light intensity incident 
on the thermocouple. If this is so, the effect will be much reduced if the maximum 
light intensity, corresponding to clear film transmission of the plate, is reduced. 
Filters were therefore placed between the lamp-house and the slit to give, with 
the slit height used, a zero-clear film excursion of approximately 16 mm on the 
tracing, compared with the value of about 100 mm normally used. (Irrespective 
of whether the reduction in the maximum excursion has the desired effect, it 
makes measurement of the tracings simple (Section 2.3) and comparatively 
quick.) 

Che height of the microphotometer slit was adjusted to give a projected 
image of o-1 mm on the plate, comparable with the diaphragm aperture of 
0°13 mm used in the earlier measurements (VI, Section 2.3). ‘The projected 
width of the slit exceeded that of the thermocouple, so that the effective width, 
0:03 mm, of the latter determines the dimension of the analysing beam parallel 
to the direction in which the plate moves; 1 mm on the plate is equivalent to 
1:5 A. A tracing-to-plate magnification of 6-89 (14) was used, and the plate 
was driven at a speed of approximately 3mmmin~'. ‘The wave-length range 
containing the lines to be measured (‘lable Il) therefore occupies a length of 
some 80 mm on the tracing and takes about 4 minutes to record. 




















iG. 2 hee produc tion of a typtal 


Scale The side of one tracing umt measures 3°71 mm on the ont inal 
Lines listed in Table II are marked The short hor ntal lin indicate approximately the 


positions at which measurement ere made on the D lines 


35* 
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‘Table LV (i1) in Paper VI shows that the earlier velocity measurements were 
made at 40 positions in the longitude range on the single plate selected for further 
investigation (Section 2.1). ‘l'racings were made parallel to the dispersion, in 
the forward and reverse directions, at precisely the same positions as before, on 
cach of the four spectra on the plate (making a total of 320). Between successive 
tracings on any spectrum, the plate was therefore moved by o-1 mm perpendicular 
to the dispersion (VI, Section 2.3). Press Contrast developer was used to develop 
the tracings; they were fixed in the usual way. Fig. 2 is a reproduction of a 
typical tracing. 

2.3. Measurement.—All_ measurements were made on a _ Hilger photo- 
measuring micromete! lhe tracing was clipped between two sheets of plane 
glass, | inch thick, and mounted on a small stage on top of the main stage of the 
machine. ‘This subsidiary stage can be moved over the main stage in a direction 
perpendicular to the direction of run of the main measuring screw. Cross-wires 
oriented parallel and perpendicular to the direction of the main screw were 
carried in the eyepiece. 

A tracing was first aligned with the horizontal grid lines parallel to the run 
of the measuring screw. For all lines except the D lines, the subsidiary stage 
was then moved until the tip of the line being measured fell on the horizontal 
cross-wire. ‘The positions of the line tip and of the vertical grid lines on either 

ide at this level were noted. ‘lhe tips of the D lines are broad (Fig. 2) and, 

owing to the small number of developed silver grains in the core of the line, are 
frequently irregular. ‘lo measure the position of maximum intensity would 
introduce a very large accidental error into the measurements, probably masking 
the true velocity shifts altogether. Consequently, measurements were made on 
the sloping sides of the lines immediately below the tip (the approximate positions 
are marked in Fig. 2); the mean of the two measurements was used for the 
corresponding line position in the reductions. Measurement of the neighbouring 
grid lines was made at the same level as that of the lines. ‘lhe recorded value of 
every measurement was the mean of three settings. ‘his reduces the accidental 
error of the D line measurements, which are difficult to make, and ensures that 
the positions of the grid lines, upon which the whole process of finding the line 
hifts depends, are known as precisely as possible. 

2.4. Reduction to velocities.—During the processing of the tracings the paper 
is distorted, as shown by the variation in the distance between successive grid 
lines which were impressed during the recording at equal fractions of a revolution 
of the screw. An analysis of the measured lengths of 360 tracing units (T.U.) 
gave a mean value 8-708 mm; a displacement of 0-001 units corresponds to 
(8-708 x 10 8/6:89) mm, 1.e. to 1°3 microns, on the plate. ‘The root-mean-square 
deviation of a single grid interval from the mean, as calculated from Peters’ 
formula, is + 0-041 mm, equivalent to + 0-005 ‘T.U. Since the microphotometer 
has been in frequent use for the past twenty years, it was feared that the two teeth 
operating the grid line mechanism (Section 2.2) might have worn unequally, 
causing the distances between alternate grid lines to vary systematically. 
Investigation of the 360 values for this effect gave 8707 +0°003 mm and 
8-708 + 0-003 mm for the average alternate distances. ‘The difference between 
the two values lies within their uncertainty and is therefore not significant. 

he following method was used to reduce the tracings. ‘The measured 


position of any line was converted to its fractional position between two grid 





No. 5, 1956 Motions in the Sun at the photospheric level 495 


lines, which were numbered consecutively. In this way the line positions are 
expressed on a uniform scale, since the grid lines mark equal fractions of a 
revolution of an accurate screw. It is these reduced values which are used in 
all subsequent calculations. 

The velocity corresponding to a given line position can be found most simply, 


and also independently of any asymmetry in the horizontal apparatus function 


of the spectroscope (3, 15), by reference to some standard spectrum, usually of 
the centre of the Sun. In this case, however, the mean position of the measured 
lines was found for each of the four separate spectra and the spectrum’s own mean 
used as the standard in each case. Conversion of line shifts to velocities then 
follows the method described in VI Section 2.2; the relevant velocity factors 
are given in ‘lable II. 

\t each point on each spectrum (160 positions in all), the mean velocity from 
the Ni line and the ‘Ti line (Table I1) was then calculated, and similarly the mean 
Na velocities. All four spectra were exposed within such a short time of each 
other (VI, ‘Table I) that the measurements made at exactly the same 40 disk 
positions on each spectrum (Section 2.2) have also been made at the same positions 
on the solar surface. A mean velocity for each of the 40 positions can therefore 
be calculated. ‘hese values are given in ‘lable II] for the Ni and ‘li lines and 
for the Na lines in the columns headed V(Ni+ ‘li) and V(Na) respectively. 
h is the distance on the solar image from the limb to the point at which the velocity 
has been found; it is measured along the E-W chord of the image defined by 
the spectroscops slit (VI, Section 1). 


‘Tasce III 


I (ke) V(Ni+ Ti) VCNa) 


km sec km a 


0°04 | 


-~ +--+. +S & 


~~ + + + 2 + SS SP DS S 
x >» ~ 
> pie & 


> + 
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Three of the four solar lines used for velocity determinations are essentially 
unblended, but the Na D, line has several water-vapour companions. One of 
the advantages of the use of the mean spectrum as standard is that the wave-length 
correction normally applied in stellar radial velocity work, to reduce the mean 
algebraic residual of a blended line to zero, has automatically been made. 

Also given in ‘Table III are the Fe line velocities, from Paper VI, for the 
corresponding disk positions. Like the new measurements, they are uncorrected 
for solar rotation and are expressed as residuals from the mean velocity for the 
40 positions; they differ therefore from the entries between heliographic 
longitudes 225 -g and 238-3 in VI, ‘lable IV (11). 

2.5. Reliability of the velocities.—\n the analysis (Sections 3.1 and 3.2) of 
the weak-line and strong-line velocities the uncertainties in the measured velocities 
are needed. It is interesting first to see how far the method of photometry has 
heen successful in removing the forward—reverse discordance (Section 2.2). ‘Uh 
average k—R of the 160 Ni-'li values is + 0°00(4)kmsec”', that of the Na values 
is —0°07(4)kmsec~. Although the discrepancy has not been entirely removed 
in the case of the strong lines, the value is small, corresponding to a displacement 
on the plate of only 1 micron. 

‘The method of assessing the uncertainties introduced into the velocities during 
measurement is fully explained in VI Section 2.4. ‘lable 


IV gives the results of 
the calculations for the three groups of lines, Ni 


‘Ti, Na and water-vapour, in the 
lhe column headed § contains the average of the 160 values 
of the root-mean-square error of a single line found from the deviations of the 
n individual lines in a group. 


present instance, 


In the column headed o is the corresponding 
estimate of the error of a single line from a hypothetical infinity of measurements, 
the uncertainty of this estimate being o/[2(n—1)]'*; the value, o(y), of the 
displacement on the plate to which o corresponds is also given. ‘The error to 
be attributed to the mean of m measurements is o,,=0//n; o,,/[2(m—1)]'* ts the 
uncertainty of this value. 


TasLe [V 


} 


5 | oO | a/{2(n 1)}' os og | o {2(n 1)}' | alps) 


| km sec ! | | km sec”! | km sec"! | km sec km sec 


rons 
| 


| 


| Ni-Ti o'11(3) + o-14(2) | + o'10(0) | + o-10(0) | + 0'07(1) 


Na o-11(8) 


WweV ° 


o'14(3) o0'10(5) o'10(5) | 07(4) 
0°12(4) | 0°04(4) 0°05(6) | 02(0) 
| 

Fe ° 


| w-v 0160) 


o0°16(8) 05/3) 0'06(g) | 02(2) 





| 
| to'14(1) t 0°04(5) tr 0°05(38) t o'01(3) 


In the lower part of ‘Table IV are given, for comparison, the values of 5, o, 
etc. from Paper VI, found for measurements on the photographic plate itself 
(VI, Section 2.4). ‘The values of # found in the present case are very uniform 
and somewhat smaller than the earlier ones. (This slightly greater consistency 
of measurement is, however, achieved at the expense of a factor of three and a half 
in the time taken for measurement.) ‘The close agreement of the 4 values for the 
Na and Ni-Ti groups of lines shows how satisfactorily the positions of the strong 
1) lines are found with the method of measurement adopted. Moreover, the 
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general similarity in the values of o obtained by the two entirely different methods 
of measurement (Sections 2.2 and 2.3 and VI, Section 2.1) strongly suggests 
that the limit of accuracy, set by irregularities in the silver grain distribution on 
the plate, is being reached. 
By a simple extension of the analysis in Paper VI one obtains for the root- 
mean-square error, é, of a velocity in ‘Table II], + 0°05(8) + 0°-03(7) km sec”! 
+t 0°06(0) + 0°03(8) km sec”! for a Ni-li and a Na entry respectively. (Equal 
reliance is thus to be placed on the two velocity sequences.) On the other hand, 


and 


the average root-mean-square deviation of the velocity at a point, as calculated from 
the dispersion in the four measured velocities (Section 2.4), is + 0°05(4)kmsec™! 
for the Ni-Ti sequence and + 0-06(1)km sec”! for the Na sequence. Hence no 
greater velocity differences exist than can be ascribed to the errors of measurement 
and the tabulated mean values should represent more accurately than the 
individual: measurements the velocities actually occurring on the solar surface 
(cf. VI, Section 2.2). 

The smaller value of e, + 0°-04(5) + 0°01(4) kmsec |, for a tabulated Fe velocity 
(VI, Section 4.1) mainly reflects the fact that only two solar lines contribute to 
a measured velocity in the present investigation, compared with six in the previous 


one (‘Table IV). 


3. Distribution of the velocity field with depth in the photosphere 

3.1. Comparison of the Fe and Ni Ti velocities.-big. 3 shows the Ni-'li, Na 
and Fe velocity sequences given in ‘lable II]; the dispersion in velocity (measured 
by the deviation [{V?/n]'*, also given in ‘lable III) is for the Ni-'Ti sequence 
+ 0'17(2) km sec~', and for the Na sequence + 0°13(1)kmsec”'. Now the weight 
to be accorded to differences in velocity dispersion displayed by the groups of 
weak (Ni ‘l’i) and strong (Na) lines depends largely on the extent to which the 
weak line velocities reproduce the values found in Paper V1 from the measurement 
of Fe lines of comparable intensity (Section 2.1). ‘The resemblance between the 


Fe and Ni ‘Ti velocity sequences in Fig. 3 is very marked and is reflected in the 


close agreement of the velocity dispersion, + 0°17(7)kmsec”', of the Fe sequence 
with that of the Ni-‘l'i. Proof of the identity of the two sequences is, in fact, 
hardly necessary, but it will provide a standard for the interpretation of the 
strong-line velocities. 

The coefhicient of correlation 7 (22a) between the Fe and Ni-Ti sequences 
is + 0-926 and one wishes to know if this shows that the two sequences have the 
same “‘pattern’’ or whether the correlation could arise by random sampling 
from an uncorrelated population (correlation coefficient p=0). For this purpose 
it is best to use Fisher’s z—¢ test (gb); z and ¢ are transformations of r and p. 
One finds that the probability of the observed correlation between the series 
arising by chance (1 in 10**) is utterly negligible. 

Having established the identity in pattern of the two sequences one must 
also show, in order to prove complete identity, that the patterns have the same 
scale. ‘This necessitates a comparison of the velocity dispersions of the two 
sequence lhe apparent velocity dispersion, [2V?/n|'*, calculated from the 
entries in lable II], is larger than the true dispersion owing to the effects of 
errors of measurement. ‘lhe best estimate of the true dispersion is given by (22b) 


nd? = 1 V* — ne*, (1) 
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bv C'omparison of the velocities displayed by weak and by strong line 


Lhscissae Distances from the limb along an E—~W. chord of the solar imawe 

Ordinates Line-of-sight velocitu 

Upper curve Velocities obtained from two Na lines, average Rowland intensity 25 

Viddle curve: Velocities obtained from one Ni and one Tj line, average Rowland intensity 2°5 


Velocities obtained from six Fe lines. average Rowland intensity % 7 
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where e is the error in a tabulated velocity, given in Section 2.5. Values of 


d, + 0:17(2)kmsec™! and d,= +0°16(2)kmsec ' are obtained for the Fe and 


Ni- Ti sequences respectively. 
‘The question is now, is this difference in dispersion significant or could th 
two dispersions have arisen by random sampling from the same parent popu 


lation? It may be answered by calculating the quantity 


P=(F—1)(n—2)"?/[4(1 — 9?) F]??, 


due to Pitman (12), where F = d,*/d,* and r is the coefficient of correlation between 
the two sequences of m members. With the present values of d and r and n= 4 

(Section 2.4), P=0'933. P is distributed like “ Student's” ¢ (ga) with n— 2 
degrees of freedom; for fairly large values of n the distribution is approximately 
normal with unit standard deviation about zero mean. ‘The deviation of P from 
the mean of the distribution is thus less than the standard deviation, which 1 

not significant; the two dispersions do come from the same parent population 
as one expt cts. 

lhe identity of the Fe and Nii sequences has thus been formally demon 
trated in a satisfactory manner, and the results of the comparison of strong-lins 
and weak-line velocities may be treated with considerable confidence. 

3.2. Comparison of the ke and Na velocities. —Following the procedure laid 
down in the last sub-section, the Na sequence of velocities will now be compared 
with the weak-line sequence. Since the velocities in the Fe series were derived 
from measurements on six solar lines, while those in the Ni-'li series depend 
only on two (Section 2.1), the Fe sequence represents more nearly the true velocity 
variation of the weak lines. Consequently, the k« equence will be used as the 
tandard for comparison. 

The coefhicient of correlation between the Fe and Na series | O'745 
Again one wishes to know whether this ts likely to be only a chance resemblance. 
Application of Fisher's z—¢ test shows that the probability of this happening ts 
less than 1 in 10°, which is negligible. Confirmation of the obvious similarity 
of the three velocity sequences shown by Fig. 3 is thus provided 

It has, however, already been pointed out (Section 3.1) that the amplitude 
of the Na velocity sequence, +0°131 kmsec™! (‘Table II1), ts smaller than that 
of either of the weak-line sequences Is the difference significant? ‘The true 
dispersion, obtained from equation (1) by substitution of the relevant value 


: Consequently, the value of Pitman’s P for a 


is d, + 0'12(0) km sec 
comparison of the Fe and Na sequences ts 3°384, and this is the multiple of th 
tandard deviation by which P deviates from the mean of the nearly normal 
P-distribution. Hence the probability that the Na velocities arise from a parent 
population having the same dispersion as that of the Fe and Ni ‘Ti sequence 
is only 1 in 1400 

Lhe conclusion to be drawn from the present investigation 1s therefore, that 
the velocity field displayed by strong lines is similar to that shown by the weak 
lines, but that its amplitude is reduced. (This fact incidentally explains why 
the correlation coefhcient found in this section | maller than the one in 
Section 3.1: the proportionately greater effect of errors of measurement destroy 
more of the correlation.) ‘The next problem ts to find the vertical distribution 


of the field implied by the observed change in amplitude 
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3.3. Location of the observed velocity field.—An expression for the mean 
optical depth 7,,, of the radiation /(o, #) emerging at an angle 4 (Section 1) from 
the solar surface has been given by H. H. Plaskett (13). It is 


| rdi I(r, aexp| | (1 +%,) sec ads’ | | 


r 


| d< I(r, 9) exp] } (14 m1) sectdr' | (2) 


“ 0) ~ 

where /(7,@) is the intensity at optical depth 7 in the continuum and 7, is the 
ratio of the line absorption coefficient at wave-length A to the coefficient of 
continuous absorption at the same wave-length; 7, varies from a maximum 
in the centre of the line to zero in the continuum. Since velocities were deduced 
trom displacements of the centres of intensity of the spectral lines (Section 2.3), 
a knowledge of the variation of /(7,9) and 7, for the centres of the lines would 
give the photospheric depths to which the measured velocities refer. 

In the present case, as the data are not available, it is only possible to estimate 
the mean depths of formation of the emergent light for a representative weak 
line and for a representative strong line, using some simple model of line 
formation 

lor an atmosphere in local thermodynamic equilibrium and stratified in 
plane parallel layers, the equation of transfer formulated by A. 5S. Eddington (4) 


o1ryve 


d I(r, #)exp ; (J psec dr’ | 


) 


[(1 + en, )B(r) + (1 €)y J | sec vexp| } (1+7 ysectdr' | dr, (3) 
where A(7r) is the Planckian source function, / is the average of /(7,4) over the 
solid angle 47 and « is the fraction of atoms de-excited by super-elastic collisions. 
(The equation implies that an absorbed photon is re-emitted without change in 
frequency.) With a linear approximation to the source function, B(r) = By(1 + Br), 
equation (3) ts easily evaluated (21) if 7, and € are assumed independent of optical 
depth; the integral is in fact the emergent intensity /(o,#). Limb darkening 
observations at A5goo A give B = 2-0 (17). ‘The value of 7,, corresponding to 
given values of », and « can then be found from equation (2) (13). 

\ppropriate values of « and 7, must now be assigned. ‘The simplest assumption 
is that for the strong lines « = 0 (pure coherent scattering) and that for the weak 
lines «= 1 (pure absorption). Simple asymptotic forms for 7, in the line core 
and in the wings are (18) 

AK exp | —(AA/b)*] (4a) 

and 


Kho \ 17 \A)*, (45) 


respectively; AA is the wave-length difference from the line centre, 6 is a species 
of Doppler width and 6 is a collisional constant. ‘The value of K is chosen to 
give a line of the correct equivalent width at the centre of the disk, 9@=0. For 
the weak lines, only the Doppler core (equation (4a)) is developed and K = 10, 
b =~ 2°63 x 10°*A (actually the Doppler half-width for an Fe atom at 6000 deg.) 
vive a line of some 50 mA equivalent width at the centre of the disk, roughly 
corresponding to the values (11) relevant to the weak ‘Ti, Fe and Ni lines which 
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have been measuied. Measured equivalent widths of the Na D, and D, lines 
are 636 mA and go6 mA respectively (18). The D, line is represented approxi- 
mately by values of K=5 x 10%, b= 4-09 x 10° A (the Doppler half-width for a 
Na atom at 6ooodeg.) and 6=6:49= 10 A (corresponding to a collisional 
damping width some cleven times the natural damping width), giving a line of 
650 mA equivalent width. With these data it is possible to find from equation 
(2) the optical depth of the emergent light not only at the centre of the lines, 
but also at any other point in the line profile. ‘This depth decreases from the 
value 0°58 in the continuum to a minimum in the centre of the line, o-o4 for the 
model weak line and o-ooo1 for the strong ling 

Ihe values of 7,, s0 found, however, take no account of the redistribution of 
light by the horizontal apparatus function of the spectroscope. Let /(A, 0,4) 
be the value of /(0,@) found from equation (3) for wave-length A and let 7,,(A) 


be the value of 7,, calculated from equation (2) for the same wave-length. ‘Then 


” 


the effective optical depth 7,(A) at which light of apparent wave-length A originates 


is given by 


r (A) r (A’)I(X’, 0, OPVA(A — A‘) dA I(A',0, #™)A(A— AX’) dX’, (s) 


vhere A is the horizontal apparatus function of the spectroscope; it has been 
determined for the Oxford spectroscope in the region A 5goo A by R. J. Bray (3) 

‘The measurements of velocity described in Section 2 refer to a mean disk 
position = cos-'o-4. Consequently, the values of 7, for the centres (Section 2.3) 
of the two representative lines were found for this value of 4 from equation (5) 
by numerical integration. For this purpose, the experimentally determined 
apparatus function was smoothed and made symmetrical; the model of line 
formation is too crude to justify a more detailed treatment. Values of 7, = 0°32 
and 7, = 0°42 were obtained for the strong line and for the weak line, respectively. 

While the treatment in the present sub-section is admittedly greatly over 
simplified, it does enable some general conclusions to be drawn from the 
observational material. In the first place, it has been shown that the information 
obtained from measurements on strong lines refers to a higher level in the 
photosphere than that derived from weak-line measurements (10b). Reference 
to a table of physical conditions in the photosphere, constructed by H. H. Plaskett 
(19) from mean limb darkening observations, shows that the difference im 
geometrical depth is only some 25 km. As the velocities given by the D line 
measurements are significantly smaller than the weak-line velocities (Section 3.2), 
the field must diminish extremely rapidly with height in the photosphere or be 
limited to a layer below a certain depth. It is also perhaps of interest that the 
ettective level of the observed velocity field lies above the granulation zone, the 
top of which occurs approximately at optical depth unity (19). However, some 
50 per cent of the light comes of course from layers below +,, so that the 
observations indicate at least the persistence of the field to very much deeper 
layers and possibly a considerable increase in velocity at these yreatet at pths 
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technique used in its solution. Dr M. G. Adam suggested the procedure used 
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wide lines, for which | am most grateful. ‘To Dr R. J. Bray | express my thanks 





Motions in the Sun at the photospheric level Vol. 116 


for providing the values of the smoothed apparatus function and for numerous 
discussions concerning line formation. 


University Observatory, 
Oxford 
19560 July 23 


References 
(x) 5. Adams, Papers of the Mount Wilson Observatory, 1, Part 1, 1911 
(2) D. Babcock, Trans. 1.A.U., 3, 93, 1928 
(3) J. Bray, W.N., 114, 540, 1954; Comm. Univ. Obs., Oxford, No. 45 
(4) 5S. Eddington, M.N., 89, 620, 1929. 
(5) Fowler, Trans. 1.A.l , 4, 00, 1932 
(6) B. Hart, W.N., 114, 17, 1954; Comm. Univ. Obs., Oxford, No. 43 (IV 


/ 


(7) B. Hart, W.N., 116, 38, 1956; Comm. Univ. Obs., Oxford, No. 54 (V1) 
(8) Houtgast, The Variations in the Profiles of Strong Fraunhofer Lines along a Radu 
of the Solar Disk (Vhesis), Utrecht, 1942 
ga) M. G. Kendall, The Advanced Theory of Statistics, Vol. 1, 4th edn., p. 239, Griffin 
London, 1948 
gb) M. G. Kendall, thid., p. 34 
10a) | \. Milne, Thermodynamics of the Stars, Handbuch der Astrophysik, Vol. U1, 
t. edn., p. 145, Springer, Berlin, 1930 
(rob) E. A. Milne, ihid., p. 157 
( IF. W. Mulders, Aequivalente Breedten van Fraunhofer-Laynen in het Zonnespectrun 
Thesis), Utrecht, 1934 
J. Gs. Pitman, Biometrika, 31, 9, 1939 
: H. H. Plaskett, W.N., g1, 870, 1931 
14) H. H. Plaskett, / 95, 160, 1934; Comm. Univ. Obs., Oxford, No. 4 
1s) H. H. Plaskett, MV , 12, 177, 1952; Comm. Umiv. Obs., Oxford, No. 33 


17) H. H. Plaskett, M/ , 114, 251, 1954; Comm. Univ. Obs., Oxford, No. 44 
(18) Hl. H. Plaskett, / , 115, 256, 1955; Comm. Univ. Obs., Oxford, No. 46 
19) Hl. H. Plaskett, Physical Conditions in the Solar Photosphere, Vistas in Astronomy 
Vol. I, p. 637, Pergamon Press, London and New York, 1955 
St. John et al., Papers of the Mount Wilson Observatory, 3, 192% 
d. R. Woolley and D. W. N. Stibbs, The Outer Layers of a Star, p. 134, Clarendot 
, Oxtord, 1953 


\ 
\ 
16) H. H. Plaskett, W.N., 112, 414, 1952; Comm. Umv. Obs., Oxford, No. 34 
\ 
\ 


Yule and M. G. Kendall, An Introduction to the Theory of Statistics, 14th edn 
, Griffin, London, 1950 
Yule and M. CG. Kendall, thid Pp 


326 





STAR FORMATION IN MAGNETIC DUST CLOUDS 
L. Mestel and L. Spitzer, Fr 
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Summary 


The paper deals with the problem of gravitational condensation in the 
presence of a magnetic field. It is shown that as long as the field is frozen 


into the contracting cloud the magnetic pressure sets a lower limit to the mass 


that can remain gravitationally bound: if the field is taken as 10°* gauss 


~ 


in regions of density 10 H atoms/cm’, this lower limit is 5 x10" 
However, if the bulk of the cloud is obscured from galactic starlight by dust 
grains, the plasma density within the cloud will decline rapidly, as ions and 
electrons attach themselves to the grains. When the plasma density is low 
enough the frictional coupling between plasma and neutral gas will be so small 
that the distorted magnetic field will be able to straighten itself, dragging the 
remains of the plasma with it, while the bulk of the cloud contracts across the 
held. With the magnetic energy so reduced to a small fraction of the 
zravitational energy, the cloud is able to break up into stars 


1. Introduction.—I\n recent years many workers have been led to postulat 
a strong large-scale magnetic field in the Galaxy. Interpretations of the observed 
polarization of starlight involve selective extinction by magnetically oriented, 
non-spherical dust grains (1, 2). Fermi proposed a magnetic mechanism for th 
acceleration of charged particles to cosmic-ray energies (3), while Chandrasekhar 
and Fermi explain the lateral equilibrium of a spiral arm as a balance between 
gravitation, turbulent pressure and magnetic pressure (4). ‘lhe field is assumed 
to arise through large-scale galactic mass-motions causing a small initial field to 
grow until approximate equipartition of energy exists between the magnetic and 
kinetic fields (5). 

The objections to the existence of such a field have been of two typ 


(1) doubts as to whether hydromagnetic turbulence will lead to equipartition of 


energy between any but the small-scale components of field and motion (6, 7), 
and (ii) concern at the difficulty of star formation in the presence of so much 


magnetic energy. It is with this second objection that this paper is concerned 
Some workers would argue that the evidence for a strong magnetic field 1 
fairly impressive, whereas our knowledge of the process of star formation is not 
sufficiently precise to yield a strong argument for or against the field’s existence. 
Others would regard a theory of star formation which successfully takes account 
of the field as removing one of the principal theoretical objections to the field. 
Whichever view one adopts, it is clear that the possible or probable existenc 
of the field should not be ignored when discussing the origin of stars. In thi 
paper the field 1s assumed, and the problem of star formation re-examined. 


* Received in original form 1956 
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2. The mechanical effect of a magnetic field.—Yor a non-magnetic mass in 
mechanical equilibrium, the virial theorem relates the thermal energy U, the 
gravitational energy (2, and K, the kinetic energy of mass-motions (8) 

2K +37(y—1)U+Q=0, (1) 
y being the ratio of the two principal specific heats. If the left-hand side of (1) 
is not zero but positive, the mass expands; if negative, it contracts. As we are 
interested in the best possible conditions for condensation, and in contrasting 
the non-magnetic case with the magnetic, we shall assume that turbulent and 
rotational kinetic energy in the cloud are small, so that K=o. We shall also put 

s/2, the value for a monatomic gas. ‘Then for a mass M of radius RK at 
temperature 7 the condition for contraction is approximately 

GM? 
R 


® and G being respectively the gas constant and the gravitational constant 
LZ } : v £ 


RMT, (2) 


Density variations in the cloud will introduce factors of order unity in (2). 
\lternatively, a mass M at temperature 7° must be of mean density 

p>3x10°7!74 M* (3) 
(M=M/.) for contraction to ensue. Conversely, the mass of the smallest 
phere able to hold itself together gravitationally decreases with decreasing 
temperature and increasing density. 

At normal interstellar densities and temperatures, M from (3) 1s much above 
tellar order. For example, with 7'= 100 deg. K, p must be at least 3 x 10°! /M?, 
© that for a solar mass to be gravitationally bound, the density must be som«e 
10° times that in a normal H1 cloud. It appears, then, that in order to build 
tars out of interstellar clouds the density must increase by contraction, without 
any great rise in temperature. ‘lhe energetics of this problem have been discussed 
by Hoyle (9) for a cloud of atomic hydrogen; his results are summarized in 
Section 4, where similar ideas are applied to the case of a cool dust cloud. 

Consider now a roughly spherical cloud with high magnetic energy, again 
neglecting turbulence and rotation. ‘lhe equilibrium virial theorem now has 
the form (8) 

2U ++ Q=0, (4) 
where #1 is the magnetic energy. If the average field within the cloud is H gauss, 
then 


J 


Sa 


= H?\ 
a~ 4 Re ( ) 1 HAR, (5) 


In a medium of high conductivity, the freezing of the lines of force into the matter 
leads to a rapid increase of field with decreasing radius* 

HR? = H,R,? (6) 
(suffixes referring to initial values), so that during contraction the magnetic 
energy increases at the same rate, « 1/R, as the negative gravitational energy. 
hus a magnetic field acts so as to dilute gravity by a factor independent of the 
density (in contrast to the centrifugal field in a cloud conserving its angular 


* The assumption that the galactic field is large-scale ’’ is implicit here; a “‘ small-scale ’’ 


tield could very well straighten itself during contraction, leading to a much smaller increase of field 


with density 
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momentum). ‘lhe smallest mass able to withstand the magnetic pressure is 
given by (4) with L/ =o, ice. 
H* 3 \? 
M <a | — : 
(6G)? \ 4xp 


As a conservative estimate for the actual galactic held, we take H eo gauss 
in regions where p = 1°66 x 10°**;_ the lower limit (7) is then = 5 x 107. A cloud 
with non-vanishing thermal energy will need even more mass to withstand both 
magnetic and thermal pressure. As long as the field is frozen into the material, 
the cloud cannot break up into fragments of mass smaller than the limit (7). 

It may be objected that the difficulties are exaggerated, in that increases of 
density without corresponding increases of field can arise through condensation 
along the lines of force. ‘Thus consider the simplified model of a spiral arm 
adopted by Chandrasekhar and Fermi (4): a cylinder of gas and stars with lines 
of magnetic force parallel to the axis. Once condensation along the field ha 
set in, it will proceed more rapidly than transverse contrection; only when the 
unstable portion of the cylinder has formed a rough spheroid, oblate to the 
direction of the field, will uniform contraction be possible, as only then will 
gravitation along the field be the same as the magnetically diluted force acros 
the field. In this way a cloud of high density and low magnetic energy could be 
formed, and stars could condense without magnetic impedance 


However, it is very doubtful whether this could occur. Jeans’ criterion to: 


the gravitational contraction of a length L of the cylinder i 


; (s 
L>e Je /m : 


where cis the velocity of sound. (Condition (2) is essentially the three-dimensional 
analogue of (8) for an isotropic medium.) But the condition for the lateral! 
equilibrium of a cylinder of mass M, radius R, and length L is (8): 
t (PRS) | GM 
R? L 
nd this yields as a condition on the density 
M 2K 7 
7R?L ~ 7GR** 


2RMT 


{10} 


Comparing (8) and (10), we see that the length of cylinder that ts just unstabl 
against longitudinal contraction is of the same order as KR. It is easy to sex 
from (7) that in order to build up densities high enough to allow star formatio; 
without magnetic impedance, we need condensation of a length of cylind« 
everal hundred times R, toa length of order R. If we take account of the decrea 
f R and the increase of H, due to the increase of density as the cylinder contract 
the length of cylinder required is even greater. It is quite implausible that 
ength so much in excess of the Jeans limit (8) should contract as a whole by suc! 
large factor without breaking up into sub-cylinders of length comparabl 
ith (8). 
It is clear, then, that a strong large-scale magnetic field will prevent sta 
fermation, as long as the field is frozen into the cloud. If a magnetic cloud | 
to break up into stars, some means must be found by which the matter may iy 


across the lines of force. 
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3. The coupling between matter and magnetic field.—'\Vhe matter in an ordinary 
H1 cloud is largely neutral hydrogen, both atomic and molecular, with a. small 
admixture of partially ionized heavier elements with their free electrons, and 
of dust grains. ‘The essentials of the problem are preserved if we write 
ny, =Nnumber density of hydrogen atoms and molecules, n,=number density 


of electrons, n number density of dust grains, and m,=number density of 


g 
ions, =n,/Z, where Z is the mean ionic charge. Alsc we write F,,= mean force 
on an electron due to collisions with ions, F,,,= mean force on an electron due 
to collisions with H atoms and molecules, F,,,= mean force on an ion due to 
collisions with H atoms and molecules. In considering the momentum balance 
for the gas, we may ignore collisions with the grains 


Per unit volume, the condition of momentum balance for the electron gas is 


( 


v,aH)\ 
en, (E4 ) +n,F,,+n,.F,,—Vp,=0, (11) 


here v, and p, are respectively the mean drift and the partial pressure of the 
vas, e the electronic charge and E the electric field. Because of the small 
clectronic mass, gravitation and inertia are neglected. For the tonic gas we 
have similarly 


vVAH ‘ dy, 
) nF, + nF, —Vp,+nmNXd = nm, 7 
( ¢ 


zen,(€ 


(12) 


# is the gravitational potential, m,; the ionic mass, and dv,/dt a total tim 


lerivative Remembering that the current density j is given by 
j =n,Zev;—n,ev, =n,e(v;—Y,), ( 
e can deduce from (11) and (12) the equation to the motion of the plasm 


iAH dv 
+ (nF, +2,F,4)—N(p, +p.) + nmX\b=nm 
f dt 
For the motion of the neutral gas, we have 


= dv 
(nF, + 2, F4)— Vpy + nygmyVb=nymy _* (1 


[hese equations illustrate how the magnetic force acts directly on the plasma, 
and affects the neutral gas only indirectly through the friction between plas: 
ind neutral gas. 

In treating collisions between light and heavy particles, it is a good approx: 
mation to assume that on the average the heavy particle acquires all the momentu: 
of relative drift of the light particle, which after a collision has a purely random 


relative motion. ‘Thus F,,, is given by 
NyF-y(V,)~M(Vyq —V_), (16) 


vhere o,,, is the cross-section for e-H collisions and (v,), the electronic thermal 


velocity. Similarly, as in an H1 region mj, <m,, F,, is of order 


NyOinlCy) yp MylVy —%))- (17) 


i 


Now in (14), both V(p; + p,) and nymV¢ must be small compared with the magneti 
force. By hypothesis we are dealing with a cloud in which the magnetic for 
is comparable with the total pressure gradient and gravitational force density. 
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But the plasma density in an H1 region is always much smaller than the total 
density; hence the equation for the motion of the plasma reduces to 


jaH ; ' dy, 
— +NyNjOy(lgq)rMylV py —Y;) 4 Z0,4(U,)7M(V iy V,)} = nym, at’ (18) 


rf 


Thus as the cloud contracts, the tendency of the distorted magnetic field to 
straighten itself, represented by the force term ja H/c, is countered only by the 
friction between plasma and neutral gas. 

As an example, consider a cloud of mass M=10* ©, and with temperature 
/'= 100deg. K; we suppose it initially to be in a state of mechanical equilibrium, 
30 that by (4) 
GM? 1 (H?R‘) 


RM R ———- (19) 


We again assume that the cloud has arisen through isotropic compression of a 
region of density n,=10, and with H=10", magnetic flux being conserved.* 
We then have 
(5:4 x 10!°)M —(3 x 10!7)M"8 = 2-4 » 
2x 104 
time of free-fall = 1/V Gp 
1s x10* 
2°5 y 10%4 
10% 
tl = 6 x 10" 


\s the currents maintaining the distorted field H flow in a region of lineat 
dimension R, the magnetic force jaH/c = H?/47R=7'5 x 10°. ‘The components 
of v, and v, in the meridional planes defined by the direction of the field must 
be nearly identical; for any meridional e.m.f. set up is rapidly cancelled out 
by slight non-uniformities of rotation across the magnetic field (10). ‘Thus 
16 and (vy), at 100 deg. K as 1-6 x 10°, we find that the friction 
balances the magnetic force if 


taking o;,—~10 


lv —¥, | = 1°5 x 10°/n, (21) 


If the cloud is falling freely under gravitation from its initial state of approximate 
equilibrium, v,—Ry/(Gp)=10°, and the drift of the plasma relative to the 
neutral gas is small compared with the inward velocity of the neutral ga 
provided n,;> 1°5 x 10°*. 

Current estimates of the proportion of metallic ions in H1 clouds give on 
such ion for 104 hydrogen atoms or molecules As in our cloud ny, =2 x 10%, 
we see that the normal proportion of plasma is too high by a factor 10* for 
any sensible drift to take place during the time of free-fall 

However, in normal H1 clouds the ratio of plasma to neutral hydrogen is 
maintained by a balance between ionization by galactic starlight, and decay of 
the plasma through a variety of causes, in particular through attachment of ions 
and electrons to the dust grains. If, then, the cloud is dense enough for the 

The compre t 1 not be strictl otropi cau f the increasing importance of baal 


npared with / ‘ 1 incre ww « ‘ t thi Dp mation 1 f no consequence 
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galactic starlight to be absorbed at the surface by the dust grains, the plasma may 
be expected to decay at the rate 
dn, 


dt 


Here a, is the area of cross-section of the grains, and it is assumed that a high 


NNO {V;)p- (22) 


proportion of the ions that hit a grain remain attached, or leave the grain as 
neutral atoms. ‘Thus n, decays exponentially with a time-constant 1/¢,(v,);n,. 
Because of their higher thermal velocity, the electrons will tend to adhere more 
rapidly still, but the rate of plasma decay must be determined by the slower rate 
of ion attachment: a negatively charged grain will repel further electrons and 
attractions. Other mechanisms of plasma decay, such as radiative recombination, 


must be less efficient, in that they are proportional to n,’. 


For a, we adopt (3x 107°)? and for n,/ny 2:5 10~ (11). ‘The mean 
6 


atomic weight of the ions we take as 10, so that (v,;),= 10°. ‘The time-constant 
of decay is then approximately 1-4 x 10! sec, which is shorter than the time of 
free-fall by a factor 10. ‘Thus, under the most favourable conditions, the plasma 
density can reach such a low value that the magnetic field is able to move the 
plasma through the freely-falling neutral gas. 

The forces acting on the plasma remain closely in balance: the magnetic 
field accelerates the plasma through the neutral gas until friction just balances 
the magnetic force. From (18), the time for this quasi-equilibrium to be reached 
is of the same order as the collision-time 1/n,0;,(v,), between ions and 
hydrogen atoms. In our cloud this time is about 10° seconds, which is far shorter 
than either the time of decay of the plasma, or the time taken by the plasma to 
move through the cloud. Hence in (12) and (18) we may ignore the inertia term 
as well as pressure and gravitation; the Joule term is also dropped, as it is well 
known to be negligible in the time available (7). ‘The electric field is then given 


by 
VAH nF, 
E 4 oO. (23) 
ti ne y 
Substituting from (13) and (18), and remembering that |,F,,, | is small compared 


with |n,F,,, |, we have 
v,AH 
E + oO. 


tf 
oH 
When combined with the induction equation cVaE , this shows that 
5 


the magnetic field follows closely the outward motion of the electrons, and 
hence of the plasma as a whole. ‘Thus the cloud magnetic field, which is just 
a distorted part of the local galactic field, straightens itself out as it forces the 
plasma through the lightly resisting neutral gas. 

As the curvature of the field decreases, the drift of plasma and field through 
the cloud slows up. If the decay of the plasma proceeds steadily, the field will 
drag the remains of the plasma out until the field strength has reached the assumed 
galactic value of 10~® gauss. On the other hand, if the extinction of the galactic 
starlight is incomplete, the plasma density will be stabilized at some value lower 
than the normal equilibrium level, and the outward drift of field and plasma 
will be halted when the field strength has reached a value somewhat higher than 
10~° gauss. However, the principal object will have been attained: an increase of 
density without a simultaneous increase of field strength at the rate given by (6). 
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The close balance between friction and magnetic force shows that hardly 
any of the magnetic energy released becomes kinetic energy of mass motion of 
the plasma. The energy equation of the problem follows by adding together 
the scalar products of (11) and (12) with v, and v, respectively : 


(4v*) + {nF,;.(¥,—v,) + 0Fin . (Vy —¥;) +1.F oy (Va —Y,)} 


i 


j.-E=nm, ai 


(mF iy +nF oy) Vn, (25) 
again ignoring gravitation and gas pressure. ‘lhe rate of working j.E of the 
electromagnetic forces measures the rate of decrease of the magnetic energy 
density; when j.€ is integrated over all space it transforms with the help of 


Maxwell’s equations to 
d _ 
—| | (H?/82)dr 
al aie i. space- 


The inertia term in (25) has already been shown to be negligible. The second 
term on the right is the dissipation of energy due to collisions between the different 
particles; the Joule term is again negligible, and the largest contribution comes 
from nF,,,.(¥,—v,). ‘The last term gives the work done on the neutral gas 
in decelerating its gravitational collapse. ‘Thus if |v, —v,| <|v,|, so that coupling 
between plasma and neutral gas is strong, the magnetic field energy increases 
at the expense of the gravitational energy maintaining the velocity v,,; the 
friction terms serve to dissipate only a little of the (gravitational) energy that 
would become magnetic energy if the field were strict/y bound into the neutral 
gas. In the critical case v,;=0, the magnetic energy remains stationary (apart 
from the effect of the small terms neglected); the gravitational energy that 
would have become magnetic energy in the frozen-in case is instead all dissipated. 
Finally, if |v,,—v;| >|v,,| and the plasma and field move back through the cloud, 
part of the magnetic energy 1s used up in decelerating the cloud, but most of it 
goes into heat.” 

4. The subsequent evolution of the cloud.-With the magnetic field within the 
cloud restored to the galactic level, and the bulk of the cloud contracting, the 
magnetic energy is soon a small fraction of the gravitational energy. We now 
summarize Hoyle’s discussion (10) of the fragmentation of an isothermal cloud 
into stars. 

During the initial contraction from a state of mechanical equilibrium, the 
gravitational energy released is almost all thermally dissipated by compression 
very little goes into macroscopic kinetic energy. After the cloud radius has 
decreased by a factor of about 3, the reverse is true, and further contraction of 
the cloud as a whole would convert nearly all the extra gravitational energy 
released into kinetic energy of mass motion. However, once densities have 
been so much increased without any accompanying rise in temperature, the 
cloud is unstable against the formation of sub-condensations—i.e. Jeans’ 


criterion (2) now applies to fragments of the cloud. Again, the initial condensation 


vork, Piddington (12) and Cowling (13) have recently pointed out the 
enormously increased dissipation of magnetic energy tn a lightly 1onized gas We are indebted 
to Professor Cowling for showing us his manuscript before publication, and for a helpful corre 
pondencs Our treatment differs from his somewhat, because in the case of most interest here the 
accelerations of neutral gas and plasma are not even approximately equal. However, as shown in 
the text, the mass-acceleration of the plasma remains small as the field straightens itself, so that the 


plasr icquire ery littl netic energ 


36* 





510 L. Mestel and L. Spitzer, Fr Vol. 116 


of the fragments dissipates nearly all the gravitational energy released. In 
Hoyle’s semi-quantitative treatment a cloud that has shrunk from R, to R,/3 
will then break up into five approximately equal masses, each of radius Rp/5 


and density 259. As long as the isothermal assumption is justified, the process 


is repeated for each fragment; it stops when the fragments become so opaque 
that in order to drive out the energy generated in free-fall, the temperature 
gradient becomes high enough to give hydrostatic support against gravity. 
‘The fragment then contracts as a whole in the Kelvin-Helmholtz time-scale, 
and no further break-up occurs. 

In his paper Hoyle treated the case of a mass of atomic hydrogen condensing 
and fragmenting during the early life of the Galaxy, with no stellar radiation 
field, and no dust grains or molecules. ‘The radiative properties of hydrogen 
keep the temperature at ~104 deg. K, so that the isothermal condition is satisfied. 
We now apply similar ideas to a dust cloud containing molecules. ‘The first 
problem is to estimate the temperature variation within the contracting cloud. 
The surface of the cloud is maintained at ~100deg. K by a balance between 
absorption of galactic radiation and radiation loss from the grains and molecules 
(14). During compression the gas is heated, but much of its energy is transferred 
to the grains by inelastic collisions. We first show that the cloud is optically 
thick in the thermal radiation from the grains and the molecular line emission. 
If Amax 1s the wave-length of maximum int nsity emitted by the grains, then, by 
Wien’s displacement law, Amax7',=0:29, 7, being the internal temperature of 
the grains. With the grain radius b= 3 x 10~°, the parameter «= 27b/A <1, and 
so the cross-sections for absorption plus scattering of radiation by the grains 
are (11) 

| 2b*(a*/3) approximately for dielectric grains | 
A(A)= 4 ? (26) 


| wh?(2-8x) approximately for metallic grains 


We shall assume that a fair proportion of the grains are metallic, so that by (26) 
much the largest contribution to the cloud opacity comes from the metallic 
grains. ‘lhe optical depth of the cloud in the grain radiation is of the order 


R, (xb? 8 ( = 1 (3°) me 27 (= (27 
(_( 70%) X 2° om, Nyty R) ae 10 y R) <7) 


where R, is the initial cloud radius. Similarly, the line-emission from the 
H, molecules is of wave-length short enough to be absorbed by the grains before 
escaping (15). We therefore expect a Planck distribution of t mperature Be to 
be built up within the cloud, in a time of orde1 


Energy/cm® of temperature radiation of 


Rate of emission by grains/cm* 


(The use of Kirchhoff’s law in the denominator of (28) is justified, since A(A) 
for the metallic grains is largely true absorption rather than scattering.) ‘This 
time is approximately 


oT! 1*3 x —(> 3 
n_c(2°82b*)\(27b/Amax)o7,* © f¥ R, 


y 
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which is far less than the time of free-fall. ‘The compressional energy of the 
cloud will build up a temperature gradient sufficiently steep to drive out this 
energy. Most of the opacity arises from the grains, the molecular line absorption 
merely making kinks in the opacity function. 
In order that the cloud should fragment, this temperature gradient must be 
mall compared with that yielding hydrostatic equilibrium; in other words, the 
thermal energy must increase much more slowly than the negative gravitational 
energy. We shall assume that the cloud stays nearly isothermal, compute the 
compressional heating and the resultant temperature gradient, and so check up 
on the isothermal assumption. At gas temperature 7’, the rate of compressional 
heating during free-fall is 
d/(M soo -  fR\ 32 
ps ) RMTE ~=yRMTy/ (Gp) =2°4 10" T'( ) (30) 
dt \ p Ro, 
on substituting from (20). When 7 and 7), differ, the rate of heating of th 
grains by inelastic collisions is given by Spitzer’s formula (15), which for out 


p 


problem gives 

(R,\3 
31 Co) T13’ 
10 (R) T'2(T — T,) 


I'he rate of loss down a temperature gradient d7’,/dR is of order 


4 a (Amax7')\ T,2|dT,|(R\°_ sgn (AT, | ( R\P 
+7 Ry (3 ser 2b ) T,,| dR R.) Ep ers dR (,) 
(22) 


Ihe equating of (30) and (31) determines the difference between 7° and Ts 
while from equating (30) and (32) we find the excess of internal 7 over the surface 
value, necessary to drive out the compressional heat. 

With R= R, and an assumed temperature 7, = 100 deg. K, the same as the 
surface value, 7’ is found to be ~105 deg. K and |d7,/dR| ~10-” with |d7/dR | 
of the same order. ‘This is much smaller than GM,/RR,?’=3 x 10", the 
temperature gradient required for hydrostatic equilibrium; this is also the cas¢ 
when R has reached R,/3, when according to Hoyle the first fragmentation takes 
place. In fact, the opacity of the grains is so low that the fragmentation process 
stops only when the cloud has broken up into masses of order ©/2. ‘The mass, 
radius and density of the mth fragments are respectively M,/k", R,/k" and pyk?" 
with k=5 (g); hence with the temperature gradient just equal to the hydrostatic 
value, the radiative flux through the fragments is 


GM, \3 
B°5 » ro ( ) k-3" =~ 10% R- (33) 


RR,) R, 


, 18 quite negligible). On the other hand, the 
rate of compressional heating during free-fall is 


u(= T V(Gpo)h (oo 
7p = y i x Co ° 
3 RM ) \ Po) 4 10 RR, ) (34) 


Fragmentation ceases when (33) and (34) are approximately equal; they yield 


(the difference between 7' and 7 


> 


k"=2 x 10%, so that the proto-stars that contract to the main-sequence are = 
in mass 

However, the successive break-up of the cloud may very well be halted before 
uch small masses are reached. We began by considering a cloud with a large 
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magnetic field but with negligible centrifugal force, and showed how the field 
may be eliminated. But rapid contraction of the cloud and its fragments, with 
conservation of angular momentum, must send up the centrifugal force. It needs 
only a moderate amount of angular momentum in the initial cloud for rotation 
to halt contraction perpendicular to the axis and lead to disk formation, long 
before the increase in opacity has stopped fragmentation. Thus the masses of 
the disks will be determined not by the equating of (33) and (34) but by the 
amount of vorticity in the cloud. 

The evolution of such a disk will depend on the means for gradual re- 
distribution of angular momentum. If the obscuration of galactic starlight is 
incomplete, coupling between magnetic field and matter may be restored, because 
of the much greater densities. ‘The strong interaction between rotation and 
magnetic field is well known. Either the bulk of the angular momentum will 
be transferred to the Galaxy, so that the disk as a whole gradually contracts into 
a main sequence star, probably with high rotation; or the cloud will break up 
into two or more masses, with most of its angular momentum going into the 
orbital motion of the proto-stars. Further discussion is beyond the scope of 
this paper. 


cs. The evolution of a cloud with a frozen-in magnetic field.—In Section 2 
s 

+ 

i 


it was shown that the minimum time of plasma decay, and hence of loss o 
magnetic energy, is shorter than the time of free-fall, but not by a very large 
factor. ‘The magnetic field does not much increase the time of free-fall; hence 
if the mechanism of attachment of ions to the grains is somewhat less efficient 
than supposed, coupling between plasma and neutral gas will remain strong, 
and the collapsing cloud will drag the field with it. However, in contrast to 
the case of a cloud with low magnetic energy, even if the contraction is isothermal 
the cloud cannot break up into smaller masses— the magnetic energy sets a lowe 
limit to the mass that can remain gravitationally bound. ‘Thus the cloud as 
a whole continues to contract, with more and more gravitational energy being 
converted into macroscopic kinetic energy after the first stages of contraction 
very little of the gravitational energy released is thermally dissipated, because 
of the impossibility of sub-condensations being formed. 

It is clear, then, that if the magnetic field is to be lost from the cloud by th« 
mechanism of this paper, something must prevent the cloud from collapsing so 
quickly, One tempting suggestion is turbulent pressure, maintained by the 
release of gravitational energy; if strong enough, it could regulate the time of 
contraction to its decay-time. However, it is doubtful if in fact the time of 
contraction will be increased in this way. Certainly, if a// the gravitational 
energy released but not thermally dissipated were to become turbulent kinetic 
energy, the virial theorem shows that the pressure would halt the contraction 
at a radius ~R,/3. But the time of decay of the turbulence is given by L/a, 
where L is the scale and a@ the mean velocity of the motion containing most of 
the energy (16). At radius R, the energy that has been released by contraction 
is somewhat less than 7(GM?*/R), where 7 is the constant factor by which the 
magnetic field reduces gravity. Hence i= /(7GM/R)=Ry/(Gpn), and so the 
time of decay is of the same order as the time of free-fall 1/\/(nGp).* In any 
case, the assumption that all the energy of the mean flow becomes turbulent 


* The suggestion that the decay-tume of the turbulence might not exceed the time of free-fall! 
vas made to one of us by Dr. F. D. Kahn 
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is probably over-generous. Experiments with wind-tunnels indicate that only 


a few per cent goes into turbulence (16); the decay-time of the turbulence is 
then longer as @ is smaller, but the turbulent energy is then too small to exert 
an important pressure. 

Thus we conclude that turbulence could be important only if nearly all the 
energy available becomes turbulent, and also if the presence of a strong magnetic 
field reduces considerably the rate at which energy cascades down to smaller 
and smaller eddy-sizes, to be dissipated by viscosity. Even so, we need to take 
account of the extra loss of energy through compressional heating and radiation ; 
this sink will be important when the kinetic energy is comparable with the thermal 
energy. 

If turbulence is negligible, there remain two other possible ways in which the 
cloud can lose its magnetic energy: 

(i) ‘The angular momentum present leads to disk formation, and the subsequent 
evolution of the disk is slow enough for the field and plasma to diffuse outwards, 
in spite of the increased densities. ‘lhe objection to this is that the strong frozen-in 
magnetic field will probably remove angular momentum too rapidly; the time 
of travel of a hydromagnetic wave across the cloud is of the same order as the 
time of free-fall, and so it is not obvious that a rotating disk will form. 

(ii) ‘The cloud as a whole falls freely until the opacity becomes too high for 
the compressional heat generated during free-fall to be radiated away. ‘The 
cloud then contracts in mechanical quasi-equilibrium, the rate of contraction 
being determined by the energy leak down the temperature gradient. Because 
of the high molecular opacity in the infra-red, this time of contraction may be 
very long, giving the magnetic field time to diffuse outwards. As the temperature 
in the cloud increases the molecules will ultimately dissociate, and the opacity 
decrease. With the magnetic energy now negligible, the cloud can fragment 
and evolve either according to Hoyle’s picture, or through formation of rotating 


disks. 


One of us (L. M.) wishes to record with thanks the generous financial aid of 
the Commonwealth Fund of New York, and the warm hospitality of the Princeton 
University Observatory. 


Faculty of Mathematics, Princeton University Observatory, 
Arts School, Princeton, N.7., 
Bene’t Street, U.S.A. 
Cambridge 
1956 July 26. 


Note added in proof.—One of us (L. 5.) has pointed out that dissociative 
recombination of electrons with molecular ions has a high cross-section, and is a 
possible alternative to grain attachment as an efficient process of plasma decay 
in an obscured cloud. 
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\ MATHEMATICAL DISCUSSION OF ‘THE PROBLEM 
OF STELLAR EVOLUTION, WITH REFERENCE ‘TO 
tHE USE OF AN AUTOMATIC DIGITAL COMPUTER 


C’. B. Haselgrove and F. Hoyle 
(Received 1956 May 10 
Summar) 


The partial differential equations describing the evolution of a star have 
been reduced to a form suitable for numerical integration by an electron 
computer. The integration has been programmed for the electronic 
computer EDSAC I and the results will be discussed in later papers 


1. Preliminary remarks.—Investigations of stellar structure have hitherto 
been concerned with a determination of the density and temperature within 
a star of prescribed composition. ‘lhe problem in such a case is governed by 
a set of non-linear ordinary differential equations. 

‘This type of treatment cannot be used to investigate stellar evolution because 
in this problem the composition cannot be prescribed except in the initial stages 
before nuclear processes become important. Subsequent changes ot composition, 
arising from the nuclear processes, must be worked out as a part of the problem. 
‘This alters the mathematical structure of the equations giving a set of non-linear 
partial differential equations in two variables, one variable being the time and 
the other describing distribution within the star. ‘lhe distance from the centre 
of the star, for instance, may be taken as the latter variable, although this is not 
the only possibility. 

‘This increased mathematical complexity would be extremely severe were it 
not that several simplifying approximations can be introduced without any 
sensible loss of accuracy. ‘hese will now be mentioned. 

2. Simplifying approximations.—(i) ‘Vhe star possesses spherical symmetry. 
‘This is an approximation because stars do not possess spherical symmetry, for 
example because of rotation, or because magnetic fields make contributions to 
the pressure. It follows that magnetic and rotary effects are being neglected. 

(ii) ‘The equation of hydrostatic equilibrium is satisfied. ‘That is to say, the 
dynamical equations of fluid motion reduce to the hydrostatic equation. ‘This 
is an excellent approximation except in cases of violent oscillation or of explosion 

as in supernovae. 


(iii) There is negligible stirring of material in the interior of a star, except 
in regions that are convectively unstable where the material may be taken to be 
well mixed (Mestel, 1953). 


(iv) In the following section the main variables of the problem are defined. 
Subsidiary variables have to be determined when the main variables are specified. 
Notable instances of subsidiary variables are the density, opacity, and the rate 
of energy generation. In all these examples the accurate determination of the 
subsidiary variables necessitates long and somewhat intricate calculations. With 
a digital computer of very large storage capacity and high speed, such calculations 
could undoubtedly be performed. But with a machine of moderate storage and 
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speed, approximations in the determination of the subsidiary variables have 
perforce to be accepted. This has been the case with all work performed at 
Cambridge with the aid of EDSAC 1. ‘The relevant approximations will be 
given later in this paper when the construction of the subsidiary variables comes 
to be considered. 
3. The main variables.—'Vhe main variables in the problem are as follows 
t, the time, 
VM, the mass inside a spherical surface within the star, 
R, the radius of this spherical surface, 
QO, the energy flux across this surface, 
P, the total pressure at radius R, 
7’, the temperature at radius R, 
X, the fraction by mass at radius R of the stellar material consisting of 
the particular element undergoing nuclear transformation. 
lhe following subsidiary variables, defined at radius R, are required 
p, the density, 
K, the opacity, 
the rate of thermonuclear energy generation per unit mass, 
the adiabatic exponent I’, of Chandrasekhar (1938), 
the potential energy per unit mass. 


‘There are two independent variables, t, and a second variable, usually taken 


as R, but chosen in our work as M. ‘The other variables, including X, are 
therefore considered as functions of M,t. It may be noted that the functional 
dependence of X on M,t is unaltered by an expansion or contraction of the 
material of the star. ‘This condition is more difficult to take account of if R is 
chosen instead of M as the second independent variable. 

\ further point concerning X needs emphasis. It seems to be generally the 
case that for given M, ¢ there is never more than one element that makes appreciable 
contribution to the energy generation by nuclear processes. ‘Thus if hydrogen 
is present no other element makes appreciable contribution. If helium ts present, 
but not hydrogen, no other element but helium makes appreciable contribution 

and so on. ‘his explains why in any particular sample of stellar material only 
one quantity X is needed to specify the concentration of the energy-producing 
clement, whatever this may be. 

As t changes, the dependence of VY on M varies on account of the nuclear 
reactions. For a particular M, the value of X decreases with increasing ¢. 
Now X cannot decline below zero. A decline to zero implies that at the values 
of M,t in question the element subject to nuclear transformation becomes 
exhausted. But this exhaustion does not demand an entire cessation of nuclear 
activity, for as ¢ continues to increase, nuclear reactions involving some other 
different element may subsequently become important. In such a case X cannot 
stay at zero value. 

To deal with this possibility we introduce a parameter f§ that takes positive 
integral values. ‘his parameter identifies the particular element subject to 
nuclear transformation. ‘lhe statement 8=1 is taken to mean that the main 
element subject to transformation is hydrogen, while B= 2 is taken to mean 
that the main element under transformation is helium, and so on for other values 
of 8 and other elements. Accordingly for any M, t we not only have X, giving 


the concentration of nuclear fuel, but also 8 specifying the identity of the fuel. 
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Should X fall to zero, the value of 8 must be changed and Y must be reset to 
take account of the concentration of a new fuel. 

A simple example will make this clearer. Suppose we start with pure hydrogen. 
Then initially X=1, B=1. As ¢ increases X decreases. Suppose X declines to 
zero, implying that all the hydrogen has been changed to helium. Then f 1s 


changed to 2 and _X is reset to 1. 

Not only is X a function of M,t, but 8 is in general a function of these 
variables. ‘This means that 8 may have different values in different parts of a 
star. 

Evidently the mathematical formula to be used for the generation of energy 
by nuclear reactions must be changed when f changes. In programming the 
calculation for a digital computer, f serves to indicate to the machine which of 
a number of routines is to be used for the calculation of the energy generation 
and also of other quantities, such as the mean molecular weight and the opacity. 

4. Procedure for determining p when T, P are specified.—\t will be noticed 
that the density p does not appear in our list of main variables. ‘The density 
can be determined, however, when 7, P are specified. ‘The following procedure 

(i) takes account of radiation pressure , 

(ii) gives the correct p when the pressure due to the material arises over- 
whelmingly from either (a) the thermal motions within the material, or 
(bh) non-relativistic degeneracy of the electrons, or (c) relativistic degeneracy of 
the electrons 

(iii) gives log,,p in the transition region between (a) and (4) to a better 
accuracy than +o-°o1; and 

(iv) gives log,,p to a comparable accuracy in the transition region between (4) 
and (c). 

Define P’ = P-}haT"‘, (1) 
a being Stefan’s constant. Define p,, py, ps as follows 

pP’ 
Py RT 
where ® is the yas constant, and 


pu re 4 | = 2, ; (3) 


P’\ 3/5 P’\34 
Pe BK, K, ’ Ps = F. x) ’ (4) 


where K,, A, are the well-known constants for non-relativistic and relativistic 
degeneracy (Chandrasekhar, 1938), and 
2 
he . B=; L 
/ 1+ ad, 


‘The formulae adopted for » in the case 8 = 1, 2, and the formula for y, for B= 1, 


~~ 9» 


are not exact, but give good approximations when, as seems always to be the 

case, the material is initially composed almost wholly of hydrogen and helium 
The required value of the density is then determined by 

p= (py? + (2? + pg”) ">"? (6) 

Still more accurate procedures than this can no doubt be obtained at the 

expense of more complicated formulae and of a longer calculation. The present 

procedure was felt to be a suitable compromise between accuracy and storage 

requirements in the case of the machine at our disposal. The values in c.g.s. units 
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of the constants required in this routine are log,,a=15°8779, log, = 7°9161, 
logy Ky = 12°996, log, K, = 15-080. 

5. General procedure for determining the evolution of a star when X is a slowly 
varying function of M.—The composition of an evolving star at a time ¢ is 
determined in part by the original composition and in part by the evolutionary 
stages that preceded ¢. For the moment, suppose that the composition at time ¢ 
is known, together with the course of the evolution preceding t. ‘Then the 
structure of the star at time ¢ can be determined by the method that will be set 
out in detail in the next section. Suppose this is done. Then X(M,t+8t), 
8(M,t+6t) can be determined, 5t being a small time increment. ‘The crucial 
point in this determination is that the structure worked out for time ¢ can be 


taken as a good approximation over the whole time increment, provided 8¢ is 


small enough. 
can be regarded as known and hence the change of X (and any change of f) can 
be computed. 

Equipped now with X(M,t+6t), B(M,t+6t) we can work out the new 
structure for time ¢+4t. ‘Then a further time step can be made, and so on. 
Accordingly the evolution can be followed provided X(M,t), B(M,t) are known 
at any one time. ‘This indeed we are entitled to suppose, since clearly the 
composition of the star at the moment of its origin must be givei:. 

It may be added that, although this step by step procedure can be made as 
accurate as one pleases by taking 4¢ small enough, this is not the best numerical 
way of achieving a given degree of accuracy. 

6. Equations for determining the structure of a star at time t.—-\n accordance 
with what has just been said X(M,t), 8(M,t) will be regarded as known. 


Differential equations 

dR I 

dM ~ 4nR%p’ 

dP GM , 

dM 47k! 

dO de 

dM ” Ot’ 

dT 3KO 

dM 647° acpT*R' 

aT’ 

dM 
_T=-1TdP. 
“Tt PdM ’ 

dT l~—1 T dP 

dM TT P dM’ 


Constants appearing in these equations 


provided the value of given by (10) is 


otherwise 


log,» G = 8-8241, log, c= 10°4768, using c.g.s. units. 
Subsidiary variables appearing in these equations 
Routines have to be given for determining the subsidiary variables p, U, 
de/ot, K, T when the main variables and A(M,t) are given. Of these the case 
of p has already been discussed. 
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(1) Procedure for obtaining ’.—The following formulae* are considered to 
give an adequate accuracy for cases B= 1, 2. 
B=1, U=[—29°05]pX*7*+[— 112-1]pXxgyT™, (12) 


where x, = fraction by mass of stellar material in the form of carbon or nitrogen. 


"3 _ x 4 
10°35 ai exp 4 oat. ot 34 ll 


T (13) 
With pin gcm~* and 7' in deg. K these formulae give the rate of energy generation 
in ergs per gram per sec. ‘lhe formula for 8=1 includes the energy generated 
both by the “proton chain’’ and by the ‘“carbon—nitrogen cycle”. Over the 
range from 7'=10%deg. K to T=5 x 10°deg. K the expression for U in (12) 
does not differ by more than 0-1 in log,, U from the accurate but more complicated 
formulae given by Borman-Crespin, Fowler and Humblet (1954). ‘The formula 
for U in (13) is taken from Hoyle (1954). If temperatures exceeding 5 x 10° deg. K 
are encountered, nuclear fuels other than hydrogen (8=1) and helium (f= 2) 
must be considered. Not enough is known at present of the energy levels of 
medium light nuclei for precise formulae to be given, but some account must 
nevertheless be taken at very high temperatures of the conversion of oxygen, 
neon, carbon, magnesium, etc. into heavier nuclei (Hoyle 1954), even though 
this has to be done in a very approximate way. Approximate formulae for 8 >2 
can be constructed as follows: 

U = Ag pX*T"®, B>2. (14) 
The index mg must be large, and can be set tentatively at some assigned value, 
30 for example. ‘Then for given £ the coefficient Ag is determined by the 
condition that U must become of order unity at a value of 7 determined by 
considerations of nuclear physics. ‘Thus for neon as a nuclear fuel the value 
of 7 in question is about 8x 10%deg. K, for oxygen the value is about 
15 x 10°deg. K, for magnesium and silicon about 2 10%deg. K. For T close 
to 5 x 10°deg. K the value of U becomes of order —1, corresponding to the 
break-down of nuclei into helium (Hoyle, 1954). ‘This is probably the last 
nuclear process that needs consideration since it seems likely that a major 
instability must arise in a star once U becomes large and negative. 

(ii) Procedure for obtaining de/ot.—The derivative 0¢/0t appears in the equations 
in order to take account of the slow changes of gravitational potential energy 
and of the heat content of stellar material that takes place during stellar evolution. 
The quantity de occurring in any small change is given by 


P’ 
de 3 52 3 d (F ) ( I 5) 


) 


provided relativistic degeneracy of the electrons is not important. 

The following procedure enables de/dt to be estimated. Denote the 
evolutionary steps by 1, 2, ..., the number 1 referring to the star when it first 
forms and reaches approximate thermal and gravitational equilibrium, and the 
remaining numbers being taken in the order of evolution. ‘Then we write 


(=) Rae - pega 
ot n t, —~ bas , 


(5 
- =O 
ot), 


* Throughout this paper we use the notation {x] for the quantity 10% = antilog,»x 
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The adoption of (16) requires d¢/dt to be omitted for the initial state. This 
is an excellent approximation since de/ot has only a very slight effect on the 
structure of the star during the early stages of the evolution. Therefore in 
determining (d¢/0t),, €,(M) is a known quantity. The time difference t,—t, is 
also known, since it is the time difference that specifies the evolutionary step. 
Hence d¢/dt is determined in terms of M, P, T for the stage 2, and so on for 
subsequent steps. 

The inclusion of d¢/df emphasizes that the basic equations are of the partial 
type. Impressions to the contrary sometimes arise from the circumstance that 
de/At is usually small compared with U. In such cases de/ot can be omitted 
as a reasonable approximation. Even when de/dt needs to be included, 
a satisfactory approximation can often be obtained from a less accurate 
procedure than (16), namely from 


(5 
ee es 
(5 
7 : oO. 
ot 71), 
(iii) Determination of K.—'The opacity K must include energy transfer by 
both radiation and conduction. We require formulae for Krag, Keona, defined 


(17) 


by 


mm. one dT 
Ora - racp1*K rad .. : ’ 
. ' 3 — a R dM 
( 17 
Ovcona —* n8acpT*Kcona - Ri 


. 3 dM’ 


where QOraa, QOcona represent the energy transferred by radiation and by 
conduction respectively. Since O ini (10) is the total flux, equal to Oraa + Ocona 
when convection is absent (if convection occurs then (10) is not required), it 
follows that : 

K 1 a e od L (19) 


We start a discussion of K by considering the case in which there is no 
degeneracy of the electrons. ‘Then Keona can be neglected since conduction 
in a non-degenerate star is very slow indeed (Eddington, 1930). Even in this 
case the determination of a highly accurate value of A is a formidable task, 
requiring the evaluation of very complicated series (Morse, 1940). In the hand 
computation of stellar models it has been the practice to use tables for determining 
K when p, 7 and the chemical composition are specified. ‘labular information, 
such as that given by Morse or more recently by Keller and Meyerott (1955), is 
not well suited to calculation with an automatic computer however. Given a 
machine of large capacity and great speed, the ideal procedure would be to 
calculate the opacity ab initio, using a programme in which the relevant atomic 
theory was incorporated. ‘This unfortunately lies far outside our resources. 
We have therefore been obliged to use a simple algebraic representation for 
Kraa. We have followed Hoyle and Schwarzschild (1955), in using 


Kraa = 0°19p(1 + X){1 + [23°287]p/T*°}, 


for B=1, and have used 
Kraa O’1gp}! r [23°287]p/T* a5 
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for 8>1. We have found that (20) and (21) give a better approximation for Kraq 
than might perhaps be expected. The approximation is tolerable for Type II 
stars on account of the small metal content of these stars. For massive Type I 
stars, electron scattering is the main source of opacity, and this is correctly given 
by (20) and (21). ‘The main source of opacity in the deep interiors of evolving 
stars of moderate mass is also electron scattering. 

At first sight it seems as if our expressions must fail badly on account of the 
neglect of free-bound transitions in the outer regions of both dwarf and giant 
stars of ‘T'ype I, but this is not so, because such regions are convective (Osterbrock, 
1953; Hoyle and Schwarzschild, 1955). Another fortunate point is that, although 
(21) gives a poor value for 8 >2 (but not for 8 =2), we have always encountered 
electron degeneracy when 8 >2, and the important mode of energy transfer is 
then by conduction. It seems then as if (20), (21) give a just tolerable approxi- 
mation in all cases where reasonable accuracy in Krag is required. A better 
approximation would, however, be desirable if a computer of greater speed 
and storage capacity were available. 

It is important to choose an expression for Aeona that gives a suitably high 
conduction rate when there is electron degeneracy, and which represents Keona 
with reasonable precision in the region of partial degeneracy. Using conductivity 
values given by Mestel (1950) we have obtained the following algebraic 
representation 

Keona = [ — §°3710](7'/p)***7!. (22 
When there is no electron degeneracy (19) gives A=Krag, with conduction 
negligible, as it should be. When the electrons are highly degenerate (19) gives 
K=Keona, which is again correct. ‘The combined expression (19) is also 
reasonably good in the range of partial degeneracy for values of 7 between 
2x 10’ deg. K and 1o* deg. K, the errors in log,, K being never more than about 
0°05. 

(iv) The determination of \'.—With P’ defined by (1), 1’ can be expressed in 
the form (Chandrasekhar, 1938) 

1 32P*—24PP’ — 3(P’)* 
3 8P?—6PP’—-(P’)? 


Boundary conditions at the centre of a star 


I (23) 


Starting values of P, T must be specified. Since R=o, O=o at M=o, 
the integration of the above equations can then be performed outwards from 
the centre. 

Boundary conditions at the surface of a star 

(i) Surface temperatures less than 6000 deg. K.—Opacity in the photospheric 
layers arises mainly from negative hydrogen ions in this case—except possibly 
for very low surface temperatures when molecular band absorption may become 
dominant. ‘The theory of the outer layers of a star must then be very intricate 


and is not considered here. We have adapted the analysis of the negative ion 


case due to Hoyle and Schwarzschild (1955) to give a routine whereby R, O, P 
are determined at 7’'=10°deg. K. Inward integration of the equations can then 
be performed. ‘To effect such a determination two quantities must be specified 
at the photosphere. ‘These can be chosen in a number of ways; the most 
convenient were found to be the effective surface temperature of the star, 7, 


and the photospheric electron pressure, P,. The choice of T= 10° deg. Kasa 
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starting point of the inward run of the differential equations was dictated by the 
following considerations. It is desirable that the starting temperature shall be 
high enough for all the hydrogen of the stellar material to be considered as ionized. 
On the other hand it is convenient to be able to neglect the thickness of the 
sub-photospheric zone from 7, down to the starting temperature of the run 
of the differential equations. ‘This means that the starting temperature should 
not be set too high. A starting temperature of 7'= 10° deg. K was thought to 
give a suitable compromise between these two requirements. 

Writing O,, R, for the photospheric flux and radius, we make the very slight 
approximation of putting O, R at T= 10° deg. K equal to O,, R,. Our problem 
is therefore to determine O,, R,, and P at T=10°%deg. K, when 7,, P, ar 
specified. Evidently we have 

O, = racR3T,!. (24) 
lhe following equations were constructed from the analysis of Hoyle and 
Schwarzschild, the meanings of the symbols being: A,, A,, Az, A, are the 
concentrations of hydrogen, magnesium, sodium and potassium, concentrations 
being in terms of numbers of atoms, not in terms of mass. These are the elements 
1» C2, Cs, C, 


are the corresponding ionization potentials (13:54eV for hydrogen, 7-61eV for 


that make the main contributions to the density of free electrons. C 


magnesium, 5°12 for sodium, 4°32 for potassium). 


4,= C/T, + } logy T, + d;— logy, P., Lt & (25) 


defines «,, where d,; is a statistical weight factor (d, 0°48, d,= 0°04, d, 0°64, 
d, 092). R,, O, are obtained from (24) together with 


— 


Lim 1 +0) R? 
the total mass of the star, M,, being regarded as known. ‘The photospheri 


‘ «A, ) 5040 GM,)\ 
2 loge log yo » ——}> —2°8 + rss ns log io ‘) oO, (26) 


density is given by 
mn & &A 
Ro,T, > — =P,. 


1 1+4,; 
The quantities p,’ and x are determined from 
1A, 


Rp,’ T 


*r+a, 


a P 


logo (% + 071373) = — 13°2783 — § logy) (GM,/R,?) + F logy, 7, 
3°48(5040/7,,) + logio(p,'/p,). 
Value of P at T= 10° deg. K 
If x<0, logy P= 16°7561 + } logiop, — } logyy 7’, — logo py. 


If x >o, logy P= 11°855 + } logigp, + 4 logy T, — 4°169x—log i 1, 
where p, = 4/(5X, +3), X, being the value of X at the photosphere. 


Values of R, Q at T=10° deg. K 
R=R, Q=Q,. (31) 
‘This treatment assumes that f=1 in the surface layers and that 
(a) the ionization of hydrogen is effectively complete at 7’ = 10° deg. K, 
(6) the ionization energy of helium may be neglected. ‘The procedure should 
be regarded as very approximate for 7',<3500 deg. K on account of the neglect 
of molecular absorption, 
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(ii) Surface temperatures greater than 6000deg.K.—In all cases we have 
encountered p/T* is a slowly varying function of M in the outer parts of the star. 
An excellent approximation is given by neglecting the derivative of p/T* with 
respect to M. It can then be shown, after some calculation, that the value of P 
at T= 105 deg. K is given by 

login P= 5 + logy, {[0°4008] + [o-g161]x/u,}, (32) 
x being the positive root of the quadratic 
{losts3} 
Ks 


‘121 M, 
Do: oO. (33) 


wk 
[-0°6977]= + 
KF, 


12430] +1 


As before we put 

R=R,, O=Q, at T=10% deg. K. 
In this case &,, QO, are the most convenient quantities to specify at the beginning 
of an inward run of the differential equations. 

It is implicitly assumed in these equations that 8 =1 in the surface layers of 
the star. 

Condition of fit 

‘T'wo parameters (P, 7) have to be specified at the centre of a star in order 
that the differential equations can be integrated outwards, and two parameters, 
(R, QO) or (P,, T,) have to be specified at the photosphere for integration inwards 
from the surface. Thus a pair of integrations, one from the centre outwards and 
one from the surface inwards require the specification of four parameters. ‘The 
condition for the existence of a possible stellar structure is that the inward and 
outward integrations can be made to give the same values of R, QO, P, T at the 
value of M at which the two integrations meet. 

It is emphasized that with X(M), 6(M), de/et arbitrarily specified there is 
no existence theorem requiring a solution. Indeed we have found cases where 
no solution exists. Nor is there any guarantee that there shall be only one 
solution—we have found cases where there are two. 

7. Equations for an evolutionary step when X is a slowly varying function of M. 
In non-convective material we have 

OX/et = ql (34) 
which can be expressed in finite differences as 
X(M, t,, .,)=X(M, t,,)-—q' U(M, t, tao — bn) + Athair tn)’, (34') 


provided A(M,t,,,,) 0, where t, --t,,,, represents the evolutionary step, and q 
is the energy released by the transformation of one mass unit of nuclear fuel 
For B=1,q=6'2 x 10"* ergs, while for B= 2, q=8-6 x 10" ergs. Ifthe right-hand 
side of (34') turns out to be less than zero, when some particular time step is 
specified, then the value of 8 must be changed and X must be reset in accordance 
with the discussions of Section 3. 

The determination of X(M, t,,,,) for convective material 1s a more complicated 
problem. Let the boundaries of the convective zone be denoted by M’, M”. 
Provided M’, M” vary slowly with time the following approximate equation may 
be used 

M(t 


M"(t,, n 
X(M, ta4,) =(M" — My ¢ X(M,t,,) dM — q~' U(M,t,)dM > 
LJ Mt) J M\ty) J 
‘This equation neglects terms involving the time derivatives of M’, M”. It has 
been found adequate in the cases we have eneountered. 


(35) 


7 
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8. Evolution when X or B is a rapidly varying function of M.—Cases in which 
X and £ vary extremely rapidly with M occur when nuclear energy is produced 
in one or more thin shells within a star. For example a star may possess a core 
with B= 2, surrounded by an envelope with B=1. Energy generation may then 
be confined to an extremely thin shell immediately surrounding the core. 
Examples have been obtained in which the shell contains only a very tiny 
proportion of the mass of the star: in one instance we have found a shell 
containing less than 10°° M,. It is then impracticable to follow the evolution 
by the method discussed in Section 5, since this would require extremely small 
time steps to be taken. 

We proceed in such a situation by an alternative method. In the case of a 
B=2 core surrounded by a B=1 envelope, in which evolution is proceeding by 
successive layers of helium being added to the core, we simply extend the core 
to larger mass, the amount of the extension determining the time interval of the 
evolutionary step. ‘That is to say, the mass of the core is used as the evolutionary 
parameter instead of the time. If (,),, denotes the mass of the core at step n 
and (Q,),, is the corresponding luminosity then the time associated with the step 
(M.),, >(M.),, 4, 18 given by 

q(Mo)n+1— (Men }/(Qs)n (36) 
with satisfactory accuracy, provided the evolutionary steps are not too widely 
spaced, 

It must be emphasized that (36) can only be used in the case where just one 
nuclear fuel is being burned in the star. If in the case just discussed, appreciable 
energy were also being generated by “‘ helium burning”’ at the centre of the star, 
as well as in a hydrogen shell, then more careful considerations are needed. 
We can still use the mass of the core as the evolutionary parameter, but the time 
must be calculated from 


qi(M naa — (Men }/On’ (37) 


where q corresponds to the burning of hydrogen, and Q,,’ is the energy generated 
in the shell alone (equal to O, minus the energy generated by the helium). The 
quantity X(M,t), 8 =2 in this case may or may not be a slowly varying function 
of M. If it is slowly varying, then the equations of Section 7 can be used to 
determine the change of X(M,t), 8 =2 in an evolutionary step, the time interval 
of the step being determined by (37). ‘This keeps the burning of the helium 
commensurate with the burning of the hydrogen. 

A still more awkward case arises if the helium core contains an inner core of 
carbon, oxygen, and neon, for in that case helium burning also occurs in a shell, 
and X(M,t), B=2 is also a rapidly varying function of M. We can then make 
an evolutionary step by extending both the helium zone and the inner core, but 
the two extensions must be commensurate with each other. ‘This can be ensured 
by using (37) for the extension of both shells, and by arranging that the time 
given by (37) is the same in both cases. When (37) is used for the extension of 
the helium burning shell the value of g for B=2 must be used, M, must be 


interpreted as the mass of the inner core and QO,’ must be replaced by O,—(Q,’. 


xn 


Still more complicated structures involving three or more shells can be 
treated in a like fashion. 

9. The conversion of the equations to a form suitable for the EDSAC I.—In the 
EDSAC I any number is treated as if it is in the range —1 to 1, and is stored 
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so that it may have a rounding error of the order of 2. Thus since the variables 
in the equations above may change during the course of the calculation by a 
factor much greater than 2®°, the numbers cannot be stored with accuracy even 
if they are scaled by an appropriate fixed factor. 
The equations themselves suggest a method of overcoming this difficulty. 
It will be observed that most of the operations involved in the calculation of the 
derivatives are those of multiplication and division and of obtaining powers 
ot quantities. The operations of addition and subtraction are comparatively 
rare. ‘This suggests that we should represent the numbers by their logarithms 
(scaled by a certain constant factor) in the machine. We therefore write 
M=log M, R=logR and similarly for all other variables and constants which 
may occur. 
[f A+ B=Candif A and Bare both positive we define addition and subtraction 
logarithms by 
C=al (A, B), 
and B=sl (C,A). 
We give the following specimen equations derived from those above 
dR 
dM 
dP 
dM 


In this representation of the problem to the machine the operations of division 


exp {M~— 3R-— e— log, 47}, 


exp |2M—4R- e+ log, G/4r} 


and of forming powers of numbers are replaced by the much faster operations 
of subtraction and multiplication respectively. It is necessary however to provid 
the machine with subroutines for calculating logarithms and exponentials, and 
the quantitics defined by al(A, B) and sl(A, B). 

For the integration we have used a programme which 1s available in the 
EDSAC library of subroutines. This programme (Gill, 1951) uses a Runge~ Kutta 
process to integrate a set of differential equations of the form 

dy,/dt =fi( V4) Vor «++» Vn) 

It was necessary to exercise great care concerning the size of the step in the 
independent variable M. It seemed reasonable to restrict the step 5M to be 
such that the corresponding changes in the dependent logarithmic variable 
were less than a certain fixed amount. ‘The machine was therefore programmed 
to halve or double the step of integration according to the size of the derivatives. 
In order to be able to print the variables at specified tabular intervals the doubling 
of the step was only permitted when the equivalent of an even number of step: 
of the current size had been performed. ‘This ensures that every tabular point 
is reached. As it was not possible to give the machine the values M=o, R=o, 
O=o (i.e. M ©, R ©, Q ©) as the starting conditions for the 
integration at the centre of the star the integrations were started at a small valuc 
of M=M, say. The corresponding R, and QO, were then calculated by assuming 
that the density and the rate of energy generation were both constant in the 
region R<R,. M, was taken to be e*M,,. 

Solutions of the differential equations which satisfy the boundary conditions 
both at the surface and at the centre are obtained by the following iterative process 

A set I,, Ig, I, of three inward integrations is performed. J, is regarded as 
basic and J, and J, have small increments of the starting parameters, say R, 


37" 
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and Q, respectively. The integrations are stopped at an assigned value of the 
mass M=M, say. A similar set of outward integrations O,, Og, O, is also 
performed from M=M, to M=M,, Og and O, having small increments in P, 
and 7, respectively. ‘The two sets of values of R, QO, P and T at M=M, are 
then compared and it is assumed that the differences between the values obtained 
by inward and outward integration are approximately linear functions of variations 
in the starting parameters. Under this assumption it is possible to find a better 
set of initial values of the variables for the inward and outward integrations, by 
solving, using the machine, a set of four simultaneous equations. 

We then integrate the equations with the new set of initial conditions. If 
there is substantial improvement in the “fitting” of the inward and outward 
solutions we may obtain a better set of initial conditions by assuming that the 
same linear dependence on the initial conditions will still hold. This gives us 
a method of first order numerical iteration for the solution we require. 

If, however, the increments are large it may be necessary to carry out a further 
set of six integrations. ‘There is no guarantee that either of these processes will 
converge but if. we start with initial conditions which are reasonably good we 
find that three or four iterations give a fit in the solutions to an accuracy of 
one part in the fourth decimal of the logarithm. 

The EDSACI performs an integration, inward or outward, in about 
15 minutes. ‘This means that the machine is theoretically capable of obtaining 
a sufficiently accurate solution to the equations in a few hours, but in practice 
the actual time depends very much on the original estimates of the initial 
conditions. 

The time taken to calculate the change in the constitution of a star and to 
estimate de/0t is negligible in comparison with the time taken for the solution of 
the differential equations. 


University Mathematical Laboratory, St. John’s College, 
Cambridge Cambridge. 
1956 May 10 
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TYPE II STARS 
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Summar\ 


Fifteen evolutionary steps have been computed for a star of mass 
2°512» 10%* grams. Comparison of the results with the colour-magnitude 
diagram of M3 suggests that this globular cluster has an approximate age 
of 6°5 x 10° years 


Fifteen evolutionary steps have been calculated by the method described in 
the previous paper for a star of mass 2512 x 10*° grams. ‘The following initial 
composition, taken to be uniform throughout the star was used 


Hydrogen concentration 09309 

Helium concentration 0:0666 

Carbon—nitrogen concentration 0'0025 
Small concentrations of other elements could be introduced, but would have no 
appreciable etfect on the results. Carbon and nitrogen are included in the above 
list because of their importance in the thermonuclear transformation of hydrogen 


into helium. ‘The present initial composition is regarded as a reasonable estimate 
of the composition of ‘Type II stars. Metals have an important influence on 
conditions at the photosphere in the later evolutionary stages of the star, but have 


TABLE | 


Solution log of log iol log iol» logjo7's 


93090 
91925 
90642 
$9241 
87256 
55356 
83126 


10°825 33°337 
10°829 33°870 
10°840 33°904 
856 33°94! 
881 33°995 
gI§ 34045 
972 34114 
‘O17 155 
‘975 "215 
‘O70 34°300 
201 34°334 
212 356 
277 497 
"335 437 78401 
‘417 34617 é “77620 
‘576 34°593 ‘02! ‘76388 


SCnNNN NNN NI 


gi817 
‘S1255 
‘$0063 
‘79809 


NNN 


79523 
78917 


Nos) sv 


R in cm 

L in ergs per sec. 

P,, the central pressure in dynes per cm‘ 
T,, the central temperature in deg. K 


X the hydrogen concentration averaged over the whole star 
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little effect in the early stages considered in the present paper (Hoyle and 
Schwarzschild 1955). 

The calculations were performed with the aid of the computer EDSAC | at 
the Cambridge University Mathematical Laboratory. Results for the radius R, 
luminosity L, central pressure P, and central temperature 7, are givenin Table | 
The mean hydrogen concentrations for the whole star at the various evolutionary 
stages are also given.* 


‘r—- 











Fic. 1 


I'he solutions 1 to 8 were found to possess a central convective region containing 
approximately 8 per cent of the mass of the star. ‘This convective core does not 
retreat substantially until essentially all its hydrogen content becomes exhausted. 
This occurs between solutions 8 andg. ‘lhus although there is little change in Y 
(the mean proportion of hydrogen in the star) between these solutions, there is a 
highly important change in the structure of the star. At stage 8 the central regions 
still contain some hydrogen and the thermonuclear production of energy still 
occurs mainly in the central regions. But at stage g the production of energy takes 
place ina shell zone surrounding a nearly isothermal inert core. ‘The mathematical 
structure of the problem is thereby altered, almost discontinuously. It is this 
change of mathematical structure that explains the marked change in the numbers 
givenin Table l. The present results may possibly depend in this respect on the 
concentrations of carbon and nitrogen used in our calculations. With lowe: 
concentrations the importance of the carbon—nitrogen cycle would be lessened, and 
this would decrease the extent of the inner convective region. It may also be 
noted that the formula we have used for energy production by the carbon—nitrogen 
cycle is based on the work of Bosman-Crespin, Fowler, and Humblet (1954), and 


* The full tables of solutions have been deposited with the Koyal Astronomical Society. 
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that we have chosen the more rapid of the two energy production formulae 
discussed by these authors, depending on there being a favourably placed resonance 
level in the nucleus of O. If there is no such level, the slower formula for 


wesw SEQuEeNct OF wd 








Values adopted in obtaimng Fig. 2 
Observational points (mean values for M3 sequence) . 
log L/L, 0°108 0°048 0'252 0°448 
log Terr 3°713 3°730 3°757 3°782 
log L/L, 0656 0864 I°1i2 1°320 
log Terr 3°777 3°745 3°720 3°715 


Theoretical points : 

log L/L, 0°260 0°292 0°326 0°363 0418 0°470 0°536 0'577 

log Tett 3°834 3°840 3°843 8 803°844 0 3°B45 3°41 | 3°B30 0 3°B17 

log L/L_  0°637 0°722 0°756 °0°779 0°830 o0’910 1'040 1°316 

log Tere 3°853 3°827,  3°770 -3°770 3°75 3°74 3°733- 3°72 
carbon-nitrogen energy production would have to be used, and this again could 
modify the details of our solutions, especially solutions 8 and 9. 

The surface boundary condition of the star was treated in the manner discussed 
previously (1956), the switch to the surface condition of Hoyle and Schwarzschild 
being made at solution 12. 
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To facilitate comparison with our results, the stars of the globular cluster M3 
are represented in a colour-magnitude diagram in Fig. 1.* From this diagram a 
mean curve can be drawn. This is given in Fig. 2, together with our results. It 
will be noticed that the scales used in Figs 1 and 2 are different. The conversion 
was made in the following way: First, apparent magnitudes mp, were converted to 
absolute photovisual magnitudes Mp, by taking the horizontal branch stars of 
Fig. 1 ashaving Mpy=o0-o. ‘This gives o-o for the absolute photovisual magnitude 
of the RR Lyrae stars, a value that is generally accepted by astronomers 
and which is not likely to be in error by more than about two tenths of a 
magnitude. ‘lhe magnitude-luminosity conversion between Figs 1 and 2 is then 
made by putting the photovisual magnitude of the Sun equal to + 4-62 (Kuiper 
1938). We then have 


L 
logio F- 0°4 {(4°62 — Mpy) + bolometric correction}. 
“0 
Bolometric corrections are given as a function of colour index in Table Il. ‘The 
conversion from colour index to effective surface temperature is also given in 


‘Table II. 
Tasce Il 


Colour index Logarithm of effective surface Bolometric correction 
‘Temperature (deg. K) 
9550 
9075 
8680 
, 8285 
797° 
7680 
7440 
720 
}°7000 31¢ 
2 OR10 3 Re 
6550 46¢ 


Discussion of Results.—-The first point requiring emphasis is that our calculated 
points ought not to fall entirely on the observed sequence, since the observed 
sequence concerns stars of variable mass presumably of essentially the same age, 
whereas our calculated points refer to a star of fixed mass at various evolutionary 
phases. 

We next note that evolution is an accelerating process. Consider a series of 
time /, ,» ++. in the evolution, ¢, representing the initial state. ‘Then we have 


¥, dX 


t,=qM | 
l q z, L 


where X,, A,, are the mean hydrogen concentrations for the initial state and for the 
stage ¢,, respectively, M is the mass of the star, and q is the energy obtained from the 
transformation of one mass unit of hydrogen into helium (q= 6-2 x 10! ergs pe: 
gram). Computations of this integral (which can be performed by simple 


* We are much in debt to Dr A. R. Sandage tor providing the observational material of this 
paper, and also fer allowing us te reproduce it 
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quadratures with the aid of the information given in ‘lable 1) show that ¢,, increases 
only slightly beyond stage 13. The steps 1-2, 2-3, 3-4, 4-5, 5-6, 6-7, require 
nearly equal time, each amounting to approximately 2:5 x 10!’ seconds. The step 
7-8 requires about 1-5 x 10!® seconds, while thereafter the steps up to stage 13 
average somewhat less than 10'* seconds. After stage 13 the evolution becomes 
progressively even more rapid. ‘Thus by the time the giant branch is reached 
significant evolution takes place in as little as 10'° seconds. ‘The acceleration is of 
course promoted by the increasing luminosity. 

This phase of the evolution therefore occupies a very short time in the life of a 
star. As the rate of evolution of a star depends strongly on its initial mass, stars 
of different masses reach this phase at fairly widely differing ages. ‘The cluster 
stars observed in this part of the sequence, supposed of equal age, must therefore 
have closely similar masses, although there will be a wider variation in the masses 
of stars observed in lower parts of the sequence.* Hence the higher parts of the 
sequence should differ little from the evolutionary track of a star of fixed mass— 
the mass being taken equal to the average mass of the stars that populate the upper 
part of the sequence. Our results, shown in Fig. 2, confirm this expectation. 

The present interpretation (which it will eventually be possible to test when 
more evolutionary tracks have been computed—particularly for stars of smaller 
masses) requires the stars of the observed sequence at and beyond point 13 to have 
masses close to the mass for which our calculations were performed 
2°512x 10% grams. Evolution beyond point 13 is so rapid that the masses of the 
highly luminous giants of M3 and of the horizontal branch stars, including the 
RR Lyrae stars, should not be much different from this. 

Our interpretation also requires the age of M3 to be close to the time of evolution 
to point 13, which when computed as indicated above is found to be 2-049 x 10" 
seconds, or 6°5 x 10° years, in good agreement with the estimate of 6:2 « 10° years 
by Hoyle and Schwarzschild (1955). 

In conclusion it is worth remarking once again that the swing left between 
solution 8 and g is probably dependent on the existence of a resonance level in the 
© nucleus. If such a level does not exist, energy production by the carbon- 
nitrogen cycle would be much weaker in these stars (approximately half the energy 
is produced by the carbon-nitrogen cycle in our computations) and this might 
well affect the form of the solutions between points 8 and 11, although probably 
not much otherwise. ‘This again is a point that can be checked by further 
computations. 

It is also worth noting that evolution from 8 to g is very rapid indeed, as also 1s 
the evolution across from g to 10. Consequently not many stars should be found 
at these particular stages, especially the step 8-9. 

While there is some theoretical doubt whether the swing left at stage 8 ought to 
occur in the globular cluster stars, there should be no such doubt for appreciably 
more massive stars in which the carbon-nitrogen cycle is certainly the main source 
of thermonuclear energy. Hence some similar phenomenon corresponding to the 
exhaustion of a convective core ought to occur in more massive stars. One might 
hazard the guess that rapid evolution taking place at this stage will move the stars 
very markedly towards the right in a colour-magnitude diagram. Such a rapid 
sequence of changes is exactly what is required to explain the existence of the 


* This point is discussed in detail by Sandage (1953) 
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well-known Hertzsprung gap. ‘This again is a possibility that can be decided 
only by further computations. 


It is a pleasure to acknowledge our gratitude to the Director of the Uni- 
versity Mathematical Laboratory, Cambridge, for the facilities that have been 
given to us, and particularly for the use of EDSAC. 


University Mathematical Laboratory, St. Fohn's College, 
Cambridge : Cambridge. 
1956 May to. 
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THE PERIOD-LUMINOSITY RELATION AND STELLAR 
COMPOSITION 


Reddish 
(Received 1956 March 12) 


Summary 


The structure of cepheids and the effect on the period-luminosity relation 
of a change in the composition of the material out of which variables are 
formed, is considered. It is shown that if a difference AMpyg is observed 
between the period-luminosity relations of two groups of cepheids, the 
maximum difference dX of the hydrogen abundances in the initial composition 
is given by the inequality 


4 I 
a. — |dX < 0:94AM 
EE Qe |: Dé, 


where QO 1 —q™' is the fraction of the mass of a cepheid contained in its 
core exhausted of hydrogen. 

It is probable that the hydrogen abundance in the interstellar gas in M31 
does not vary by more than 15 per cent across the galaxy, and the variation 
may actually be very much less than this. 


Investigations by Baade and Swope (1) into the period-luminosity relations 
for cepheid variables in various regions of the Andromeda Nebula M31 suggest 
that any difference in the zero points of the relations are small, say less than 0-25 
magnitude (2). 

These observations lead to the conclusion either that there are only small 
differences in composition between the cepheids of a given period or that any 


differences which do exist displace the period-luminosity relation parallel to itself. 


The considerations which follow show that a change of only a few percent in the 
hydrogen abundance in the material out of which the cepheids formed would 
cause a measurable displacement in magnitude of the period-luminosity relation. 
It is concluded that any differences between the hydrogen abundances in the 
interstellar gas near the centre and towards the edge of M31 is less, probably much 
less, than 15 per cent. 

I heories of the structure of stars lead generally to an equation for the luminosity 
in the form L = Ap* MR” 2" (1 +Xy", (1) 
where yp denotes molecular weight, M and R mass and radius, Z and X abundances 
of heavy elements and hydrogen respectively. With opacity due to photoelectric 
absorption, usually k= 7-5, /= 5:5, m= —ors5 and n 1, while if opacity is due 
solely to electron scattering k= 4,/=3,m=n=o0. A varies with the structure but 
is a constant for all stars of similar structure. 

‘The theory of pulsation gives the equation for the period 


P= BR** M~"*, 


where B also is a constant for all stars of similar structure 
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Observations of cepheid variables give a relation between P and R (3) which 
may be written 
a C is a constant. Poth, (3) 
Eliminating M and R between these three equations gives 
log L=log A + 2llog B—(3l/q + m/q) log C + {l(3/q— 2) + m/q} log P 
+ klogy+niog Z —log(1 + X) (4) 
Neglecting the effect of variations in Z and X we have 
2/Alog B+ Alog A + kA logy = Alog L — {1(3/q — 2) + m/q} A log P. (5) 
lhe empirical mass-luminosity relation (§) gives /~ 3-5, close to the theoretical 
value for electron-scattering opacity. If the luminosity is almost independent of 
radius (as in the Hoyle-Lyttleton type of giant models (4)) then m Zo and 
7A log B+ Alog A +kA log wp = Alog L — 3-5 (3/g—2) Alog P (6) 

With the exception of the two points of largest P, the values of log P and log R 
given by Allen (3) are fitted well with a line of slope g=1-0. If the two highest 
points are included and if all the points had the same errors, a line of slope as high 
as 1°18 could be fitted. However, it is well known that the error of determining 
effective temperatures and thus stellar radii increases with higher luminosity class 
and thus in this case with increasing R. Consequently the slope will not be fat 
from 1-0 (providing there are not unknown systematic errors in the data which 
\llen used). 

\s the effective polytropic index mer: increases, B decreases. An increase in 
Nett, 4N increase in the central concentration of mass, usually (but not necessarily) 
implies higher luminosity, a more advanced stage of evolution and greater average 
molecular weight. Consequently A log A and A log » usually have opposite signs 
to Alog B 

With Alog P=1°6 (P increasing from 24 to 100 days), \ log L 21-6, and 
q= 1°0, if Alog B is negative 

7A log B 40, Alog B 0°57, 


or B decreases by a factor of more than 3-7 as P increases from 2} to 100 days. 
With g=1°18 7A log B. 14, Alog B 0°2, 


or B decreases by a factor of more than 1-6. ‘Thus as P increases trom 24 to100 


days B decreases by a factor between 1-6 and 3:7. ‘This suggests that convective 
transport of energy becomes more important as the period shortens. Epstein has 
shown (6) that B decreases by a factor 2°5 as the effective polytropic index negr 
increases from 1-5 to 4; this seems to be about the largest probable increase, and 
if that is the case then with g <1°18 it follows that 
AlogA+kAlogyu<t-4, 

or the increase in luminosity due to increasing central concentration (A) and 
increasing molecular weight is by less than a factor of 25. If the luminosity 
increases with radius so that m is positive, Alog A+kAlogy becomes smaller 
still. An increase in » implies a decrease in X and sometimes an increase in Z. 
With nin the range oto — 1, Z"(1 + X)~! may increase or decrease ; but the change 
is almost certainly fairly small as there is no evidence of large differences in Z 
along the cepheid sequence and nis probably close to zero. 

The alternative of Alog B positive would imply (with k=4 for electron- 
scattering models) Alog A + 4A logy < — 4:0. 
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is unlikely to decrease more than by a factor 3/8 (all He to all H) so that 
Alog A < —2°3 


An increase in B would imply an even larger decrease in A than this factor of more 
than 200. Such large changes in A and yw suggest enormous differences in the 
structure, composition and masses of the cepheids along the sequence, and do not 
appear to be as probable as the comparatively small changes which result from 
assuming that B decreases as P increases. 

A decrease in B by a factor between 1-6 and 2:5 implies an increase in effective 
polytropic index (6) of between 1-3 and 2°5. ‘The larger increase would imply an 
increase in central concentration which would largely offset the effect of decreasing 
mean density to give similar central densities along the sequence of cepheids. If 
that is the case, the connection between variability and particular central conditions 
would suggest that the onset of degeneracy may be a cause of pulsation. But 
changes in central concentration determined from the above equations depend 
critically on q through mers and B, and the figures merely suggest a possibility. 

The problem now to be considered is the effect on the period-luminosity relation 
of a change in the composition of the matter from which the cepheids are formed 

At the low densities which exist in these extended stars, the opacity is due to 
electron scattering except in the outer layers which contain only a small proportion 
of the mass and have little effect on the luminosity. Hence in equation (1) n~o 
and k~4. 

Now p= 4/(§X + 3) for small Z, and if a fraction O of the mass of a star is 
contained in a hydrogen-exhausted core, and X was the hydrogen abundance in 
the gas from which the star formed and is equal to the value in the envelope of the 
star, then writing g =(1 — Q)~! we have 


A (klog p — log (1 + X))= 0°43 (4/(X +0°6q) — 1/(X +q)} dX, (7) 


where bars denote present mean values in the star. ‘Thus the displacement of the 
period-luminosity relation in bolometric magnitude resulting directly from a 
change in the initial composition of the stars is 


AM,,) = 1° (4/(X + 06g) — 1/(X + q)} dX (8) 


An increase in the molecular weight of the matter from which the variables 
are formed means that the ratio of the molecular weight in the envelope to that in 
the hydrogen-exhausted core will be increased and the radius will decrease (4) 
Thus the period will shorten, Alog P will be negative and the displacement of the 
period-luminosity relation will be increased.* The surface temperature will rise 


* It may be that variability occurs at different stages in the evolution of stars of different ma 
Any change in the composition of the envelope may change the stage at which pulsation occurs. If 
pulsation occurs before maximum extension of the radius, it may be that an increase in the mole 
cular weight in the envelope (which will cause it to contract) will retard the stage at which pulsation 
takes place to the point where more of the core has been burned out and the resulting further 
expansion of the envelope produces a similar distribution of gravitational potential energ 
which would have existed at the pulsation stage had the composition of the envelope 
changed. ‘The overall molecular weight, and the luminosity, will have been increased still further 
provided the envelope is not more extended at the pulsation stage as a result of the change in com 
position (this is evidently unlikely) the preceding equations will remain valid. Conversely, if 
pulsation occurs after maximum extension of the radius, the stage at which pulsation occurs may be 
advanced. Proper consideration of this matter will have to await a more complete knowledge of the 
structure of cepheids 


to that 


not been 
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and the star will become photographically brighter on that account, 1. 
A Mog - AM,,,, 


Any increase in the abundance of heavy elements besides helium accompanying 
the decrease of hydrogen abundance would affect the effective radius of the star 


and the depth of the outer convective region (and thus the effective polytropic 


index of the star). Hoyle and Schwarzschild have considered such an effect in 
Population II red giants (7). An increase in the opacity in the outer layers due to 
increased abundance of heavy clements would have little effect on the luminosity 
but would tend to raise the photosphere to a higher level. With constant lumino- 
sity L= R*? 7,4 requires that 7’ R-!? in the outer layers; but 7 falls off more 
quickly than this. Consequently the surface temperature would fall too quickly 
as the photosphere was raised and the star would be unable to radiate the energy 
coming from below. ‘The density of the outer layers must decrease to allow radia 
tion to escape from deeper in the star—that is to say the effective polytropic index 
vould increase I'he constant B in equation (2) would decrease and the period 
shorten. ‘The luminosity would decrease slightly and the radius increase 
somewhat. In the case of polytropes, both uniform and composite, the effect on 
the period of the change in Busually outweighs the opposite effect of the change in 
radius, and the difference is enhanced by the fact that the outer half of the star is 
most important in determining B (6). Although the cepheids will differ in struc 
ture from simple polytrope models it is likely that these changes will be similar in 
effect. ‘The effect of an increase in heavy elements accompanying an increase in 
helium content would therefore probably increase rather than decrease the dis 
placement given by equation (8). 

The conclusion is therefore reached that if the period-luminosity relations of 
two groups of cepheids differ in magnitude by AMpg, the maximum difference of 
the hydrogen abundances in the material out of which these stars condensed is 
given by the inequality 


4 (X +4 o'6q)—1 (X 4 q)| dX o'9g AM.,,, (9) 


where q is probably about 2, (Q~}). 
With a maximum difference of 0°25 in Mpg between the cepheids near the 
centre and those half way to the edge of M31, the corresponding maximum 


difference in hydrogen abundances is 15 per cent. 


Department of Astronomy, 
University of Edinburgh 
1956 March 9 
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Summary\ 


‘The star p Velorum is a visual binary, one component ot which is a 
spectroscopic binary. Ina period of four years, commencing at the beginning 
of 1952, a total of 56 spectra has been obtained with the Cassegrain spectro- 
graph attached to the Radcliffe reflector. ‘The measures by the author have 
been combined with the earlier observations by Sanford. Using the solution 
for the visual orbit given by van den Bos, it has proved possible to eliminate 
the effects of motion in the visual orbit from velocity observations of the 
spectroscopic binary. A new solution for the spectroscopic orbit is obtained 
Three equations are obtained giving the values for the individual masses in 
terms of the parallax. Using spectra obtained at epochs of large velocity 
difference the spectral types of the components of the spectroscopic binary 
are e.timated. Microphotometer tracings of al] suitable spectra have been 
examined in an attempt to detect lines due to the visual companion Ihe 
failure of this attempt enables the spectral type of the visual companion to 
be estimated. The spectroscopic data are compared with the photometric 
observations of colour and magnitude obtained at the Cape. A value for 
the parallax not inconsistent with the trigonometrical and spectroscopic 
results is adopted and the masses of the individual stars inferred. ‘There 
seems no escape from the conclusion that the visual companion, the most 
massive of the three components, is markedly under-luminous 


Introduction.—The star, p Velorum, R.A., 10" 33.1, Declination, — 47°42 
(1g00), is a close visual binary, for which an orbit has been given by van den Bo 


(1), (2). He quotes the magnitudes as 4-49 and 5:3, the spectral types as F2, A3, 
finds a period of 16-00 years, a semi-axis major of 0"-316 and an eccentricity of 


0-70. ‘lhe spectrum shows the characteristics of a spectroscopic binary, and, as 
long ago as 1918, Sanford (3), found a period of 10°210955 days and an eccentri 

city in the spectroscopic orbit of 0-541. ‘The system is, therefore, a triple on« 

with one component of the visual binary itself a spectroscopic binary. In his 
analysis, Sanford failed to take account of the variation in the velocity of the centre 
of gravity of the spectroscopic binary due to motion in the visual orbit, and, a 

van den Bos has pointed out, this variation is too large to be negligible, so that new 
observations and a fresh discussion became desirable. Sanford mentions the 
existence of the visual companion, but, lest it be thought that the foregoing 
remarks imply any criticism of him, it should be stated at once that he was the 
victim of a double observational coincidence. ‘The two observations of the visual 
double then available were 18 years apart, which is so nearly the visual period that 
there had, apparently, been only a small change in relative position during this 
interval. His radial velocity observations fell into two groups, one near 1905 
with a three-prism spectrograph, the other near 1914 with a two-prism instrument, 
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At both epochs the velocity corrections for motion in the visual orbit chanced to 
be rather small. Sanford did find a fairly small systematic difference between the 
two groups of velocities, but he evidently put it down to systematic differences 
between the two spectrographs, and applied a correction of —2:84km/sec to 
correct the second group to the first. In handling Sanford’s results we now 
remove this correction. 

As the result of the suggestion by Dr van den Bos, special attention ha 


been paid to p Velorum in the radial velocity programme undertaken by 


TABLE | 
Sanford’s Observations 


Julian Date : . Corrected 
Velocities fk : Phase 
2 400 000 Velocities 


P P S 

Three-prism spectrograph 
16 463°345 35°3 30 
16 517°B53 38°. 33° 
16 %46°797 >] ; 29°2 
16 371°744 25 4 
16 $72 Fol 39°2 

S61 13" 
] Sic 33 
177 613°714 o's 
17 654°716 5° 


T'wo-prism spectrograph 
20 120°3837 42°7 
853 39°3 
851 38°6 
S19 27 
$04 | hy 
252°507 40°7 
20 255609 40 
20 257°4605 26 


20 260°66% 20°2 ; . 1175 


7 
20 261 ! 13 ‘ 2157 


20 205°572 41 I 597% 
20 266 2 69960 
20 207°57 22°6 7936 

SOS 

1909 

2769 
20 204 , $243 
20 296 620 
20 308 . 5: 1947 
20 309 63 C ; $952 
20 340°47° ’ 65 5 ° 9335 
20 349 . S140 
20 350°473 64 sii 9° 62 g129 
20 351° + 83 . 39° +S1 0095 
a1 260° ’ : 4°6 0729 
at 280° -33° +87 ° -38°0 83 0246 
al 290°719 ‘4 105 ‘295 44°5 +101 0000 
a1 321 °662 2 85 297 33°4 8 5205 


Sanford regards this observation as uncertain 
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Taste II 
The Cape Observations 


Julian Date ‘ , Corrected : 
Velocities F : ; Phase 
2 400 000 Velocities 


S S 
34 069 
34 196 
34 351° 
34374 
34 427° 
34 436° 
34 445 
34457 
34 475° 
34 457° 
34 499° 
34 508° 
34 519° 
34 547 
34740 


"306 ' 5432 
145 5 9520 
‘460 35° 1735 
553 4227 
‘734 ’ ; 6014 : 
‘757 +43" 4344 > 3a 
"783 + 33° 6514 2a 
S00 t . §312 2a 
‘822 36° 2875 2b 
‘831 39°3 4072 3a 
835 + 34° 6393 4b 
"836 + 39°; 5173 gb 
‘836 , 5910 3b 
827 38° 3282 + tb 
661 34° 2656 2b 
34.755 647 ‘ 7388 } 2b 
34 770 634 7 2009 2b 
34738 618 2 +2 9935 2b 
791 614 4°" 2545 2b 
5 108°s8:z 4°77 +O°402 3127 7 ib 
5 123°5! ‘ 396 7519 , ib 
156 2° 8 380 gQs2 tb 
171 374 46017 ' tb 
1912 2 363 s° 4130 tb 
191°453 24° 363 4290 ib 
196°27 2 361 go1o & tb 


ooooooo0o00000 0000 


~ = 
oO © 


° 


3 


196° 2 58°" 361 5 gIol t tb 
196 d 71 361 *4 76°3 9205 Pas ib 
441 ¢ 69 “20% 2 ie 9277 4 tb 
433 °538 82 +O'252 26°7 87 0356 22 4b 


AuwuuUWNN aoa 


Www ww Ww 


‘This observation was obtained with the c’’ camera, {/3°7, giving a dispersion of 


A/mm It is discordant and is omitted from the discussion 


the Cape Observatory using the equipment of the Radcliffe Observatory, 
Pretoria. 

The observations.—TVables | and II above give complete schedules of the 
observed radial velocities. In the case of Sanford’s observations, results from 
the three-prism spectrograph are counted with double weight. ‘The high reso- 
lution of this instrument, which had a dispersion of 10-4 A/mm at Hy (4), enabled 
him to resolve the stellar lines near the minor maximum of the velocity difference 
between the components of the spectroscopic binary. ‘This has not been possible 
with the two-prism Pretoria instrument. ‘The Pretoria spectra have been taken 
either with the “a’’ camera, {/8, giving a dispersion of 21 A/mm at Hy, or with 
the ‘‘b”’ camera, //6, giving 29 A mm (5). Spectra on each of thes dispersions 
have been given unit weight. 

Comments on the tables.—The results used in ‘Table | are extracted from 
Sanford’s paper, and have been changed only by the removal of the systematic 
correction which he applied, ‘The velocities in this form are given in the second 
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and third columns of the table. In ‘Table I] are the Cape results, all the measures 
being by the author \ discussion of the systematic errors of the Cape radial 
velocities, to be published elsewhere, indicates that, to correct velocities measured 
by D. S. E. to the Lick three-prism system, the following corrections should be 
applied: ‘‘a’’ camera, +1°7km/sec; “b”’ camera, +0-7km/sec. These 
corrections have already been applied to the data exhibited in ‘Table II. It 
frequently happened that several spectra of p Velorum were obtained on one plate. 
In such a case a mean result for the several spectra has been used. ‘This is 
indicated in the last column of ‘Table I1, where a designation such as 3b indicates 
the mean from three “‘b’’ camera spectra. Means have not been taken when 
the velocities were changing rapidly. 
It will be noted that the modern observations do not agree as they stand with 
those of Sanford. ‘This is interpreted as being due to orbital motion in the 
visual orbit. ‘lhe two sets of observations are capable of being reconciled within 
the observational uncertainties in the following way. Van den Bos (2) has given 
a table of the function pV, p being the parallax in seconds of arc, and V the radial 
velocity in km/sec of one component relative to the other in the visual orbit. 
For the epoch of each radial velocity observation the corresponding value of pV 
can be found by interpolation from van den Bos’ table. ‘These are the values 
denoted by F which appear in the fourth columns of ‘Tables Land I]. ‘The upper 
signs in van den Bos’ table have been used throughout. By trial and error 
some multiple of these numbers, /’, must be found, so that when these multiples 
of Fare applic d to the observed velocities, the old and the new values agree. Atthe 
ame time the period must be redetermined on the basis of all the available 
observations. ‘lhe solution was reached by trial and error. ‘The period is very 
closely determined from the epochs of large velocity difference. Adopting the 
period so determined the phases can be computed and a first value of the multiplier 
of F which best reconciles the old and new observations can be found. Appli- 
cation of these corrections leads to changes in the velocities near maximum 
velocity difference, requiring a slight revision of the period, followed by a slight 
revision of the value of the multiplier. ‘These procedures lead to a value for the 
period of 10°21040 days (compared with Sanford’s value of 10°210955 days) and 
the phases given in the seventh column are computed using this period and the 
epoch J.D. 2421290°719 as corresponding to zero phase. ‘The period is not 
introduced as an unknown in the computations which follow. ‘The value of the 
factor adopted in the correction for motion in the visual orbit is +14. It is not 
very closely determined, but it 1s thought that it could not be in error by as much 
as 20 per cent either way. Certainly values of +12 and +16 when tried gave an 
accordance noticably worse than the value adopted. ‘The corrected velocities in 
the fifth and sixth columns are the values obtained after the application of this 


correction. ‘The element which has been removed is the velocity of the centre of 


gravity of the spectroscopic binary relative to the centre of gravity of the triple 


ystem. If we assume for the present (the opposite assumption will be discussed 
later) that the spectroscopic binary is the visual primary, then the fact that the 
correction factor is positive identifies 43° as the longitude of the ascending node in 
the visual orbit. Moreover, we clearly have the relation 


I mM, 


— 14 
p m, +Myg+m, 
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where m, and m, are the masses of the spectroscopic components, and m, the mass 
of the visual companion. 

The data of Tables | and II with these corrections applied are plotted in 
Fig. 1. ‘lhe probable error of a single observation is usually between 1-o and 
1's km/sec, and, since the observations are numerous, one or two deviations from 
the mean which amount to several times the standard error occur, but the con- 
sistency is quite in accord with the usual statistical behaviour of radial velocity 
measures. 

From this figure it will be seen that there are parts of the empirical velocity 
curves where the velocity difference between the components is changing at the 
rate of roughly 6km/sec per hour. With a trigonometric parallax according to 
the Yale Catalogue of 0-033 + 0”:007 and a semi-axis major in the visual orbit of 
0”:316 according to van den Bos (1) the light-time across the relative visual orbit 
approximates to 2) hours. ‘The light-time from the centre of gravity of the 
spectroscopic system to the plane though the centre of gravity of the triple system 
parallel to the plane of the sky is some fraction of this, but is not so small as to be 
negligible on the steeper parts of the velocity curves. | am indebted to the kind- 
ness of Dr van den Bos for providing a correction curve tor this effect. ‘The scale 
of this curve is determined by the empirical relation of equation (1) for the same 
combination of symbols occurs in the formulae. ‘lhe total range of light-time 
involved is 50 minutes. From Dr van den Bos’ curve the corrections of Column 
8 in the tables are derived. ‘These are corrections to the phases in units of the 


fourth place of decimals. 


The determination of the elements.—The following table gives the normal 


points used in the solution for the elements of the primary he corrections to 


velocities and phases described above are included, 
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Taser III 
Normal Points 
Observed Computed Residuals No. of 
Phase (¢) Velocity Velocity Oo-C Observations 
1005 11°54 11°2 0°29 5 
2011 32°00 31 0°47 5 
2848 30°70 37°5: 0°83 14 
4236 40°70 54 16 6 
5935 39°49 14 65 15 
6142 37°79 37°43 "36 
6953 33°33 33°35 "02 
793° 25°35 ‘06 ‘29 
8908 ° 37 1°42 
9369 13°34 18°95 o’61 
0247 33°21 32°98 0°23 


The computed velocities are derived from the final elements given below. They 
were obtained by the customary least-squares procedure for ihe derivation of 
corrections to the initial elements stated. ‘Ihe value of y found from the least 
squares solution has been increased by 0-03 km/sec to give the adopted final value. 
In carrying out the solution it was found best, after some trials, to give each of the 
normal points the same weight. ‘The strictly correct statistical procedure 


produces a solution which passes very exactly through the normal points near phase 


o’5 at the expense of rather large residuals for the points near zero phase. Since 
the former group exercises much less influence on the determination of the 
elements than the latter does, given the large eccentricity and the orientation of 
the orbit, the policy of uniform weights seemed best, even though the formal 
probable errors of the solution are probably somewhat increased in consequence. 


‘TABLE LV 
Solution : Primary 
Initial Final Sanford 
21°22 km/sec + 21°30 km/sec 19°25 km/sec 
42°70 km/sec 43°33 + 0°76 km/sec 42°34 km/sec 
0°56 0°5597 + 00055 0°54! 
186 187°°32+0 °71 184° "62 
0'0044 0'0017 + 0'0010 
equivalent to J.D 2 420 259°452 > 420 259°381 
r 10°21040 days 10°210955 
a sini 5°042 « 10° km 4°981 - 10° km 


m sin’ 2 0°304 * sun 0°251 sun 


The average probable error of a single radial velocity observation of unit weight 
is+ 2:0 km/sec. ‘lhe formal probable errors deduced from the interagreement of 
lines for the Radcliffe plates range from + 0-7 km/sec to + 2:1 km/sec. Phases in 
our solution are computed from the arbitrary epoch J.D. 2 421 2g0°719 and ¢, 
denotes the phase of periastron passage in this system. 

For the spectroscopic secondary we merely have to find that value of K’, w 
being increased by 180 , which fits the observations best. As Sanford remarked 
in his paper he found some discrepancies, and we naturally meet difficulties of the 
same type. ‘The best value of A’ is 53°6km/sec. This is the value which makes 
the sum of the squares of the residuals, ‘* Observed minus Calculated’, a minimum, 
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‘The value which makes the sum of the residuals, with their proper signs, zero, 1s 
52;2km/sec. If, from ‘lables I and II, we extract the observations of the 
secondary in order as given there, we have the following table of residuals and 
data for the secondary. 


TABLE V 
Data for spectroscopic secondary) 
Residuals; Observed minus Calculated 
Sanford’s Observations Cape Observations 
16°7 


/ 


13 


Elements 
Sanford’s values 
AK 53°6 km sec 52 km/sec 
6°237. 10° km 6°107 10° km 


0°245 sun 0°235 sun 


The foregoing tables constitute the ordinary solution of the spectroscopic 
orbit. ‘The probable errors of the elements are rather disappointingly large, 
but there seems no point in attempting a closer fit. In any case, as can be seen, 
a high formal accuracy does not necessarily imply that the solution will give 
correct predictions. So far as this solution goes we can probably assume quite 
safely that the mass ratio in the spectroscopic primary could not possibly be in 
error by as much as 10 per cent, while the figure in equation (1) can hardly be as 
much as 20 per cent in error. 

Discussion of the components of the triple system.—\t was initially hoped to 
discover from the spectrograms evidence of the spectrum of the third star, the 
visual companion, and, to that end all suitable plates were put through the 
Radcliffe microphotometer. Special attention was paid to the Balmer lines in 
the hope that a careful study at different phases of the spectroscopic binary 
might lead to a determination of the relative positions of the lines of the third star, 
and thus to a velocity measurement and a parallax determination. ‘This hope 
has not been realised. A careful study of the spectra and tracings has failed to 
show any metallic lines due to the visual companion. However, this negative 
result allows us to make a positive deduction concerning the spectrum of the 
companion. ‘lhe spectroscopic secondary has a spectrum rather similar to the 
spectroscopic primary, and, speaking at this stage, roughly, must be of the order 
of one magnitude fainter than the primary. ‘lhe visual secondary is roughly 
half a magnitude to one magnitude fainter than the visual primary. In other 


words, the visual companion must be intermediate in brightness between the 


other two stars, and, if it had a spec trum containing strong metallic lines, or if it 


were similar to the spectra of the other two stars, we hould have no difficulty in 


detecting it. lhe only possible solution ould seem to be to designate the visual 
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companion as an early A-type dwarf, the spectrum of which consists of nothing 
but broad hydrogen and calcium lines which are too broad to be distinguished in 
the joint spectrum. ‘The joint spectrum does have broad hydrogen lines—a 
fact noted so long ago as the compilation of the Henry Draper Catalogue—and 
these we ascribe to the visual secondary. ‘he hydrogen lines of the other two 
stars appear as narrow absorptions, readily separable at high velocity differences, 
on top of these broad lines. 

‘The spectra at maximum velocity difference have been classified according to 
the MK system (6), a task of some difficulty in view of their complexity. What 
emerges is that the two components of the spectroscopic binary are of nearly the 
same spectral class with the fainter, possibly, very slightly the earlier and the 
brighter, one luminosity class higher than the fainter. It does not seem possible 
for the types to be earlier than F3 or later than k'5. Provisionally they are denoted 
as 3 —51V and k3—5¥V. In terms of a calibration to be published elsewhere, we 
have, on the S system for stars on this part of the Russell-Hertzsprung diagram, 
absolute magnitudes near + 2:5 and + 3-5 and colours near +0-12 and +0°17 for 
these stars. ‘hese values of the magnitudes give for the spectroscopic pair 
together the absolute magnitude of 2°14. ‘lhe visual companion is certainly 
fainter than the primary by an amount variously quoted as 0-5 and 0-8 magnitudes. 
If we take a value of 0-65 the absolute magnitude of the companion is + 2°79. 
‘The observed joint colour of all three is + 0°02, so that the colour index of the 
companion on this set of assumptions must be — 0-21, the colour of an Ao V star. 
If the difference is taken as 0-5 magnitudes, the companion must have absolute 
magnitude 2-64 and colour o-1g, which is that of an At V star. 

The result of this discussion is more closely determined than the uncertainties 


of the data might suggest. ‘lhe components of the spectroscopic binary certainly 


have positive colour indices. ‘The visual companion is fainter than the primary 
and the joint colour index is small and positive. Hence the visual secondary 


must have a negative colour index, which is what we infer both from this argu- 
ment and from the appearance of the joint spectrum. ‘The relative luminosities 
of the components of the spectroscopic binary are also fixed within limits, at 
roughly the value adopted, by the observed line-strengths in the two spectra. An 
important Consequence of this discussion is that the visual companion 1s under- 
luminous for its type, according to the foregoing figures, by about a magnitude 
and a half. 

It might seem more plausible to make the visual secondary the spectroscopic 
binary to avoid this conclusion of under-luminosity. ‘laking rough figures we 
may argue the hypothesis through as follows :—‘The magnitude difference between 
primary and secondary (visual) is, say, half a magnitude. ‘lhe difference between 
the two components of the spectroscopic pair is one magnitude. If the visual 
primary is of magnitude + 2-0, the visual secondary is + 2°5, and this now being 
assumed to be the spectroscopic binary 1s composed of stars of magnitudes 
and 3-9. ‘Thetwobrighter stars are together equivalent to a star of about magni- 
tude +1°5. ‘The lines of the spectroscopic companion are therefore seen against 
a total light from other sources of a star of this magnitude, and the minimum 
possible residual intensity even for a perfectly black line is over go per cent. ‘This 
does not seem to accord with the fact that the lines of the secondary are, though 
fairly faint, perfectly easily visible in the joint spectrum. 
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We therefore are confirmed in our original assumption that the visual primary 


is the spectroscopic binary and we have in consequence the following system of 
data: 


Spectroscopic pals 

K3-s51V M 

F3-5V M 
Visual companion 

Ao—1V M 
‘Total joint M= + 1°65 
Observed m= + 3°89, observed C.J. = + 0°02 
Spectroscopic parallax = 0":036, from above data. 
lrigonometric parallax = 0"-033, from Yale Catalogue. 
If we adopt a value of 0"-035 for the parallax we have, from equation (1), 

O'S51 Mz = 0°49 (m, + m,) 

and trom the spectroscopic solution 


m,/Ms = 0°304/0°245 241. 
Van den Bos’ data give, 
m, + Mm, + mM, = (0°316)9/(0°035)*(16°0)* = 2°875 
I‘hese lead to the following values for the massest: 
m,=O'S1 m,=0'05 m,=1°'41, 

and thence we have easily the values of sin z in the spectroscopic orbit, and of the 
semi-axes major of the orbits of the spectroscopic components. 

l'wo points should, however, strike the reader One is that the 


deduced 


values of the masses are rather low: the other is that the whole of the foregoing 


argument depends on the adopted value for the parallax in a rather sensitive way 
Values which are quite indistinguishable in the field of trigonometric deter 


minations of parallax produce a considerable change in the mass values, Even 
the method of spectroscopic parallaxes is not quite definitive enough. If 
leave the parallax disposable and write top 


terms of X: 


we 


1/X, we can express the masses in 


m, = 0'0682X% — o-0g5gX* 
mM, = 0°0550X* — 0:0772X" 
Ms =0'1720X* 
giving the following system of values 
X 2°90 
Pp O50 
mM, O17 
O'l4 79 
oO 69 1°55 
It is regrettable that we should have to exercise a choic 
avoiding this would be to find and measure the 


and the only method of 
velocity of a narrow line or lines 
in the spectrum of the visual companion. 

In seeking grounds on which the best choice of parallax may be based, we must 
first take note of an inescapable anomaly ‘| he ordet ot de reasiny brightness 


of the stars 1s I, 2, 2, but this is only the order of the inferred masses for values of 
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X from about 4:0 to 5-0. ‘Thus the mass luminosity relation is certainly not 
satisfied by all the stars outside this range. Inside this range the assumption 
that they all satisfy the mass-luminosity relation leads to a contradiction. For 
example, take X= 4:0. ‘The bolometric magnitudes of the components inferred 
from the mass-luminosity relation are about +0°5, +1°4, +0°5, which when 
combined and compared with the apparent bolometric magnitude give a parallax 
corresponding to X = 6. 

If we decide that the one star which, according to current ideas, ought to satisfy 
the mass-luminosity relation, does in fact do so, we can obtain a very precise 
values of the parallax, ‘This is the F3—5 V star. For it we have 

M,,., = 5 — 10 log (mass) = m,,,,— 5 log X 


or, using the expression for the mass derived above 


O°5 — OvI m,,, = (3/2) log X + log (0-05 50 X — 0:0772). 


‘The value of X which satisfies this equation depends on estimating the apparent 


bolometric magnitude of star No 2, but is not critically dependent on the value 
adopted. If we take m,,,=5:°7 we have X= 3°6. Slightly varying assumptions 
as to the magnitude differences between the components permit small variations 
in this value, up to about o-1. We end by adopting X = 3°5, corresponding to a 
parallax of 0-029. If we maintain the magnitude difference between the Class 
V and Class IV stars at exactly 1-0 magnitudes this decision gives a total absolut« 
magnitude for the system of + 1°20, and absolute magnitudes (photographic) for 
the components of 2-0, 3:0, and 2:2. On the other hand, as is explained in the 
following paper, an alternative is to hold star No 2 strictly on the main sequence 
and to solve for the magnitude difference as an unknown. ‘lhe data do in fact 
hang together so well that this does not make differences which are important from 
the standpoint of the location of the stars on the Hertzsprung~ Russell diagram. 
Adoption of 0”-o29 for the parallax gives, for the spectroscopic pair, 1 = 34°00’, 
a, =g'02 x 10°km, a,= 11°15 x 10°km. 

Acknowledgments.—I\ am indebted to colleagues at the Cape, Radcliffe and 
Union Observatories for helpful discussion: most particularly to Dr van den 
Bos for drawing our attention to this system and for wise counsel in the treatment 
of the more involved difficulties, and to Dr Halton C. Arp for his help in discussing 
the photometric data. 
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p VELORUM AND STELLAR EVOLUTION 


Halton C. Arp and David S. Evans 


(Communicated by H.M. Astronomer at the Cape 
(Received 1956 Ma 16) 


Summar\ 


The data derived for p Velorum are in good agreement with the present 
theory of stellar evolution. ‘The age of the system is about 10° vears. It is 
concluded that, in that time, the faintest component has remained essentially 
unchanged, the brightest component has completed the first stage of its evolu- 
tion, and the most massive component has advanced far in its evolution 
The position of this third component on the Hertzsprung Russell diagram is 
an indication of the probable course of the late stages of evolution of star 
of this type. 


Introduction. he data for p Velorum deduced in the preceding paper are 
unusually complete, and provide an opportunity for discussion from th 
evolutionary standpoint. ‘I'wo of the component stars lie off the main sequence 
and their masses do not satisfy the mass—luminosity relation. Present evolutionary 
theory demands just such divergences, and p Velorum provides an opportunity 
for observational test of the theory by comparison with data for a group of stars 
which have a common age. 

Observed and derived parameters.—Vhe data for the three components are 
listed in ‘lable |. Photometric data have been converted from the 1953 S system 
used in the preceding paper to the B, V system. ‘lo the accuracy needed for this 


discussion this means that corrections of + 0-2 have been applied to the magnitudes, 


+024 to the F star colours and + 0°28 to the A star colour. ‘The step which 


Tasie | 


Adopted 
parallax : 0” 029 Joint My 112 Joint B— I 
Star No. 1 Star No. 2 Star No. 3 
Mass 1°74 +0°08 1*40 2*11+ 0°40 
Uy 1'goto'l 3°30+ 073 240+0°2 
B—I +0'°37+ 0°03 0°41 + 0°03 0°06 + 0°05 
Sp F3—5I\ r3—5\ (A) 


fixes the masses given in this table is the choice of 0-029 for the parallax. ‘his 
point has been discussed at length in the preceding paper. It is the choice which puts 
the F3 — 5 V star in the position where it ought to lie on the Hertzsprung- Russell 
diagram, that is on the main sequence. ‘The main sequence chosen is the main 
sequence for unevolved stars in galactic clusters determined by Johnson (1). 
The colour of star No. 1 was then chosen as that of an F3 —51V star, and the 
difference in visual magnitude between star No. 3 and Nos. 1 and 2 together 
was chosen as 0°8 magnitudes. A solution is then made for the absolute visual 
magnitudes of Nos. 1 and 3 and for the colour of No. 3. 





545 Halton C. Arp and David S. Evans Vol. 116 


In Fig. 1, the filled circles represent the present positions of the stars on the 


absolute magnitude-colour diagram. ‘The two open circles represent the 
positions on the main sequence from which stars Nos. 1 and 3 must have come 
if they have evolved without appreciable change of mass. Lines through these 
points show the range of uncertainty which represents the lack of definition of 
the empirical mass-luminosity relation over small segments. 


The component of Pp Velorum on the Hertz prung Russell diavram 


) original positions 
* present position 


Kestrmated lromits of data indicated by horizontal and vertical lines 


Accuracy of the parameters.—'Vhe uncertainties associated with the quantities 
listed in ‘Table I were estimated in the following way: the adopted parallax 
fixes the mass of star No. 2. ‘The uncertainty in the mass of No. 1 is then the 
uncertainty of the ratio, A/A’, of the velocity ranges in the previous paper, and 
this is fairly exactly known, whereas mass No. 3 comes from equation (1) of the 
previous paper, and the right-hand side might be uncertain by as much as 
20 per cent. ‘The accuracy with which star No. 2 may be placed, in absolute 
magnitude, on the main sequence, from its spectral class is estimated as 
+o-3 magnitudes. ‘lhe uncertainties in the B—V values correspond to the 
uncertainties in spectral classification. 

Since No. 2 is the faintest of the three stars, the derived magnitudes of 
Nos. 1 and 3 are relatively insensitive to any error in the adopted magnitude 
of No. 2. ‘The colour derived for No. 3 depends primarily on the adopted value 
of the magnitude difference between the visual primary and secondary. 
Van den Bos estimates that the difference is 0-8 and that it cannot lie outside 
the range from o-5 to 1-0 magnitudes. In Fig. 1, the diagonal line drawn 
through the plotted position of No. 3 represents the range of solutions corre- 
sponding to the extreme assumptions of o-6 and 1-0 magnitudes difference. 
The estimated uncertainties of placing the remaining two stars on the diagram 
of Fig. 1 are shown by the lines through them. 





No. 5, 1956 p Velorum and stellar evolution 549 


It should be emphasized that, although uncertainties exist, they are not so 
great as to cause any significant uncertainty in the place of the stars in the 
colour—absolute magnitude diagram. 

In addition there is a check afforded by the spectrographic observations. 
Spectroscopically No. 1 is a luminosity class IV star and also, from the relative 
line strengths on the microphotometer tracings, a little more than a magnitude 
brighter than No. 2. No. 3 has broad hydrogen lines and metal lines which 
are either absent or very weak. In the preceding paper a solution was made, 
not on the assumption that No. 2 was a main sequence star, but instead on the 
assumption of 1-0 magnitude difference between No. 1 and No. 2. ‘The derived 
magnitude of No. 2 agrees with the value used here within a few tenths of a 
magnitude. 

It may be remarked that p Velorum lies at galactic latitude +10 roughly 
in the direction of 7 Carinae. At a distance of only 35 parsecs, however, it ts 
presumed to be unreddened. 

Evolution.—A star, initially on the main sequence and chemically homogeneous, 
does not alter its luminosity or effective temperature rapidly until an appreciable 
proportion of the hydrogen in its interior has been transmuted into heavier 
elements. Schénberg and Chandrasekhar (2) have shown that the star then 
begins to evolve in one of two ways. 

In the case where the diminution of hydrogen in the interior halts convection, 
the hydrogen-exhausted core is isothermal, and the star brightens by about 
one magnitude remaining at about the same effective surface temperature. It 
thereafter goes on to more rapid evolution through non-equilibrium configura- 
tions. In terms of Fig. 1, this type of evolution means a movement away from 
the initial position on the main sequence in an upward direction. ‘This case has 
been applied by Sandage and Schwarzschild (3) and Sandage (4) to stars in 
globular clusters which are initially fainter than absolute magnitude + 4. It ts 
evident that this case does not apply to p Velorum, for, although star No. 1 has 
not become brighter by as much as a magnitude, it has become considerably 
redder. 

In the case in which the stellar core remains convective, the luminosity 
increases somewhat in the early stages of evolution, but the radius increases 
rapidly, so that the predominant change 1s reddening of the star. ‘his cas« 
fits the observed position of star No. 1 very well. During the time in which the 
ratio of the mean molecular weight inside the core to the mean molecular weight 
in the rest of the star increases from 1 to 2, the luminosity increases by a 
factor 1°41, and the effective temperature at the surtace decreases by a factor 1°18. 
The actual evolutionary path of star No. 1 is roughly indicated by the solid line 
in Fig. 1. ‘The dashed line immediately beneath it represents the approximate 
path computed by Schonberg and Chandrasekhar, and terminates at the 
above-mentioned limit. Schonberg and Chandrasekhar point out, however, 
that in a real star the reduction of the hydrogen content of the core will causé 
the central temperature to rise, and hence counteract the tendency of the radius 
to increase. ‘They predict that the actual path should be similar to the dashed 
line in Fig. 1, but should lie somewhat above it. Roy (§) computes a similar 
initial evolutionary path. Since this corresponds closely with the observed 
path, the observational data for p Velorum confirm rather exactly the predictions 
of the theory for stars with initially convective cores. 
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It should be noted that this does not necessarily imply that stars of fainter 
absolute magnitude do not evolve by way of the initially isothermal core case. 
‘The rapidity of mixing would be expected to increase with the initial brightness 
of the star. 

Since convection is impossible in an exhausted core, the growth of such a 
core slows and finally stops convection. ‘The star then goes over to the isothermal 
core evolution and thereafter the principal change is an increase in luminosity. 
Although star No. 1 has, apparently, started its evolution with a convective core, 
it is not possible to decide whether it has yet gone over to the isothermal core 
case. For this star, an increase in luminosity by one magnitude from the initial 
conditions would bring it to My=1°5, and for the convective core case to 
something brighter than My=2:1. From the observed value, My= 1-9, it 1s 
concluded that star No. 1 is approximately at the end of its initial evolutionary 
track in either case. 

At the end of the initial evolutionary track, 12 per cent of the mass of the star 
is contained in the core (3). Placing star No. 1 at that limit enables us to compute 
the age of the system from the equation given by Sandage and Schwarzschild (3): 


Lt = 0-007¢?X ,q,M. (1) 


‘Laking the fraction of the mass in the core, g, =0°12, the abundance by weight 
of hydrogen outside the core, X, = 0°596, the mass of star No. 1, M= 1°74 x Sun, 
and its initial luminosity, L, to correspond to an absolute bolometric magnitude 
of 2-4, we have 


t = 1°4 x 10" years = age of p Velorum. 


‘This approximate age for the system enables us to compute an upper limit for 
the mass in the exhausted core of the main sequence star, No. 2. Since the 
evolution process proceeds faster as the star becomes brighter, the mass and 
luminosity of star No, 2 indicate that it must have less than 6 per cent of its mass 
in this core, and, possibly, considerably less than this. ‘lhe star has thus not 
had enough time to evolve appreciably away from the main sequence because 
of its low luminosity. ‘This is consistent with the original assumption that it 
still, in fact, lies on the main sequence. 

Finally, star No. 3 may be discussed from the standpoint of its evolution. 
An observational extension of the theory now becomes possible from a con- 
sideration of this sub-luminous A star. In Fig. 1 the approximate initial path 
of its evolution is shown by the line of crosses, which is drawn approximately 
parallel to the observed path for star No. 1, taking into account the bolometric 
corrections necessary at these colour indices. ‘The star rises from its initial 
position on the main sequence to about My = + 1 at the Schonberg—Chandrasekhat 
limit. From this point onwards it must evolve rapidly. Its luminosity increases 


by a large factor while its effective temperature decreases somewhat. ‘The 


beginning of this path, above My = +1 is roughly indicated by the extension, 
upwards and to the right, of the line of crosses in Fig. 1. 

It is computed from equation (1) that this star requires 8 x 10" years to reach 
the Schonberg—Chandrasekhar limit. Hence only 6 x 10° years have been spent 
in the evolution through the upper parts of the colour-magnitude diagram down 
to its present position. ‘This is in good agreement with the relative track times 
computed by Sandage (3) in the globular cluster M3. 
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Conclusions drawn from the observed colour-magnitude diagrams of globular 
clusters indicate that the stars in them reach M, 3 and B—-V= +1-6in then 
evolutionary development. ‘Thereafter they become less luminous and move 
toward bluer colour indices in the vicinity of the horizontal branch (region of 
M,=o and B—V=o-'0 to 0-2). Observations of globular cluster stars also 
indicate that, so far as those stars are concerned, the limit of evolution so far 
observed brings the stars to a region of the diagram lying to the left and below 
the region defined above. In M13 (6) the limit is in the region My = + 2°5 to 
+370 and B—V O°4. 

Star No. 3 has, evidently, come down from the more luminous regions of 
the diagram in a way analagous to the behaviour of the globular cluster stars. 
The track is suggested by the downward line of crosses. ‘The star, however, 
appears to be somewhat redder than corresponding stars in globular clusters 
at the same absolute visual magnitude. 

Hoyle and Schwarzschild (7), however, have recently developed in detail 
the theoretical evolutionary behaviour of stars to the limit of the giant branch. 
On the basis of the higher metal to hydrogen ratio, believed to be characteristic 


of type I stars relative to type II, they also obtain the observed galactic cluster 


giant branches. ‘The latter branches are a factor of about 10 less luminous. 
From consideration of the age and spectral characteristics of p Velorum it is 
probable that the component stars are of type I. ‘Therefore the subsequent 
evolution from the giant branch is not necessarily exactly the same as that 
observed for type II, and in any case they start this phase from fainter luminosities. 

The observational data thus suggest that, for a type | system of age around 
10° years, the far-evolved stars pass through a stage analagous to that observed 
for globular cluster stars. ‘They are, however, redder at the same absolute 
magnitude, and are best characterized at this stage of their evolution as A type 
sub-dwarfs. 


Indiana University Observatory, Royal Observatory, 
Indiana, U.S.A.: Cape of Good Hope. 
1956 April 
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A STUDY OF ‘THE FOUR ECLIPSING BINARY SYSTEMS: 
RW MONOCEROTIS, RW GEMINORUM, U CORONAE BOREALIS, 
AND ‘TY PEGASI 


Alan H. Batten 


(Received 1956 June 22) 


Summary 


\ new attempt at the determination of the photometric elements of the 
eclipsing binaries RW Monocerotis, RW Geminorum, U Coronae Borealis 
and ‘TY Pegasi has been made from visual photometric light curves obtained 
by the late Drs Dugan and Pierce and published posthumously by Wood. 
The first three systems appear to exhibit total exclipses, and ‘TTY Pegasi is 
a partially eclipsing system. A notable difficulty in the solutions for RW Mon 
and RW Gem has been to reconcile the “ shape ’’ and “ depth’”’ ratios of 
the radii k of the two components ; and, in the latter part of the paper, an 
attempt has been made to demonstrate the possibility that this isa consequence 
of the considerable distortion of the secondary componet.ts which, in each 
case, probably nearly fill the largest closed equipotential capable of containing 
their mass 

The photometric observations are combined with spectroscopic evidence 
available for three of these systems, and the contact hypothesis is used to 
estimate the mass-ratio of the systems which leads to the determination of 
their absolute dimensions 


1. Introduction._Vhe immediate object of this investigation has been the 
analysis of light curves of the four eclipsing systems RW Monocerotis, 
RW Geminorum, U Coronae Borealis, and ‘TY Pegasi. It has also been 
concerned, however, with a more general consideration of the problem presented 
by the light curves of highly distorted binary stars. ‘These four systems were 


all chosen because it was suspected that their secondary components occupy 


completely the space enclosed by the largest closed equipotential surface capable 
of containing their mass. ‘Throughout this paper such systems are termed 
‘semi-detached systems”’ (1), and the secondary components are described 
as being at their respective Roche limits (2). Stars at their Roche limits are 
necessarily highly distorted, and do not possess symmetry about their axes of 
rotation. ‘The question of the stability and evolution of these systems is at 
present of considerable interest, and new evidence of the dimensions and physical 
characteristics of the secondary components of additional systems of this type 
is very desirable. 

2. Previous investigations of the systems: source and description of the present 
investigation. All of these systems have been the subjects of previous analyses 
for their geometrical elements on the basis of less substantial photometric 
evidence. ‘The best known of the four is U Coronae Borealis which has been 
investigated by, amongst others, Fetlaar (3), Baker (4), Shapley (5), and Dugan 
and Wright (6); this last investigation is also concerned with the evidence of 
variability of the period. All these analyses resulted in solutions corresponding 
to partial eclipses. On the other hand, the catalogues both of Schneller (7) 
and of Kukarkin and Parenago (8) list the system as exhibiting a total eclipse 
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at primary minimum. As will be seen, this latter alternative has been confirmed 
by the present investigation. A second notable feature of the system is the 
discrepancy between the spectral type of the secondary component, as computed 
from the depths of the minima of the light curve and as deduced from the 
spectroscopic observations. Sahade and Struve (9) argue that since both spectra 
of the system have been observed they must be closely similar in type (spectrum 
of brighter star: Bs). Previous light curves, however, have shown very shallow 
and barely detectable secondary minima which lead to estimates of the spectral 
type, of the fainter star, as late as G or early K. Petrie (10), also, has determined 
the difference in visual magnitude, between the two stars, as 1-83 +0"-5 by 
his spectrophotometric method, implying a late type for the secondary component. 
The curve used in the present investigation is notable for its deep and well-defined 
secondary minimum. ‘The significance of this fact will be discussed later. 

Analyses of light curves of RW Monocerotis have also been made by Shapley 
(11) and Fetlaar (12). ‘The former is based on observations by Haynes and 
Seares (13), and the latter on the visual observations of Nijland. ‘There is some 
divergence between the two sets of elements, but no greater than would be 
expected. 

RW Geminorum has also been investigated by Shapley (14) and Fetlaar (12). 
\ preliminary analysis has also been made by Pierce, and has been published 
posthumously by Wood (15), together with the observations used here. ‘The 
various solutions are in tolerable agreement. 

Only two analyses of light curves of ‘TY Pegasi have been published. ‘The 
first is by Dugan (15), which is based on the observations used here; the other 
by Dugan and Wright (16), based on observations made to investigate the 
variation of the period of this system. ‘here is some divergence between the 
two sets of elements. Dugan, in the first of these analyses, met with some 
difficulty in fitting the calculated light curve to the observations, near the shoulders 
of primary minimum—the computed curve falls systematically below the observed 
curve. (‘The analysis was made on the assumption of zero darkening for the 
primary component.) No such difficulty is reported by Dugan and Wright. 

The observations used in the present analyses were obtained with the visual 
polarizing photometer of the Princeton Observatory by the late Drs R. S. Dugan 
(T'Y Pegasi, ( Coronae Borealis) and N. L. Pierce (RW Geminorum, 
RW Monocerotis), and published posthumously by Wood (15). ‘The analysis 
followed the iterative least-squares procedure described by Kopal in Sections 
3.5-3-9 and 3.13 of “The Computation of Elements of Eclipsing Binary 
Systems’”’ (17) for the totally eclipsing systems, and Sections 3.14-3.16 and 
3.18 for TY Pegasi. Owing to the comparatively limited accuracy of visual 
observations no attempt was made to go beyond the “intermediary” solutions 
for the elements of the several systems. ‘lhe errors were calculated according 
to Section A4. Intrinsic weights were evaluated according to Kopal’s 
equation 3.344 for RW Monocerotis and RW Geminorum, and equation 3.376 
for U Coronae Borealis and equation 3.63 for TTY Pegasi. For TY Pegasi only, 
were empirical weights included in the solution. For RW Geminorum, 
rectification of the light curve for ellipticity and reflection effects was found 
necessary, and Kopal’s formula 6.7 was used, assuming B =o. 

For three of the systems, the assumption of a total eclipse at primary minimum 


converged to sets of elements representing the observations satisfactorily. For 
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the fourth—TY Pegasi—however, such a solution did not appear to converge 
to any specific value of the ratio of the radii, k, and a solution was eventually 
obtained on the hypothesis of a partial eclipse caused by an occultation of the 
early type component at the time of primary minimum. 

‘I'wo of the light curves show unusual features. ‘I'Y Pegasi displays what 
Dugan described as ‘‘ an embarrassing series of observations from phase 2¢-6876 
to 24-8022’. ‘These points, which are just before the beginning of secondary 
minimum, fall consistently below maximum light as determined by the remainder 
of the points between minima. As the difference between the two values of 
maximum light obtained by including and excluding these observations in the 
weighted mean amounts only to o™-o1, they were included, although naturally 
this increases the error of the determination. It might also be remarked here 
that although Dugan rectified the light curve before analysing the primary 
minimum, such a course does not seem to be jutsified by the relatively low 
accuracy of observations between minima. 

Irregularities of an even more serious nature appear in the light curve of 
U Coronae Borealis. About three-quarters of the way through primary 
minimum, near phase + 0-15, there is a discontinuity in the slope of the curve. 
On either side of the secondary minimum there appear “spurious minima’”’ 
of much the same depth as the real minimum, but of shorter duration. ‘These 
irregularities look like nothing so much as a wholesale miscalculation of phases, 
but they must be considered in conjunction with others found between minima, 
especially in the latter half of the period. Here, the observations fall into two 
or three well-defined groups, according to phase. Not only is there an 
exceptionally wide scatter in the magnitude of the individual observations—as 
there is throughout the cycle--but the mean magnitudes of the groups are also 
widely scattered, and not in a way that can be explained by ellipticity and 
refleetion effects. Since neither the author nor the editor added any explanatory 
remarks, it is probable that neither of them realized the full extent, or even the 
existence, of these anomalous observations. Under the circumstances, we have 
no option but to omit these recalcitrant observations from our discussion. 
Inevitably this gives rise to some degree of arbitrariness, particularly between 


minima, but the general trend of observations is not violated. Accordingly, 


maximum light was determined from observations 1—-16, 31-48, 62-127 and 
144-155. Observations 17-20, 51-61, 134-143, 176-183, 346-387, 407-481, 
530-613 and 686-703 were used in the analysis of the primary eclipse, while 
secondary minimum, as in the other cases, was estimated directly from the light 
curve 

3. The results and their comparison with previous investigations.—'Vhe results 
obtained from the analysis of the light curves alone are presented in ‘lable I. 
All errors quoted in this table are mean errors. ‘The suffix a refers to the primary 
component: 1.e, the star of greater surface intensity. ‘The sufhx 6 refers to the 
secondary component. ‘lhe symbol r denotes the stars’ radii expressed in 
units of their separation; k denotes the ratio of the two radii r,/r,; ¢ is the 
inclination of the plane of the orbit to the line of sight; p, 1s the maximum 
geometrical depth of eclipse during mid-primary minimum (becoming 1 for 
total eclipses). ‘lhe ratio /,/.J, is that of the mean surface intensities, L is the 
fractional luminosity of the star, in terms of the system’s maximum luminosity 
as unit. Spectral types in brackets have been computed, 6, @” are, 
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TaBie | 


U Coronae 
Borealis 
233+0 O'175+0°014 to-o14 
303 + 0°00" 0°274+0°012 + 0'004 
@: 
° 


RW Monocerotis RW Geminorum TY Pegasi 


0°203+ 0°0I10 
0°311 + 0°007 
0°514+ 0°007 530+ 0°449 + 0005 37 t 0'O13 
0°65 +0°05 a 0°64 +0°096 t+ 0°06 
Sh 1t1 oo 7tt ' 84°°5t1°'1 t 


° 
O° 
° 
°° 


o''4 
I 1 . + 0006 
17°9+3 20°1 + 1°3 44t13 6+4 
2. Sot 1°2 11°47 2° 37°7t073 
Primary Spectrum Ao Bs Bs 
Secondary Spectrum (gG4) Fs (As) 


> 
“~ 


(gG1) 
0°904 + o'016 
Ly o'117 + 0°002 O'134t0 0°350+ 0°005 0'096 + 0016 

j 31 ‘oto 4 a3°'°3 70"° 26°'2+0°''5 
D o4-328+ 0004 of -si4t+o 01-503 +0°010 
3 ‘3+0'5 3 6+0'°6 1'*5+o'4 
0-940 + 0'006 0-057 + 0°010 . 
Period (assumed) _1"-go6102 2"-8654970 


| 0°383 + 0-002 0'366+0 0°644 + 0°005 


25 970 8 
d 


0" "445 + O'OI2 


310922335 


respectively, the phases of external and internal contact. D is the duration of 
the primary minimum, and d that of totality. P is the period of the system 
The two values given for /,/.J, are derived from the formulae 


a 


(1) 


Y(k, - 1) + —* (2) 
h I Ap 
where A, , are the depths of primary and secondary minima, respectively, and 
\(k, —1) is a slowly varying function of k. In each case, it has been assumed 
that the limb darkening coefhicient of the primary component is 0-5, and of the 
secondary component 0-6. 

The results of ‘Table I are not widely divergent from most of the previous 
results cited above. ‘The most notable exceptions to this rule are Fetlaar’s 
elements for U Coronae Borealis and RW Monocerotis and these were derived 
from light curves based on visual estimates. ‘The present analysis of the light 
curve of ‘l'Y Pegasi did not encounter Dugan’s difficulty in fitting the curve to 
the observations at the shoulders of primary minimum. It is true that the three 
(O—C) residuals in this region are of the same sign, but they are all small, and 
can hardly be said to represent any systematic trend. 

In the analyses of the light curves of RW Monocerotis and RW Geminorum, 
however, an outstanding difficulty was the initial estimation of the ratio of the 
radii. For totally eclipsing systems, Kopal (17) gives two ways of determining 
this quantity. It may be estimated directly from the observed depths of both 
minima of the light curve, or determined from the shape of the light curve during 
partial stages of the eclipse in the course of the iterative solution for the elements. 
The two values so obtained are usually called the “depth k” and ‘shape k”’ 
respectively. For good light curves with deep minima, these two values of k 
should agree within the limits of observational error, For the two systems 


3Y 
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named above, however, the depths of the two minima lead to “ depth”’ k’s greater 
than one—values which are ruled out by the light curve. ‘The “‘shape”’ k’s are 
given in ‘Table 1. For U Coronae Borealis a smaller discrepancy exists which 
is in the opposite sense. ‘The relatively large error of the “shape” k for 
U Coronae Borealis, however, renders the significance of the discrepancy in this 
system rather doubtful. It is noteworthy that Miss Gordon (18) has recently 
encountered the same difficulty in her analysis of a much more accurate light 
curve of R'T Andromedae, which is another system of this type. ‘The present 
writer hopes to demonstrate, shortly, that this discrepancy in the two values of k 
is a consequence of the extreme distortion of, and peculiar light distribution 
over, the secondary components. 

A particular effect of the discrepancy is the impossibility of any precise 
determination of the limb darkening coefficients of the primary components 
of the four systems, for it is known that a determination of the “ shape”’ k depends 
crucially on the coefficient of limb darkening of the smaller component (which 
undergoes eclipse during primary minima), while the “ depth”’ k is completely 
independent of it. In consequence, an equalization of the “‘ shape” and “ depth”’ 
k’s for well-behaved eclipsing systems should ordinarily specify the limb- 
darkening of the smaller star within narrow limits. ‘The impossibility of 
equalizing the two k’s, therefore, eliminates this possibility of determining the 
limb-darkening of the primary component. ‘These coefficients had, therefore, to 
be estimated on the basis of the spectral types of the respective stars for the 
effective wave-length of visual observations (A5300). ‘The values adopted, for 
all four systems were u,,=0°5 and u, =0-6 (19) 

A further result of the discrepancy is the divergence between the values 
of J,/J,, as computed from the two formulae. ‘The first of these involves the 
‘shape’ k explicitly, and the second involves the ‘“‘depth’’ k implicitly. For 
RW Geminorum, consideration of the reflection effect provided an independent 
determination J ,,/J,=10°9 + 0°8, and this is the value used in determining L, 
and L, and the spectral type of the secondary component. For the other 
systems, the value (1) corresponds to the quoted values of L,, L, and the 
secondary spectrum 

The value of /,//, quoted for U Coronae Borealis is of some interest. It 
corresponds to a difference in the visual magnitude of the two stars of about 
o™-65 and this is significantly less than Petrie’s value of 1°83 +0°-05 quoted 
earlier. If Sahade’s and Struve’s estimate of the secondary spectrum (g) is 
right, however, the ratio J,/./, should be even smaller. ‘This light curve of 
Dugan’s is the only one with such a deep secondary minimum, and all previous 
observers have found the secondary minimum to be barely noticeable. ‘There 
seem to be strong arguments in favour of both hypotheses about the depth of 
secondary minimum, and further observations are needed to settle the question. 

4. Combination with spectroscopic observations and discussion.—Spectro- 
scopic observations are, at present, only available for U Coronae Borealis, 
RW Geminorum and ‘TY Pegasi.* For U Coronae Borealis, the spectra of 
both components have been observed, three determinations of the orbit and 
absolute dimensions have been made, by Plaskett (20), Pearce (21), and Sahade 


* According to the latest Finding List, RW Monocerotis has been observed spectroscopicalls 
by Millman and Hogg at David Dunlap Observatory 


Dr Heard, Director of the observator, 
informed us in June 1956 that 


the spectrum is Ao with very broad H lines, AA 3933, 4481, and 
a few diffuse metal lines ’’, and that the total velocity range is not likely to be less than 170 km se« 
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and Struve (9). The orbit by Struve and Sahade is characterized by an 
appreciable eccentricity (e=0-13). ‘The secondary minima of Dugan’s light 
curve are quite symmetrical with respect to successive primary minima, and the 
system, therefore, provides yet another instance of a system with photometrically 
circular and spectroscopically elliptic elements. ‘The orbit by Pearce has an 
eccentricity about half the value found by Sahade and Struve, and was, therefore, 
chosen as the most suitable with which to combine the photometric elements. 
The combination of the photometric and spectroscopic elements was made 
without taking the orbital eccentricity into account. ‘lable I] contains the 
absolute dimensions of the system; the symbols m, R, a and p denoting the 
masses of the stars, their radii, the semi-major axis of the relative orbit and 
their mean densities respectively. ‘hese quantities are all expressed in terms 
of the appropriate solar unit. ‘The table also gives the absolute bolometric 
magnitude, M, of each star. 


Tasie I 


Coronae 


RW Geminorum TY Pegasi 
Borealis 


Spectroscopic Mass Function ‘082 0°378 + 0008 ‘0038 
(mass-ratio) 
Contact Hypothesis Mass-Ratio ‘45+ 0°01 0°31 + 0°05 


mM, 8S Ss 
mM» "SS 
a 9 


R, 8 


1 


R 6 


b 
M., ‘Q2 41 


Vi, 3 go 
Pa ‘O92 14 
P» ‘O18 O14 


Notes.—The contact hypothesis mass-ratio of RW Mon is 050+ o-o1 

For U Coronae Borealis the spectroscopic mass-ratio is used to derive the other 
quantities in the table. For the other two systems, only the contact hypothesis mass- 
ratio is available 


In addition to the value of the mass-ratio obtained from the spectroscopic 
evidence, a value may also be obtained on the assumption that the secondary 
component is at its Roche limit—in the sense defined in the Introduction, 
\ table facilitating this has been published by Kopal (2), and the two values 
of the mass-ratio are given in ‘lable I]. All mass-ratios given in ‘lable II are 
defined to be less than unity. ‘laking account of the errors involved, an 
inspection of the two values for U Coronae Borealis reveals that the assumption 
that the secondary component is at its Roche limit is quite compatible with the 
observational evidence. 

For RW Geminorum, information is more limited \lthough both spectra 
of the system have been identified (the observations and calculations are in 
agreement about the secondary spectrum of this system), measurements of 
radial velocities are only possible for the early-type spectrum. ‘lhe only existing 
set of elements has been obtained by Struve (22). ‘The spectroscopic observations 


39” 
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yield only the mass-function, but the mass-ratio, and the absolute dimensions 
given in ‘Table Il have been derived on the assumption that the system is a 
semi-detached one, with the secondary at its Roche limit. ‘The outstanding 
characteristic of this system is the very low masses of the two stars, for their 
luminosity. ‘The system is thus seen to be one of the R Canis Majoris type 
characterized by abnormally small masses and gross overluminosity of both 
components. It should be noted that if the secondary component is smaller 
than its Roche limit, the total mass of the system would necessarily be still less 
than the value quoted here. ‘The abnormally low masses of these stars may 
therefore be regarded as established. 

For ‘TY Pegasi, the only available spectroscopic data are unpublished 
observations by Joy (23) which indicate that the semi-amplitude of the primary 
component's velocity curve is not less than 23 km/sec. According to private 
information from Dr Joy, this figure has been derived from only six plates on 
which “ the lines are few and more plates would be necessary to form any definite 
idea of the velocity curve’. ‘This being so, the figure of AK,=23 km/sec 
provides only the lower limit of the mass function. ‘The mass-ratio has again 
been estimated on the assumption that the system is semi-detached, with the 
secondary component at the Roche limit, and the absolute dimensions so 
obtained are consistent with a primary component which lies fairly close to, 
but just above, the main sequence. 

The main interest in these results is the light they throw on the hypothesis 
that the secondary components of these systems are all at or near their Roche 
limits. In only one instance, that of U Coronae Borealis, is a direct and rigorous 
test possible, and this shows that the secondary component, though it may not 
yet have reached the limit, is certainly not far off. In the case of ‘TTY Pegasi, the 


fact that the assumption that the system is semi-detached leads to a main sequence 
primary component is also encouraging, but the limitations of the available 
spectroscopic evidence must be clearly borne in mind. 

‘The absence of any published spectroscopic observations of RW Monocerotis 
preclude a definite decision in this case, but the closely similar system of 
RW Geminorum provides some of the most interesting, though not the most 


certain evidence. As has already been pointed out, the assumption that 
RW Geminorum is a semi-detached system gains plausibility from the resulting 
low masses of the component stars. ‘his system is one in which the difficulty 
of reconciling the “‘shape”’ and ‘‘ depth” k’s was encountered. ‘The ‘‘ shape” k 
is determined from the primary minimum only, but the “‘ depth”’ & is determined 
from observations of both minima, ‘Thus, it follows that if the secondary 
component were distorted in such a way as not materially to alter the depth of 
primary minimum, but so to alter the observed depth of the secondary minimum, 
such a discrepancy as is observed should result. A star at its Roche limit fulfils 
just these conditions. ‘The cross-section presented to the observer at both 
eclipses is circular to within a few per cent. ‘The shape of the star, however, 
is such that the distribution of light over its surface is greatly affected by gravity 
darkening. ‘The apparent disk presented at the time of primary minimum 
should be of fairly uniform intrinsic brightness, and should be affected only by 
normal limb darkening. ‘The disk which itself undergoes eclipse at secondary 
minimum, on the other hand, is so affected by gravity darkening that it should 
appear to have a completely dark centre and to grow brighter towards the limb. 
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Superimposed on this, there will be the effect of a limb darkening law which is 
more rapid than the normal cosine law. ‘The brightest region of such a disk 
should, therefore, be an annulus at some distance from the centre of the disk, 
and in extreme cases the disk eclipsed at secondary minimum might be almost 
completely dark. 

This last fact provides further qualitative support for the hypothesis that 
the secondary component of RW Geminorum is at its Roche limit. It has already 
been mentioned that the light curve of this system required rectification for the 
reflection effect. As is well known, the effect of such a rectification is to augment 
the total observed light of a system, at any given phase, by the difference between 
the light reflected at that phase and the light that is reflected when both 
components are at “‘ first quarter’’ with respect to each other—1.e. at quadratures. 
In particular, the observed depth A 
fractional luminosity of the “dark”’ side of the secondary component, should, 


, of primary minimum, which gives the 
when rectified give the corresponding quantity for the “bright” side. For 
RW Geminorum the result A, (unrectified) > A, (rectified) was obtained. ‘This 
is clearly impossible for any normal spherical, or even spheroidal star, as it 
implies that despite the light from the primary component reflected by the 
secondary component, its ‘‘bright’’ side is still darker than the “‘dark”’ side. 
It is evident, however, that such an etfect is, qualitatively at least, not incompatible 
with the hypothesis that the secondary star is at its Roche limit. 

It is possible, then, that the discrepancy between the “‘ shape’ and“ depth”’ k’s 
is, in itself, evidence that the systems concerned are semi-detached. In the case 
of U Coronae Borealis, which may not yet be at its Roche limit, the discrepancy 
in the two k’s is much less significant. 

The above discussion has been purely qualitative. It has been sufficient, 
however, to make it clear that a satisfactory method of determining the elements 
of semi-detached systems cannot be evolved until a detailed study of the 
distribution of luminous intensity over the surface of their highly deformed 
secondary components has been made. It will then be possible to predict 
quantitatively the effects discussed more generally here. Such a study is 
currently being undertaken in Manchester by the writer, and, indeed, the present 
investigation forms merely a preface to it. 
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Summary 


An instrument is described for the direct recording of astronomical 
spectral profiles. It consists of a prism monochromator whose output is 
recorded photoelectrically as the wave-length transmitted is varied. Com- 
pensation for atmospheric seeing effects is essential when narrow slits are used, 
and is provided by recording, as the compensated profile, the ratio of the 
monochromator output to the undispersed light collected after the entrance 
slit. The performance of this compensation, and the limiting effects of 
shot noise on the accuracy of the instrument, are discussed 

The instrument has been used with a fairly high resolution (about 0-6 A) 
to record stellar absorption line profiles, and with a low resolution, as a 
slitless spectrograph, to measure the relative intensities of emission lines of 
planetary nebulae. 


t. Introduction.— Astronomical spectra are usually investigated photo 
graphically, but there are advantages to be gained by using a photocell to 
measure spectral intensities directly. ‘The reduction of the results is much 
simplified because the photo-emissive effect is linear, the precision ot the 
electrical measurement involved is potentially high, and the larger intensity 
range over which precise measurements can be made is advantageous, particularly 
in the case of line emission spectra. In addition, the relatively high quantum 
efhciency of available photocathodes increases the speed with which a single 
wave-length element can be measured; and, although this is offset by the need 
to measure successive wave-length elements in time sequence, whereas a photo- 
graphic plate records the whole spectrum at once, the photoelectric method 
becomes advantageous when time changes of a limited wave-length region are 
to be investigated. 

Our instrument consists essentially of a monochromator whose output is 
recorded photoelectrically as the wave-length transmitted is varied. Jacquinot 
and Dufour (1) have discussed the optical properties of instruments of this 
kind and have pointed out their flexibility; a single instrument may be converted 
from slow high-resolution to fast low-resolution performance simply by increasing 


the monochromator slit widths, whereas similar conversion of a photograph 
spectrograph requires the reduction of the focal length of the camera lens. 
Numerous examples of such photoelectric spectrophotometers for laboratory use 


have been described, and examples for astronomical use have been reported by 
Hiltner and Code(2), Guérin and Laffineur (3), and Dobronravin and Nikinov (4). 

The principal difficulty encountered when such an instrument is used for 
stellar spectrophotometry is that fluctuations in the amount of light passing the 
entrance slit are superimposed on the recorded spectral profile. ‘These are due 


partly to changes in atmospheric transparency, but mainly to changes in the 
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amount of light stopped by the slit jaws as the stellar image moves about under 
the influence of atmospheric turbulence and imperfect guiding of the telescope. 
If the entrance slit is made so wide that the whole tremor disk is always included, 
the effect of image movement disappears, but with our instrument such a slit 
would yield only a low wave-length resolution. For many problems a narrow 
slit is necessary, and the fluctuations become so large (Fig. 4) that some method 
of compensation is essential. In the present instrument compensation for the 
fluctuations is achieved by allowing a small proportion of the undispersed light 
which has passed the entrance slit to fall on a second photocell. ‘The ratio of 
the outputs of this photocell ar: the one which scans the spectrum is taken 
automatically; fluctuations common to both cells do not appear in the ratio, 
which is recorded as the compensated spectral distribution. A similar method 
of compensation has been described by Hiltner and Code (2), but for compen- 
sation purposes they used dispersed light in a small fixed wave-length range 
on either side of the wave-length range being scanned, whereas we have used 
light in which the whole spectrum is represented. 


2. Description of the instrument. (a) Optical.—Fig. 1 is a diagram and 
Fig. 2 a photograph of the arrangement. ‘The monochromator 1s auto- 
collimating and of aperture f/6. It has two flint glass prisms, one of 51° and one 
of 30 , and a simple doublet lens of 25 cm focal length mounted on a micrometer- 
controlled slide to permit focusing. ‘lhe linear dispersion is about 25 A/mm at 
Hy, and the maximum resolution at this wave-length is set by optical imperfection 
at about o-4A. ‘The entrance and exit slits are both bilaterally adjustable, and 
the plane of the former ts inclined to the optical axis so that light reflected from 
a glass cover slip held in contact with the slit jaws can be used to guide the 
telescope. 


The light passing the exit slit is focused, by means of an optical system 
discussed later, onto the cathode of the signal cell, which is a photomultiplier. 
Provision is made for cooling this cell in order to reduce its dark emission. It 
is enclosed in a brass tube which is surrounded by a brass jacket containing the 
cooling agent, either liquid air or a mixture of solid CO, and acetone, and the 
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The monochromator, with photomultiplier mountings and head amplifiers attached, 
mounted at the Newtonian focus of the 120 cm reflector at Asiago, Italy 
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Fic. 5 Repeated scans of a stellar absorption line (Vega X 4481) showing (a) the comp nsated 
record, (b) the uncompensated signal record, and (c) the comparison record, Slit width 25 p, corre- 
sponding to a resolution of about o°6 1 The time taken was 6 min per scan. 
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whole is lagged with expanded P.V.C. Light is admitted to the cell along a 
tube of Inconel, an alloy of low thermal conductivity, which also serves to attach 
the cell mounting mechanically to the monochromator. Another Inconel tube 
carries the leads from the photomultiplier base connections to a box containing 
the dynode resistor chain. ‘The air surrounding the photomultiplier and its 
connections and resistors is dried by means of silica gel to prevent the condensation 
of moisture which would cause electrical leakage and misting of the cell window. 

‘The comparison cell, also a photomultiplier, is mounted on the side of the 
monochromator and receives light which is inevitably reflected from the first 
prism face. It is not necessary to cool this cell, since it receives much more light 
than the signal cell. 

‘The wave-length passed by the monochromator is set initially by means of 
a micrometer screw which turns the 30° prism. ‘lhe spectrum is scanned by 
traversing it across the fixed exit slit, and two alternative methods of doing this 
are provided. A short scan of fixed range and rate is achieved by tilting a plane 
parallel glass plate in front of the exit slit, by means of a lever, cam and synchronous 
motor. Alternatively, the wave-length-setting micrometer screw can be turned 





Fic. 6 Block diagram of the electrical system The parts to the left of the dotted line are mounted 
on the telese ope 


by means of another synchronous motor; in this case a selection of pulleys and 
motors provides a wide choice of scanning rates, and the scan may be of any 
length up to the whole range of the instrument. Adjustable stops can be set at 
any pair of wave-lengths to provide a reproducible starting point for the scan 
and an automatic stop at the end. 

(6) Electrical.—¥ig. 6 is a block diagram of the electrical equipment. Part 
of this is carried with the monochromator on the telescope, while the remainder, 
of which Fig. 3 is a general view, is mounted on a rack on the floor of the 
observatory. ‘lhe two parts are connected by a cable 20m long. 

The light entering the monochromator is chopped at 30 c/s by a rotating 
sector disk immediately in front of the entrance slit. ‘The alternating output 
voltages of the two photomultipliers are amplified by separate a.c. amplifiers 
of conventional design, with some negative feedback in each stage to improve 
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the linearity of response. Each amplifier is in two parts: the first part is attached 
to the monochromator and the second part is with the main electrical equipment 
on the observatory floor. The level of signal carried by the long connecting 
cable is thus made so high that electrical interference picked up is negligible. 
\ calibrated logarithmic attenuator precedes the second part of each amplifier. 
The signals from the photomultipliers are accompanied by noise fluctuations, 
sary to reduce the bandwidth 
of the measuring system (which is equivalent to increasing its time-constant). 
We achieve this by using phase-sensitive detection (§), which involves multiplying 
the output waveform of each cell by a reference waveform of the same frequency 
and phase as the signal. ‘The signal then emerges as a d.c. component whereas 
other frequencies present as noise fluctuations emerge as difference frequencies ; 
these can be removed to any desired extent by a low-pass filter, which determines 
the effective bandwidth of the system. 

In our instrument the output of each amplifier channel is effectively multiplied 
by a square reference waveform, as a high-speed relay is used to connect each 
amplifier output to earth for alternate half-cycles at the chopping frequency 
of 30 c/s. ‘This reference waveform is locked in phase to the signal by using 
the chopper to control the relay. Light from a small lamp is interrupted by the 
chopper and falls on a low-grade photomultiplier—the reference cell—whose 
output is applied to the grid of a power valve in series with the relay. For 
simplicity we use only one relay, and this switches the two channels in anti-phase. 
‘lo have the d.c. outputs of the channels of the same sign it is therefore necessary 
to have the a.c. outputs of the amplifiers also in anti-phase. ‘This is achieved by 
having one valve stage fewer in the comparison amplifier, which is permissible 
since the comparison cell normally receives much more light than the signal 
cell. 

Simple resistance-capacity filters are used in each channel with a choice 
of values giving time-constants ranging from 1 to 45 sec; it may be shown that 
the corresponding equivalent bandwidths, from which the total noise may be 
calculated, are from | to 1/g0 c/s. The filters in the two channels are carefully 
matched. ‘The use of a square reference wavetorm, which contains all the odd 
harmonics of its fundamental frequency, gives the output stage additional pass- 
bands centred on these frequencies, and it can be shown that in these bands the 
signal-to-noise ratio is lower than in the fundamental ; the amplifiers are therefore 
designed to have a low response to these higher frequencies. 

The ratio of the d.c. outputs of the two channels is taken automatically and 
recorded by means of a modified Leeds and Northrup “‘ Speedomax ”’ recorder (2) 
which is essentially a self-balancing potentiometer. ‘The signal channel output 
is connected to the normal input terminals, and the battery which originally 
supplied a constant potential to the potentiometer wire is replaced by the output 
of the comparison channel. The pen then records the ratio. For convenience 
in testing the equipment a separate 3-pen recorder is used to record the outputs 
of the two channels separately alongside a duplicate of their ratio. A cathode-ray 
oscilloscope monitors the a.c. output of either amplifier and can also be used to 
test the contacts of the high-speed relay. 


The high voltage required by the signal and comparison photomultipliers 
is supplied by a stabilized r.f. power pack. The h.t. supply for the amplifiers 
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is stabilized, and the heaters of the amplifiers and output stages are supplied 
from a rectifier stabilized by a large accumulator. 

‘To facilitate the use of the equipment by a single observer, meters showing 
the positions of the pens on the 3-pen recorder are attached to the monochromator, 
and several of the important electrical controls can be operated remotely. 
A push-button which stops the wave-length traverse and recorder charts 
simultaneously is also provided; this enables the scan to be stopped should an 
interruption occur, and subsequently re-started without showing a break in the 
record. 


3. Compensation technique.—The first tests were carried out with the 
instrument installed at the Newtonian focus of the Smith—Clarke 18 in reflector 
at the Jodrell Bank Experimental Station of Manchester University. ‘The signal 
and comparison photomultipliers were R.C.A. 1P21 tubes, and following the 
usual astronomical practice the optical system was arranged to image the telescope 


objective on the cell cathodes. ‘The compensation with this arrangement was 
unsatisfactory, for reasons which will now be discussed. 

When the telescope objective is viewed through the slit, dark regions with 
their edges parallel to the slit jaws encroach from time to time on the edges of 
the image of the objective. ‘This effect, which we observed visually when the 
telescope was guided on a bright star, is due to the slit jaws acting as Foucault 
knife-edges. ‘The compensation will be imperfect if identical views of the 
objective are not presented to the two cells, or if the distributions of sensitivity 
over the two photocathodes are different. ‘lhe first effect was removed by 
ensuring that any obstructions which were necessary in one optical path were 
copied exactly in the other. ‘The second effect, which is particularly important 
for 1P21 photomultipliers whose photocathodes are obscured by a closely-spaced 
wire grid, was made much less harmful by placing cylindrical lenses just in 
front of the photocells. ‘hese served to compress the illuminated region on 
each cathode to a narrow strip parallel to the slits. With this arrangement the 
compensation was found to be satisfactory. Fig. 5 shows an example, taken 
from the 3-pen recorder, of the uncompensated signal cell output and the 
compensated result, when narrow slits were used. 

Recently the instrument has been transferred to Asiago, Italy, where it is 
being used at the Newtonian focus of the 120 cm reflector of the University of 
Padua. Various mechanical modifications were necessary and we took the 
opportunity to try a different kind of photomultiplier. With R.C.A. cells, used 
as described above, the effectiveness of the compensation was found to be 
critically dependent on the positions of the photocathodes. Laboratory tests 
on E.M.1. ‘Type 6094 photomultipliers, which have unobstructed photocathodes, 
showed that although there were appreciable local variations of sensitivity over 
the cathode surface there were small areas of substantially constant sensitivity 
which were large enough to include displacements of the star image under average 
seeing conditions. We therefore redesigned the optical systems to image the 
entrance slit on the photocathodes. ‘The compensation was as good as before, 
but no better, and the positioning of the photocathodes was found to be just as 
critical. The E.M.lI. cells were, however, retained for other reasons: the ratio 
of dark current to sensitivity for our best E.M.I. cell at room temperature was 
found to be as low as that for our best R.C.A. cell at the temperature of solid CO,. 
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As explained in the next section this means that whereas it was desirable to cool 
the R.C.A. cell when used on all but the brightest stars, cooling the E.M.I. cell 
is not profitable and some inconvenience is avoided. A further advantage of the 
E.M.I. cells is that they are much less microphonic. 

4. Performance of the instrument. (i) Noise tests.—TVhe recorded output of 
4 photomultiplier which is exposed to a constant light flux shows noise fluctuations. 
‘These are due mainly to the shot effect of the photocurrent and the dark current, 
and the ratio of signal to r.m.s. noise is given approximately by 


i 


S/N 


a. or 

[(t,+2 
where i,,= photocurrent at cathode, 1,= equivalent dark current at cathode’, 
and Af = equivalent noise bandwidth of the recording system. In our instrument 
the comparison cell receives much more light than the signal cell, so that S/N 
for the comparison channel is always much higher than S/N for the signal channel. 
‘The above proportionality. may therefore be expected to hold for the compensated 
record if 7, and 7, refer to the signal cell. 


iG. 7 


By laboratory tests with a lamp of variable brightness we have measured S/N 
as a function of 7,, for our instrument with an E.M.1. photomultiplier, which was 
selected for its unusually low dark current, as the signal cell. ‘The relation 
between S/N and V1, is shown as line B in Fig. 7 for a bandwidth of } c/s, 
which corresponds to a time constant of 1 sec. 


S/N x 1/V Af. 


Our tests have also verified that 
Comparison with stellar records has enabled the V a scale 


* ‘The actual dark current may not all originate in cathode emission: by i, we mean that fictitious 
cathode dark current which would produce the observed dark noise 
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to be calibrated in terms of stellar magnitude: the calibration refers to Ao stars, 
at a wave-length near Hy, using 50, slits (giving a resolution of 1-25 A), and on 
a night of average seeing with the 120 cm telescope. 

The dark current of a photomultiplier may be reduced considerably by 
cooling (6), and line A in Fig. 7 shows the ideal performance which the instrument 
would have if the dark current of the E.M_1. cell were removed completely. ‘The 
vertical separation of lines A and B thus represents the maximum improvement 
that could be produced by cooling. ‘The effect of dark current becomes 
increasingly important as the amount of light, and therefore the photocurrent, 
is reduced: at the point where 7,,=1,, the improvement produced by removing 
the dark current would be V 2, and this photocurrent is represented by limit | 
in Fig. 7. Roughly speaking cooling produces a worthwhile improvement to 
the left of limit I. If 7%, at room temperature is very small, 7, at limit I is also 
very small, and it may not be profitable to use the instrument at all in the region 
where cooling would produce an improvement, because an unreasonably long 
time constant would be required to produce a useful accuracy. ‘This we found 
to be the case for our instrument using the E.M.1. cell: limit | occurs at about 
magnitude 6-5 for an Ao star, using the 120 cm telescope, and a resolution of 
1:25 Anear Hy. With a bandwidth of } c/s, corresponding to a time-constant 
of 1 sec, S/N is then about 1-2, so that to make S/N = 15, say, it would be necessary 
to use a time-constant of about 1 minute. For fainter stars or greater accuracy 
a still longer time-constant would be needed, and the time to record a spectral 
profile would be inconveniently long. 

(ii) Linearity. ‘The photometric linearity of the system was tested in the 
laboratory, using the inverse square law. A diffusing screen illuminated by a 
small lamp was placed in front of the slit, and the output of each channel was 
measured with the lamp at various measured distances from the screen. No 
departures from linearity were revealed by this test, which was accurate to about 
I per cent. 

(iii) Seeing compensation.—\he fluctuations in the cathode current of the 
signal cell due to seeing are proportional to 7,,, whereas the shot noise is propor- 
tional to Vi,,. Any seeing fluctuations which remain in the compensated record 
will therefore be larger in relation to shot noise fluctuations for brighter stars. 
The relative contributions of these two sources of fluctuations have been 
investigated by comparing compensated records from stars with records from 
a steady light source. Part of such a test, made with the 18 in telescope on an 
average night, isshownin Fig. g. ‘This shows recordings of several stars of different 
magnitudes, taken without scanning in wave-length, and with a slit corresponding 
to a resolution of about 1A at Hy. ‘The amplifier gain was adjusted to give a 
convenient deflection in each case. Beside each stellar record is the record of 
a test lamp whose brightness was adjusted to give about the same deflection at 
the same amplifier gain, so that the amount of light reaching the signal cell from 
the lamp was the same as that from the star. ‘The lamp records show only noise 
fluctuations inherent in the signal cell photocurrent (the dark current noise is 
insignificant), while the stellar records contain additional fluctuations due to 
imperfect seeing compensation. Only for the brightest stars are the fluctuations 
in the stellar record appreciably greater than those in the corresponding lamp 
record, 
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From similar tests with the 120 cm telescope it is found that on an average 
night with resolution 1-25 A at Hy the residual seeing fluctuations are equal to 
the photocurrent shot noise fluctuations for an Ao star of magnitude 2°5. ‘This 
condition of equality is shown in Fig. 7 as limit I], and curve C represents the 
actual performance of the instrument with this imperfection of compensation 


taken into account. Some of the imperfection may be due to scintillation in 
colour, but we are not satisfied that the effects of non-uniformity of the cell 
cathodes have been completely eliminated. 

(iv) Variation of sensitivity with wave-length.— he variation of the sensitivity 
of the instrument with wave-length is usually unimportant when single absorption 
lines are scanned, but it must be determined when intensities are to be compared 
over a large wave-length range. ‘The instrument in this context includes the 
telescope and the atmosphere. In the case of the nebular work described below 
the wave-length response was calibrated by scanning the continuum of a star 
of known colour temperature at the same altitude as the nebula. 

The useful wave-length range of the instrument is roughly 3950 A to 6000 A. 
‘The fall-off in sensitivity at the long wave-length end is due to the photocathode, 
and at the short wave-length end to absorption by the glass of the prisms. 

5. Some applications of the instrument. (i) Observations of profiles of stellar 
absorption lines..-A programme of observations of time changes in absorption 
line profiles of variable stars was begun at Manchester and is being continued 
at Asiago. Fig. 8 shows a typical record of the Mg tt A 4481 line in the spectrum 
of Algol, taken at Manchester with a resolution of o-6A, which is near the 
maximum resolution of the instrument. ‘lhe record is on the large chart of the 
Speedomax recorder, and several successive scans of the line were superimposed 
by turning the chart back to the same point at the beginning of each scan. ‘The 
noise is predominantly that inherent in the photocurrent. In this respect it is 
interesting to compare this record with Fig. 5 (a), which shows the same line 
at the same resolution in the spectrum of a much brighter star, Vega. ‘here the 
fluctuations are predominantly due to imperfections of compensation. 

It is worth noting that for stars which are so faint that it would take too long 
to record satisfactory line profiles, records suitable for the precise determination 
of equivalent widths may be obtained with wider slits. 

(ii) Measurement of relative line intensities in the spectra of planetary nebulae. 
‘The instrument is particularly suited to the determination of relative intensities 
of emission lines in the integrated light of planetary nebulae, and a programme 
of such measurements is also in progress at Asiago. For this purpose the entrance 
slit is made wide enough to pass the whole of the image of the nebula, and an 
equal or slightly larger exit slit passes each monochromatic image in turn as the 
wave-length drum is turned. ‘lhe resolution of the instrument used in this way 
depends on the size of the nebula. 

Part of a record from the planetary nebula NGC 6543 1s shown in Fig. to. 
\ complete analysis of observations of this nebula will be given elsewhere, but 
it is interesting to note the large intensity range conveniently covered by the 
instrument. ‘The intensity of Nl, for example, is more than 300 times that of 
He tt A4686; to accommodate such differences the calibrated attenuator in the 
amplitier was adjusted between lines to give convenient deflections on the 
recorder, 
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SURFACE PHOTOMETRY OF THE GLOBULAR CLUSTERS 
47 TUCANAE AND OMEGA CENTAURI 
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(Received 1956 March 16) 
Summary 


Photoelectric measurements of surface brightness, colour and integrated 
magnitude have been made for the globular clusters w Cen and 47 Tuc. 
For w Cen the integrated photographic magnitude is 4°25, mean colour 0°68 
and absolute magnitude — 10°45, or 2°0 x 10® suns. For 47 Tuc the corres- 
ponding figures are: P, 4°7; P-V, 0°67; absolute magnitude 10°0, or 
1*3 10® suns. The two clusters do not follow the same law of luminosity 
distribution, and the light profile of 47'Tuc resembles quite closely that 
found by de Vaucouleurs for elliptical galaxies. 47 Tuc is redder at the 
centre than in the wings 


1. In recent years there has been a marked revival of interest in the dynamical 
theory of globular clusters and the closely related one of elliptical galaxies, papers 
on the subject having appeared from Camm (1), Kurth (2), Belzer, Gamow and 
Keller (3), de Vaucouleurs (4), Woolley (§) and others. A primary requirement 
of such theories is that they account satisfactorily for the observed variation of 
surface brightness across a cluster. It is therefore important that data be available 
against which they can be tested. All previous work in this field has been photo- 
graphic, the most recent having been published in 1936 (6). Since that time the 
superiority of photoelectric methods has become clearly established, and there 
can be no doubt of the advantages to be gained by their application to this type 
of problem in particular. A second of the factors which has led to the present 
programme is that 47‘luc and wCen, the two brightest and largest globular 
clusters, and hence the most suitable for this kind of investigation, are accessible 
only to southern observers. Finally these objects are of considerable intrinsic 
importance, and are badly in need of better colours and magnitudes. 

I‘wo difficulties are inherent in this kind of work. First, because a cluster 
consists of a finite number of stellar sources its light distribution is highly 
discontinuous, and to admit any sort of interpretation must be smoothed. While 
the smoothing may be achieved easily enough— photographically, by using extra- 
focal plates, o1 photoelectrically, by scanning through a hole of finite aperture 


it is difficult in practice to decide on just what degree of smoothing to adopt, 


and each cluster tends to become a particular case. Secondly, the cluster is seen 
against a background, or more accurately through a foreground, of field stars, 
night sky emission and general diffuse radiation. ‘The fluctuations in this 
foreground are difficult to take into account, and set a limit to the faintness of the 
outer regions which can be detected, and to the accuracy of the photometry as 
a whole, 
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‘The observations described in this paper were made with a 1P21 multiplier 
and photometer, on loan to one of us from the Lick Observatory. The multiplies 
current was amplified in a d.c. amplifier, made by Mr D. G. Thomas to a design 
of Dr G. E. Kron, and recorded on a Brown recording potentiometer (except for 
Fig. 1(c), which was recorded on an Esterline-Angus). Blue and yellow filters 
were used to transmit spectral ranges corresponding closely to the International P 
and V systems. ‘The yellow filter was 2mm of Schott GGu1, and the blue 
1mm of Schott BG 12 together with 1 mm of Chance OY 10, the function of the 
latter being to exclude the ultraviolet. ‘he photometer was used either at the 
Cassegrain focus of the 30-inch Reynolds reflector, or with the g-inch Oddie 
refractor. ‘The clusters were measured through holes of different diameters, 
and were also scanned along a diameter from west to east, by stopping the telescope 
drive or altering its rate. North south scanning was not practicable 

It is convenient to collect here some of the principal data relating to 47 ‘Tuc and 
wCen, and also to summarize some of the main results of this paper. The 
distances and moduli are from a compilation of Shapley’s (7); both are uncertain, 
especially for 47'Tuc. ‘The distance to the galactic centre has been taken a 
&4kpe, ‘The distance moduli are apparent (i.e. not corrected for absorption) 


‘TABLE | 

$7 ‘Tuc 
Right ascension (1900) 0! 22™ 
Declination (1900) 
Longitude (/) 
Latitude (b) 
Distance from sun 6°8 kp« 
Distance from galactic centre ‘ 6°4 kpe 
Distance from galactic plan 1°7 kp 
Concentration class Vill 
\pparent distance modulus (m — .\7) 14°7 14 
‘Total magnitude (P) 
Absolute magnitude (P) 
Luminosity 
Mean colour (P-l’) 
Central brightness (P mag/sq. sec) 
Effective radius > pe ‘ 10°0 pr 
Maximum radius 2 40’ — 79 pr 


2. Observations of 47 Tuc.—Vhe following observations were made of 47 Tuc 

(1) measurements in blue and yellow of the light received through apertures 
of various diameters, from 41” upwards. ‘These are detailed in ‘lable Il. ‘The 
observations on the g-inch were made on three nights. 

(ii) two-colour measurements on the 30-inch reflector through a 58’ hole 
at about a dozen points, distant up to 12’ from the centre. ‘These observa- 
tions were also made on three nights. 

(111) scans in blue and yellow through a 58” hole on the Reynolds at diurnal 
speed (4'-55 per minute of time). 

(iv) scans in yellow at 02745 diurnal speed (1'-25 per minute of time), made 
on the Reynolds through 16”, 26” and 58” apertures. ‘I'wo of these are repro 
duced in Fig. 1. 

(v) a scan in yellow at diurnal speed through an aperture of g’ diameter on 
the g-inch refractor, 


40 





S.C. B. Gascoigne and E, 7. Burr Vol. 116 


The measurements were reduced to zero air-mass and referred to the P, V 
system by observations on the same night of stars in one of the Harvard Standard 
Regions C12, D2, Eg (8). (In this paper Pand V mean IPg and IPv as defined 
by the Cape Observers), ‘The connection between IPv and IPg and Johnson 
and Morgan's V and B is: 

V=IPv; B=o-g1 IPg+o-0g IPv+o0-16. 
On a few occasions the star 8 Oct was used as intermediary, although because 
of the wide difference between its colour and that of the cluster it was not 
very suitable for this purpose. Our magnitude and colour for this star agree to 
within o-o1 mag with those found at the Cape (9). 

To combine the above data into a light profile for the cluster we have to take 
into account field stars, irregularities and lack of circular symmetry in the cluster 
itself, and distortions introduced by the finite size of the scanning aperture and 
by the response time of the amplifier and recorder. ‘The nature of some of these 
difficulties is illustrated by Fig. 1, which shows scans at 1’-24 per minute through 


Fic. 1 (a) Scan of 47 Tuc made in yellow with a 16" diameter hole at 0°2745 diurnal pee 7(1°+25 
arc per minute); (b), as (a), with a 58” diameter hole, and (c), a scan in yellow at diurnal speed through 
ao’ diameter hole (c) was made with an Esterline—Angus recorder In each record the broken 
horizontal line indicates the sky background, the full line zero tllumination The magmtude scales 
hole diameters and distances from the centre are also indicated 
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holes of 16” and 58” diameter, and at 4'-55 per minute through a g’ hole. The 
smoothing introduced by the larger holes is very noticeable. The distortion 
which arises from the finite response time of the recorder is less obvious 
Particularly in the central region where the brightness is changing rapidly, the 
recorder pen is unable to keep up with the light signal at the scanning speeds 
available. It was decided in consequence to derive the brightness within a radius 
of about 60” from the data in Table II. If the surface brightness in units of 
a tenth magnitude star per square minute of area is denoted by f, and the light 
emitted within a radius r by /(r), we have 


a 


Kr)=an| rf(r).dr. 


‘Table II gives the /(r) for various values of r. ‘hese values were plotted against 
r*, as shown in Fig. 2, and a smooth curve drawn through them and differentiated 
graphically to give f. In particular the brightness at the centre corresponds to 
the slope of the curve at the origin. 


TABLE I] 
Light received from ly Tuc and w Cen through circular apertures of various diameters, ¢ “press / 


in terms of that of a tenth magmtude star Sky corrections have been applied 
Lele scope Diameter | P P-| 


luc 


9 inch refractor 16°2 


sk 


20-inch reflector 


20-inch reflector fy I oO 49 h 


The 47 Tuc observations with the 9-inch are the means of three night Other observatioy 
refer to one mght only The observations on the 20-inch through the 116" hole was made at ti 


Vewtoman focus, the others at the Cassegrain focu 


Unpublished star-counts of 47 Tuc to 16 mag, kindly made available by 
Dr de Vaucouleurs, indicate that the cluster is elliptical, the ratio of the axes being 
0°92 + 0°02, and the position angle of the major axis 37° +1 \s this is also the 
angle between the minor axis and the diameter along which our scans were made, 
it seemed unnecessary to make corrections for ellipticity when deriving the 
integrated magnitude of the cluster. No evidence was found for E-W asymmetry. 

In the outer regions the values for the surface brightness were combined 
from the various scans. It is clear that no one scan can cover more than a limited 
region of the cluster. ‘The irregularities must be almost entirely inherent in the 
cluster. In Fig. 1 (a) for example there is only a ten per cent chance of a field star 


40" 
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brighter than 14 mag occurring in the strip measured (16” wide by 15’ long). 
The profile actually adopted is shown by the dotted lines. The arbitrary nature 
of these lines seems unavoidable. In the case of Fig. 1(c) a correction for the 
finite size of the scanning aperture shows the true curve to lie slightly within the 


00! T T T T T 
on _ 


Ya 


80} 


ay 
Oliv | | 
0 6 10 
miné 


Fic. 2 The light received from 47 Tuc (full line) and w Cen (broken line) through circular aper- 


tures of various diameter The light, expressed in terms of that of a tenth visual magnitude star, 1 


plotted against the area of the hole in square minute 


observed curve (10). According to this scan the brightness of the cluster drops 


off very rapidly at the edge, and it cannot be traced photoelectrically beyond 
a radius of 22’, when its brightness is about five per cent of that of the sky. ‘This 
radius is almost certainly an underestimate. Shapley’s value, derived from 
microphotometer tracings, is 27’, and de Vaucouleurs’ star counts indicate a radius 
of at least 40’ and possibly 60’. ‘The sky brightness was found to be 
1 3°2 mag/sq min, or 22°1 mag/sq sec (pg); this is reasonable for a latitude of — 45 

Information on the colour of the cluster and on its variation with radius can 
be extracted from (1), (ii), and (111) above. ‘Thus from the data in ‘Table II, we 
find the colours in successive annuli to be: 


" | 
Annulu oO -20 R14 315 -50 


Colour 0°77 74 o'60 o’6s , 0°O55 


I‘hese results are plotted in Fig. 3, together with the colours derived from measure- 
ments at various points (ii), and from two-colour scans (iii). Some of the plotted 
points therefore refer to whole zones, and others only to small regions. Even 
o the agreement is not good. ‘The annular measurements are the means of 
three nights and should be reasonably free from observational error. ‘The 
effect of field stars should be negligible. 

The major effect here is probably that the hole used (area 0-72sq min) was 
too small to take a good sample of the cluster. Sandage (11) found that 0-20 
of the light of M3 came from stars with M, brighter that —1-8. If this is true 
for 47 Tuc, a fifth of the light will come from stars of apparent magnitude brighter 
than 12°5. 5° from the centre the surface brightness is equivalent to 4:1 stars 
of mag 12°5 per square minute (‘Table III). ‘Thus there is about a sixty per cent 
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chance of a star brighter than 12:5 mag appearing in the hole used here, while 
the total light measured through the hole is only that from about three such stars. 
The corresponding fluctuation in colour could easily be a tenth of a magnitude. 
A corollary to this is that unless holes of large angular size can be used, star counts 
will give a more satisfactory picture than will a photometer of the outer structure 
of globular clusters. 


FIG. 3 Variation of colour with radius in 47 Luc Filled circles, mear fannuli; openci 


measurements at points along the diameter of position angle go The hor mtal lines indicat 


widths of the annuli. The broken line was derived from two-colour scans, and the full lin that 
adopte d 


But in spite of the large scatter there seems little doubt that 47 Tuc is blue 
in the wings than inthe centre. ‘This effect had been foreshadowed by Stebbins, 
who in the course of a description of some photoelectric measurements of globular 
clusters said: ‘‘ Probably most of these clusters have varying integrated colour 
from centre to outside. We have found as much as o-17 difference on th 
V—I scale, (0:06 on the P-V scale) inner redder than the outer region, but a study 
of this difference will involve a check upon the possible instrumental effects ”’ 
(12). In the present work the unexpected feature is therefore not so much the 
existence of the effect as its magnitude. 

Colour differences of this sort may be significant. Animmediate interpretation 
is that the blue stars in globular clusters are to be found farther from the centre, 
and hence are less massive than the red stars; this is at variance with current 
ideas, which suppose that almost all the light of globular clusters comes from 
stars of substantially the same mass. 

Final results are collected in ‘Table III. Columns 4 and 5 show ((r), the 
amount of light within a radiusr. It will be shown later in the paper that, except 
near the centre of the cluster, the logarithm of the surface brightness decreases 
very nearly asr™*, If we assume that this law continues to hold, the contribution 
to the integrated magnitude beyond the observable region (22’) is o-o5 mag; 
this figure is probably an upper limit. ‘The integrated visual magnitude of 
47 Tuc is then 4-01 + 0°05 say, and the mean colour is 0°67. ‘This gives a mean 
photographic magnitude of 4-68: Shapley’s estimate was 4°5 (7). Half the light 
of the cluster is found within a radius of 2'-7, or 6-0 parsecs at the distance assumed 


in ‘Table I. 
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‘TABLE II] 


Photometric measurements of 47 Tucanae 
(3) (4) (5 (6) 


03534 


0417 
0243 


0099 


0050 027 
00297 ‘016 
oo1si 0095 
OOlI2 006! 

2*0007 | 0034 
00044 "0024 

"0003 OOS 


0°007 
sky 0°09 
} ) 
( olumms (1) and (3) are distances in angular measure and in parsecs., olumn (2) 18 t/ 


ivface brightness in tenth visual magnitude stars per min* Column (3) is the P-V colour, and 


the total light inside a radius r is shown in column (4) as a number of tenth magnitude stars, and 
im column (5) in stellar magnitudes (both visual) 
nitude stars per « 


Column (6) is the space density in tenth may 


ubte minute and column (7) the space density in stars of zero absolute magmtude 
per culne parsec 


) 


3. Observations of w Cen.—The following observations were made of w Cen: 
(1) measurements in blue and yellow of the light received through apertures 
of various diameters. ‘These are detailed in ‘Table II. 
scans in yellow through a 140” diameter hole on the 30-inch reflector at 
a speed of 10-2 arc per minute of time 
scans in yellow through a 26” diameter hole on the 30-inch ¢eflector at 
a speed of 2'-81 arc per minute of time (one of these is reproduced in 
Fig. 4). 

‘The measurements were reduced to zero air-mass as before and related to the 

P, V system by observations of stars in E6 (8). 
The reductions proceeded much as before. ‘he surface brightness of 
w Cen varies less rapidly but is rather more irregular than that of 47 Tuc. As 
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with 47 Tuc the bumps on the wings are almost certainly inherent, and the profile 
was drawn through them, as shown in Fig. 4. The dip in the centre is real, and 
must be due either to a deficiency of cluster stars or to a highly localized patch of 
absorbing matter (cf. Plate 2 of reference (13)). Such irregularities are not 
uncommon in globular clusters. In default of a better procedure it was decided 
to draw a smooth curve through the dip, and to determine the brightness in the 
central regions as before, by graphical differentiation of the data in Table II. 
In this computation it was necessary to take into account the well-known departure 


Fic. 4 1 scan of w Cenin yellow through a 26" diameter hole at 0°2745 diurnal speed (2'*%1 a) 
per minute of time). The full horizontal line indicates zero illumination, the broken line the sh 
background. The dotted line is the adopted curwe 


from circular symmetry of wCen. From a study of Schilt’s photographi 
measures (14) and Shapley’s star counts (15), it was decided that the isophotes 
could be represented well enough by a series of concentric ellipses with major 


axes all directed east and west. Combining as much data as possible to smooth 
out local irregularities, the ratio of the axes b/a was estimated as a function of the 
semi-major axis a from Schilt’s Fig. 1 (a<g’) and from Shapley’s star counts 
(4°°5<a<14’). ‘The resulting points and the adopted curve are shown in Fig. 5 
The curve has been extrapolated on the assumption that the isophotes are circula: 
when a ts very small or very large. 

The results are collected in ‘Table IV. As has been pointed out, the figures 
for the innermost minute represent a considerable degree of smoothing. ‘The 
colours of the annuli, derived from Table II, are as follows 


Annulu Oo-100 


Colour #) 76 0'60 
These figures are again much more irregular than would be expected, but 


as they stand hardly suggest any systematic colour variation in wCen. More 
observations are needed to clear up this point. 
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Lillipticit wv 1 The filled circles are from Schilt urface photometi 
from Shapley’s star counts 


‘ABLE I\ 


jtometric measurements of a Centaur 


(4) (6) 


I 


' 
28 


/ 
ha 


73 
51 
"38 244 y ‘0266 
29 
221 272 O145 
160 
116 , a 87 0070 
00S 102 ae ( 0030 
044 312 2 00154 
034 
0265 00074 
0167 00043 
oll 00024 
32 007 |! ‘00016 
36 004! folele) fe) 
40 002 363 360 
sky 20 
Columns (1) and (8) are radti in angular measure and in parsecs. Column (2) ts the surface 
brightness in tenth visual magnitude stars per min*. The light inside a radius r its shown in 
column (3) as the equivalent number of tenth magnitude stars, and in column (4) in stellar 
magmtudes. Column (5) ts the ratio of the axes, column (6) the space density in tenth magnitude 


stars per cubic minute, and column (7) the space density in zero absolute magnitude stars per 
cubre parsec 


40’ from the centre the brightness of the cluster has fallen to about one per cent 
of that of the sky. ‘The resulting diameter of 80’ is rather larger than the 65'-4 
found by Shapley (7). ‘The sky was found to be about 0-8 mag brighter than in 
the vicinity of 47'Tuc. It is to be emphasized that all our sky measures are the 
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means of a number of areas, chosen to avoid either abnormally dark patches or 
conspicuously bright field stars. Out to 40’ the integrated visual magnitude is 
3°60. ‘To extrapolate beyond this point we assume that the logarithm of the 
surface brightness continues to decrease as r'4 (cf. Section 4). ‘The contribution 
from the unobserved part is then 0-03 mag and the total 3:57. ‘The mean colour is 
0°68. ‘This gives a mean photographic magnitude of 4:25, appreciably brighter 
than Shapley’s 4-7 (:). Half the light is concentrated within a radius of 5'-05, 
or 1O parsecs. 

We note here that half the variables are contained within a radius of 67 arc. 
‘rom the open nature of this cluster it seems unlikely that many can have been 
missed, even in the centre. ‘This may be an indication that the variables ar 
less massive than the stars which produce most of the light of the cluster. 

Schilt’s measurements have been referred to above. After an adjustment has 
been made to the zero of his calibrating stars, his figures agree with ours to about 
0-05 mag in the mean. In view of the unsatisfactory magnitudes Schilt had to 
use for these stars, this agreement is as close as could be expected. 

4. Discussion.—It has been supposed for many years that globular clusters 
obey a common law of luminosity distribution, although the only real evidence 
for this statement seems to be a photographic investigation by Lohmann (6), 
who, after allowing for scale factors, found substantially identical distributions 
for the light in the clusters, Ms5, Mis and Mgz2. Star counts have also been 
used as basis for this comparison, and recently a very good series has been published 
for Mg2 by ‘layler (16). Our results are compared with Lohmann’s and ‘Tayler’: 
in Fig. 6, where we plot log f/f, against r/r,, f, being the brightness (or star count) 


Fic. 6 Relative brightness as a function of relatwe radius for a nonber of globular cluster 
The full line is 47 Tuc, the broken line w Cen The circles are Lohmann photographic results from 
U5, Mis, and Moz, the crosses Tayler’s star counts for Mg2 
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at the effective radius r,, i.e. that which contains half the light of the cluster. 
From Lohmann’s work the effective radii for Ms, Mis and Mg2 were found to 
be 100”, 58” and 75” respectively. 

In the figure the full line is 47 Tuc, the broken line w Cen, the points Lohmann’s 
(from three clusters) and the crosses Tayler’s. It is clear that 47 Tuc and w Cen 
differ significantly both from each other and from Lohmann’s and ‘Tayler’s 
clusters. On the other hand Lohmann’s and ‘layler’s data confirm each other 
well, as they should, for the stars ‘Tayler counted supply almost all the light of the 
cluster, and he found their radial distribution to be largely independent of 
magnitude. 

Three comments may be made. First both 47 ‘luc and w Cen are abnormal 
in being at least four times brighter intrinsically than any other globular cluster. 
Secondly, there is a steady gradation in concentration class from III for 47 Tuc 
through IV and V for Lohmann’s clusters, to VIII for w Cen: as will be shown 
shortly the reality of this gradation is supported by measurements of absolute 
surface brightness. Finally our method of smoothing the centre of w Cen may 
not be free from objection. However, the differences in Fig. 6 seem well in 
excess of any observational error. 


| 1 


1 
r | 


Log surface brightne avainst the fourth root of the radi for 47 Tu 
(dots) 


¢ (crosses) and w Cen 


\n interesting comparison can be made between the brightness distributions 
in 47 Tuc, w Cen and the elliptical galaxies. Of the various empirical laws which 
have been proposed for the latter we select that of de Vaucouleurs (4). ‘This is 


log f/f, 3°33i(r/r,)h* 1 


f, and r, defined as above.* In Fig. 7 we plot log f for the two clusters against 
r!4 yin minutes. Over a brightness range of 3000 the 47 Tuc points lie within 
+ 0°03 of a straight line, and the deviations from the line are in the same sense as 


de Vaucouleurs (private communication) considers the figure preferable to the 3°25 


vhich he original) suggested 
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those of elliptical galaxies (cf. Fig. 2 of (4)). But according to de Vaucouleurs 
the central brightnesses of elliptical galaxies is about 160 times the brightness 


at the effective radius, whereas the corresponding figure for 47 ‘T'uc is only about 4o. 


Thew Cen observations also fit a straight line, over a brightness range of about 
1000, though the deviations are greater than with 47 Tuc. 


For distances from 
the centre greater than about 3’ the two clusters follow the same brightness 


law to within better than 20 per cent, and differ only in that w Cen is relatively 
much fainter in the centre. ‘That the two lines have the same slope is of course 
only a coincidence. With an elliptical galaxy this slope would correspond to an 
effective radius of 2'-61; with the clusters the true effective radius is larger, 
because of the lower central intensity. 

The absolute central brightness of 47 luc is considerably in excess of that 
of other globular clusters. The following table, which includes some figures for 
elliptical galaxies, illustrates the point: 


‘TABLE \V 
Central brightnesses of some globular clusters and galaxies, in nu f tenth photograph 
magnitude stars per square minute 
Globular Clusters 
Cluster 47 Tuc Mis VI 92 w Cen 
(Concentration class It] LV 1\ Vill 
Central brightness 30 8°5 6°55 d 3°4 


Elliptical Galaxies 
NGC numbers 


221 3115 3379 4494 4049 
‘T'ype 


k2 E47 Eo ko I 
Central brightness 100 12 21 10 


/ 


Note, inter alia, the tendency in the globular clusters for central brightness to var 


regularly with concentration class 


We derived ¢, the light emitted per unit volume, from the equation 


1% —f'() 
ory=-+{ LO) ax, 
aw), (x*—7r*)\¥ 
where / is the surface brightness and R the cluster boundary (17). ‘10 avoid the 
infinity in the integrand we replaced x* by r? + 2; the integral then becomes 
1 (cf 
d(r) = g(V 98 + 2*) dz. 

We have written g(r) for f’(r)/r, and the upper limit (R* — r*)'* is denoted by 7 
lhe function g, which is always finite, was found by numerical differentiation 
of log f, and plotted on semi-log paper, the requisite interpolations being made 
graphically. ‘This scheme, which worked well, seems 
convenient than others which have been proposed (18,19). ‘The contribution 


to ¢ for z larger than about 2r were small, and for practi al purposes the upper 
limit could be taken as infinite. 


simpler and more 


For points near the observed boundary it was 
necessary to extrapolate. For 47 luc we did this by assuming that the surface 
brightness continued to fall according to de Vaucouleurs’ law; with w Cen 
the experimental law log f = — 0-62-0°048r was better for large » lhe values of 
¢ for these points may therefore be rather uncertain 

The space densities are shown in column 6 of ‘lable [II for 47 ‘luc, and 


in the same column in Table IV for w Cen, in units of a tenth visual magnitude 
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star per cubic minute. In column 7 in each table we have converted these to 
zero absolute visual magnitude stars per cubic parsec. This step requires 
a knowledge of the cluster distances, but this is not critical, a half-magnitude error 
in the modulus producing a change of only twenty per cent in the conversion 
factor. Averaging over stars within ten parsecs we find that in the region of the 
Sun the light emission is equivalent to that from 0-0005 stars of zero absolute 
magnitude per cubic parsec. That at the centre of w Cen is about 30000 times, 
and at the centre of 47'Tuc about a million times this figure. ‘To an observer 
at the centre of 47’Tuc the light of the stars would be equivalent to that from 
several thousand full moons, and the resulting illumination quite comparable 
with terrestrial twilight. 
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A NOTE ON THE SPECTRA 
OF SOME VARIABLE STARS IN THE MAGELLANIC CLOUDS 


M. W. Feast 


(Communicated by the Radcliffe Observer 
(Received 1956 June 7 


Summary 


Spectra (86 A/mm at Hy) have been obtained of a number of variables 
in both Magellanic Clouds. Four long period cepheids in the Small Cloud 
and six in the Large Cloud show normal supergiant spectra of late F and 
early G type near maximum. ‘There appears to be no spectroscopic evidence 
for their being other than normal Population [I cepheids. HV 2882 is a 
foreground Me variable. ‘The R CrB variable W Men is a member of the 


Large Cloud and has an absolute magnitude at maximum of sig 


The long period cepheids.—Vhis note gives the results of a spectros¢ opic study 


of some variable stars in the two Magellanic Clouds. ‘The majority of star 
considered are long period cepheids. ‘These stars are of interest since it has 
been suggested that they may differ significantly from galactic Population | 
(classical) cepheids. Such a conclusion would be of far-reaching significance: 
both for the interpretation of the Magellanic Clouds and also for all determinations 
of extragalactic distances based on the cepheids. ‘The position has been 
summarized by Gascoigne (1). ‘The discussion of the problem has up till now 
rested mainly on the results of photoelectric observations of Magellanic Cloud 
and galactic cepheids. ‘The results given in Gascoigne’s Fig. 15 show the cloud 
variables to be apparently bluer than the galactic ones. ‘This led Gascoigne 
and Kron (2) to suggest that the cloud variables were in fact Population II 
(W Virginis) cepheids. Such a suggestion appears to be ruled out by the 
change in the distance of the Clouds (3) which indicates that the cloud variables 
have about the same luminosity as galactic Population | cepheids of similar 
periods. ‘The discrepancy in the colours remains however and its interpretation 
as due to the underestimation of the reddening of galactic cepheids is not entirely 
favoured by Gascoigne (1). It is clear that spectroscopic work on the cloud 
variables might help to clarify the position. 

The results on some Harvard variables classified as cepheids are given in 
Table I. There are four variables in the Small Cloud and six in the Large Cloud 
as well as one star which turns out to be a foreground object. ‘The variables are 
arranged in order of increasing period and an approximate phase of the 
observation is given if available. ‘The latter depend mainly on unpublished 
epochs kindly supplied by Dr Shapley and Mrs Nail. ‘The possibility that 
some of these calculated phases may be in error must be borne in mind in 
discussing the observations, especially in view of Gaposchkin’s work on HV 2447 
(4) which indicates that the long period cepheids may be somewhat less regular 
in their behaviour than those of shorter period, ‘The spectra from which th 
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types were estimated were all obtained with the f/2 camera (86 A/mm at Hy) of 
the 2-prism spectrograph at the Cassegrain focus of the 74-inch Radcliffe 
reflector. In view of the faintness of the stars the spectra were not widened 
as much as is the usual practice in using the MK system and a few spectra were 
less well exposed than might have been desirable. However the spectral types 
should not be in error to any appreciable extent. ‘The average exposure time for 
the spectra listed was about 270 minutes. ‘The spectra were classified on the 
MK system by comparison with standard stars drawn from the list of Johnson 
and Morgan (5) 


TABLE I 


Cepheid Variables in the Magellamec Clouds 


Small 
or 
H\ Large Period Phase Spect- Magnitude Colour 
Cloud rum P 


LMC ‘ 1956 Jan F8 I | 150] 
LM¢ $3 1955 Dec 3 7 I(a) 6-15 °3 
SM¢ ; 1985 Dec Go | ‘6—-14'2] 
L™M¢ S°2 1956 Jan 2 0°35 late | 2 14°2| 
early G 
SMC - J 1955 Sept. F8 la \ an 
1955 Oct. } F8 Ia f 
SM¢ “1955 Aug. Go I 12°O:-13°4 
LMC 1956 Jan. 52 G21 [12°3-13°6| 
1. Me 1956 Mar ‘34 Go I 12°8-14°3 
1953 Dec. ‘QI Go I 
1955 Oct Gril { 
883 LMC 1956 Jan °22 G2 Il 12°5-14°3 
882 Fore- 2 1956 Jan. 64 Me [12°3-15"1] 
ground 


821 SM¢ 12°4-13°5 


(1) It appears doubtful if the variable was as faint as would be indicated by the calcu- 
lated phase 

(2) K7e (HA100), Mb (HB 754) HDE 270419. 

(3) Ks (HAroo) HDE 269362 


Photoelectric magnitudes and colours from Gascoigne (1). Those in brackets are 
Harvard photographic 


I'he radial velocities of the stars were measured in a Hartmann spectro- 
comparator against the supergiant Go star 8 Aqr and were all found to agree 
well with the expected velocities for cloud members. In cases where more than 
one spectrum was available the agreement of the velocities is good. It is not 
considered desirable to publish the individual velocities at this stage since furthe: 
work is contemplated on these and other variables. The velocities averag: 

133 km/sec for the Small Cloud and +274 km/sec for the Large Cloud 
HV 2882 has a low velocity indicating that it is a galactic, foreground, object. 
(his is confirmed by the spectrum which shows TiO absorption and H 
emission. It seems to be a normal Me variable. 


The Magellanic Cloud variables appear to be normal late F and early G tyjx 
supergiants. No abnormalities were noticed and in particular no emission lines 
were found which would have indicated them as possibly Population I| 
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(W Virginis) cepheids. It is instructive to compare the present results with 

those found by Code (6) for galactic cepheids. He found that over the range 

of periods studied by him there was no change in spectral type at maximum witl 
the period, Code’s results and those of this note are shown in Fig. tf. 
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Fic. 1 Relationship betiween spectral type at maximum and 


Code (6) (filled circles) and Magellanic Cloud cepheids (open cu 
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It should of course be borne in mind that a spurious dependence of type on 
period could result from the types used not referring exactly to maximum light, 
there being a strong dependence of type at minimum on period. However, the 
dependence of type on period which appears when the longer period stars are 
considered is borne out by the colours of the cloud variables (1) which also show 
a slight dependence on period. 

It would appear that, although the Magellanic Cloud variables studied here 
have exceptionally high luminosities for obvious reasons of observational 
selection, there is no spectroscopic suggestion that they are other than normal 
Population I cepheids. In considering the problem of the difference in colour 
of the cloud and galactic cepheids Code’s spectroscopic observations deserve to 
be given some weight for they indicate that for the range of periods studied by 
him the intrinsic colours should show little dependence on the period whereas 
Eggen’s colours of Code’s stars (7) show a marked increase of reddening with 
period, a phenomenon shown by galactic cepheids in general—see Gascoigne (1), 
Fig. 15. ‘he constant spectral types make a good case for attributing the changes 
in colour to an underestimation of interstellar reddening, the effect being 
dependent on the period since the long period stars are brighter and on the 
iverage more distant. In fact a very recent paper by Stibbs (8) does indicate 
that galactic cepheids suffer much greater space reddening than hitherto supposed. 
It appears that the case for considering the Magellanic Cloud cepheids as normal 
Population I cepheids is now sufficiently strong for the intrinsic colours of the 
cepheids to be obtained from observations of these stars. ‘This would probably 
lead to a somewhat greater space reddening of galactic cepheids than is found 
by Stibbs who used the space absorption deduced by Oort, van Rhijn and 
Parenago from B star colours, but this conclusion does not seem out of the 
question. 

The RCrB variable W Men.—The variable W Men [05" 27:1-71°14 (1950) 
13-8-16™ pg] was discovered by Luyten(g). As will be seen from its coordinates 
it lies in the direction of the Large Magellanic Cloud. Its importance lies in the 
fact that it is an R CrB type variable and if it is actually a member of the Large 
Cloud it is the only star of the class known in an extragalactic system and the 
only one with a known distance and absolute magnitude. The only other means 
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available at present for determining the absolute magnitude of this class of 
variable are spectroscopic criteria. ‘The stars show supergiant spectral 
characteristics but, in view of the considerable abnormalities also found 
(particularly strong C, weak H), little reliance has in the past been placed on 
the spectroscopic parallaxes. 


In view of the faintness of the star it was necessary to use a short trail to 


obtain a spectrum with the f/2 camera, and the spectra are necessarily not of the 
best quality. ‘They are however good enough for certain definite conclusions 
to be reached. ‘The best spectrum was obtained on 1953 December 1 and 
another spectrum (in all respects similar) on 1954 February 10. It had been 
intended to obtain more spectra in 1956 but examination of the field on 
January 18 and March 8 showed that the star had faded out of reach of the 
spectrograph. On both occasions it was probably fainter than 16™ visual. ‘This 
circumstance at least has the advantage that one can be certain of the identification 
of the star observed with the variable, a matter of some concern with so faint a 
star in such a rich field. ‘lhe radial velocity of the star is about +260 km/sec 
in good agreement with the expected velocity for a cloud member in this 
position. ‘The spectral type is necessarily somewhat uncertain; it appears, 
however, to be a supergiant of type F5—Go, most probably F8 with Hy weak 
for the class. ‘The observations therefore suggest that W Men has a typical 
R CrB type spectrum though the spectra are not suitable for determining whether 
or not C and C, are present. We may conclude that the supergiant characteristics 
shown by the R CrB variables are real indicators of high luminosity in these stars. 
Using a distance modulus of 19-2 for the clouds (3) we find an absolute 
magnitude at maximum for W Men of — 5™-4. 

Acknowledgments.—\1 am much indebted to Dr H. Shapley and Mrs V. 
McKibben Nail for kindly supplying unpublished Harvard data and maps on 
some of the variables discussed in this note. I should also like to thank 
Dr A. D. Thackeray for his encouragement and for many helpful suggestions. 
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RED SUPER-SUPERGIANTS IN THE LARGI 
MAGELLANIC CLOUD 


M. W. Feast and A. D. Thackeray 


(Received 1956 December 7) 


Summar\ 


Five members of the Large Cloud have been discovered with absolute 
visual magnitudes about 9°O, spectral types from FS to Gs and with 
colours between 0°37 and 1°45 Four of these stars are among the 30 
brightest known members of the Cloud. ‘The census of the brightest blue 
stars must be regarded as complete, but owing to the difficulty of distinction 
of foreground stars more red super-supergiants probably remain to_ be 
discovered. 


he Magellanic Clouds offer a unique opportunity to study individually the 
brightest stars in an external system. ‘The high degree of resolution guarantees 
in most cases that the slit of a spectrograph will isolate the light from single stars. 
In apparent magnitude there is a gain of 5 magnitudes over corresponding objects 


in M31. Consequently the number of objects which can be reached spectro- 


scopically with existing equipment is very large. If one adopts 13°6mag. as a 
practical working limit for spectra at 86A/mm _ with the Radcliffe reflector, 
Shapley’s luminosity function for the Large Cloud suggests that there are between 
2000 and 3000 members within range of this equipment 

The Henry Draper Extension covering the region of the Large Cloud contains 
some 200 objects classified as O, B or Con, mostly brighter than 12-0 mag. From 
a sample of 25 such objects observed spectroscopically with the Radcliffe reflector 
(1) we can be sure that the great majority, if not all, are members of the Cloud. 
The HDE “ B” or ‘*‘Con”’ stars prove to be mostly late B or early A types with 
very narrow hydrogen lines. ‘The presence of a true absorption O star still 
remains to be established. It is probable that the HDE ‘‘O”’ objects (some of 
which are as faint as 15 mag.) are mostly Wolf-Rayet type. 

The HDE Catalogue also contains over 3000 stars classified as A or later. 
Viost of these are foreground stars. But it is of the utmost importance to dis 
entangle the minority of true members of type later than A. One such star was 
included in the previous tabulation (1), and as a result of further investigations 
the following bright stars have been found to be Cloud members from velocity 
shifts and supergiant characteristics. 


Red Super-supergiants in the Large Cloud 
HDE 
Star Sp MK Type SP¢ , ; Remarks 
68757 Mo Gs: la 11°83 3 , near NGC 1743 
269723 Ks5 Go la 10°8 near NGC 2014 
269953 Ko 2 F8 la 7 2 near NGC 2085 
271182 Ko FS la 
30 Dor ... F8-Go la 
No. 50 
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‘The magnitudes and colours in the fifth to seventh columns have been deter- 
mined photoelectrically by Dr A. J. Wesselink on the Cape system using the 
Radcliffe 74-inch reflector and Cassegrain photometer. ‘These results are to be 
regarded as provisional and in any case, in view of the large discrepancies from 
the HDE (photographic) magnitudes given in the third column, it appears likely 
that these objects of exceptionally high luminosity are variable in light. 

‘The absolute visual magnitudes listed in the eighth column are based on an 
assumed distance modulus of 19:2 for the Large Cloud (2). It seems to be un- 
desirable to apply the term “‘supergiant”’ to objects covering a range of some 7 
magnitudes, and we would suggest the term “‘ super-supergiant ”’ to be appropriate 
to objects with My brighter than — 7-0. 

With the exception of 271 182, which is in a loose grouping of stars containing 
some known Cloud members, all the above stars are near nebulae and clusters. 
30 Dor No. 50 is not in the HDE and its precise position will be defined in a later 
paper devoted to the 30 Dor complex. Its magnitude has not been measured 
but it is estimated as between 10™-5 and 11™-opg so that its absolute magnitude 
is of the same order as that of the other stars. 

Spectral classification.The MK types in the fourth column have been 
obtained from Radcliffe spectra with dispersions of 86 and 49 A/mm at Hy. ‘The 
MK standards 6 C Ma (F8 la) and HR 2974 (Go Ia) match well the last four 
stars, except that in 30 Dor No. 50 the H lines cannot be used owing to the super- 
position of strong nebular H emission ; also in this star ail lines appear to be some- 
what weakened, perhaps due to their being partly filled in by a nebular continuum. 

Judging by the strength of the G band and the weakness of H, 268757 is 
clearly somewhat later in type than HR 2974 (Go Ia), but certainly not as late as 
M (HDE). It is understood that the star has been classified as K from an ADH 
objective prism spectrum (3). Accurate classification is difficult on account of 
the fact that there are no MK Ia standards between Go and M1. No good 
match is found with Go Ib, Gs Ib, K2 Ib or K3 lab. A fairly good match is 
provided by « Aqr (G2 Ib), but Sr 11 4215 and some other lines are stronger in 
268757. ‘lhe classification as Gs: la is .ubject to considerable doubt. 

Radcliffe spectra of long period cepheids in the Clouds taken at 86 A/mm have 
been classified at F'7 la to G21(4). At this dispersion, the super-supergiants 
reported in this paper show no ceitain differences spectroscopically from the 
cepheids which are about 2 magnitudes fainter. We must conclude that the MK 
criteria are insensitive to luminosity among the very brightest stars. 

Colours.— Detailed discussion of the colours and magnitudes of these stars 1s 
postponed to a later communication when it may be hoped that further instances 
of such red super-supergiants will have been discovered. But it will be noticed 
that the colours, for the first three stars at least, do appear to be redder than for 
normal F to G giants, although they accord well with the colours of the cepheids 
(4). ‘This fact probably accounts for the HDE classifications as K or M. 

\ comparison of the colours and spectral types suggests a rather abrupt 
increase in the colour index setting in at about F8 to Go. A similar effect is also 


shown by the cepheids as will be seen by comparing the colour-period (§) and 
spectrum-period (4) relations. 


\s is well known, G to K dwarfs are bluer than giants of the same spectral 
type, the same degree of ionization being produced by the lower pressures and 
lower surface temperatures of giant atmospheres. In the super-supergiants of 
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this paper one may expect even lower surface temperatures to be required to 
reproduce a Go spectral type. 

Owing to doubts concerning the temperatures and bolometric corrections 
appropriate to these objects calculations of radii are most unreliable, but for 
268757 a radius of the same order as that of the red component of VV Cep 
(2400 ©) 1s suggested. 

The possibility of interstellar reddening has to be borne in mind particularly 
for stars involved in nebulosity, but the absorption cannot be put high without 
having to assign unacceptably high intrinsic luminosities. Moreover, in the cas¢ 
of 271 182 near the northern border of the Cloud, there is no associated nebulosity 
and no appreciable absorption can be expected in this region. 

Frequency of red super-supergiants in the Large Cloud.—Mrs Nail and Shapley 
(6) have found that star counts in the Large Cloud indicate the presence of red 
supergiants brighter than 14 mag. ‘lhe five member stars reported here indicate 
that even in the brightest range of magnitudes (M g and brighter) red super- 
supergiants can exist in appreciable numbers. Among stars of apparent magni- 
tude 10-4 or brighter, there are some 18 known blue Cloud members (excluding 
a few cases of clusters formerly regarded as individual stars) and at least 4 red 
Cloud members with types F8 or later. We can be reasonably sure that the 
census of the brightest blue members is complete, but we are very far from being 
able to claim the same for red members; in fact the probability of finding more 
instances seems to be fairly high. We can conclude that among the 30 brightest 
members of the Large Cloud not less than 20 per cent are likely to be red F8 and 
later types, few are earlier than B 5 and none earlier than Bo. ‘The conclusion is 
rather surprising since the brightest Population I stars are normally regarded as 
blue. 

A detailed list of proved members and foreground stars will be published at a 
later stage. But as a guide to the extent of the search for bright members of types 
later than B we give below the frequency distribution of HDE magnitudes of stars 
that to date have proved to be foreground objects 


Number of proved foreground stars (LMC) 

Range of 8:0 8°5 re) 9°5 10°O 10°5 11° 
m. 5°4 5°9 9°4 9°9 10°4 10°9 11’4 
Number , 6 7 4 12 ; 
The investigated stars brighter than 10-0 mag. included many type A star 
from the HD catalogue because it was hoped that objects like HD 33579(1) 
might be found. But it now seems reasonably certain that any such object would 
have been marked out by the peculiarity of exceptionally narrow lines and that 


no Cases have been missed 


Rad liffe Ohbsert atory, 
Pretoria 
1956 November 30 


Note added 1957 February 8.—A star (HDE 269810, 12-1, Sp. Go) near 


the nebula NGC 2032 has been proved to be the first instance of a true absorp 
tion O member of the Cloud. It is more than 1 mag fainter than the limit for 
inclusion among the 30 brightest members. 
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